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The production of uniform shear flow in a wind tunnel 


By P. R. OWEN and H. K. ZIENKIEWICZ 
Mechanics of Fluids Department, University of Manchester 


(Received 3 May 1957) 


SUMMARY 


A nearly uniform shear flow was obtained in the working 
section of a wind tunnel by inserting a grid of parallel rods with 
varying spacing. 

The function of such a grid is to impose a resistance to the 
flow, so graded across the working section as to produce a linear 
variation in the total pressure at large distances downstream 
without introducing an appreciable gradient in static pressure near 
the grid. A method of calculating a suitable arrangement of the 
rods is described. Although this method is strictly applicable 
only to weakly sheared flows, an experiment made with a grid 
designed for a shear parameter as large as 0-45 gave results in 
close agreement with the theory. ‘here was no evidence from 
the experiment of any large-scale secondary flow accompanying 
the shear—a danger inherent in an empirical attempt to grade the 
resistance of the grid—nor was any tendency observed for the 
shear to decay with increasing distance from the grid. 


1. INTRODUCTION 


A simple method of obtaining a two-dimensional shear flow u(y) with 
vorticity w. in a wind tunnel is to insert near the beginning of the working 
section a grid of parallel rods in which the distance between adjacent rods 
varies in the y-direction, as sketched in figure 1. ‘The effect of such an 
arrangement is clearly to produce a grading in the resistance of the grid 
which in turn gives rise to a non-uniform (irrotational) flow upstream of 
the grid and a non-uniform distribution of total pressure downstream of it. 
If it were possible to select a resistance grading such that the gradient of 
total pressure far downstream of the grid were constant, the vorticity » 
would also be constant. 

In principle it should not be difficult to devise empirically a grid which 
satisfies approximately the above condition on total pressure gradient 
(a remark on the feasibility of realizing it in practice is made in §5), but 
such an approach carries with it the danger of introducing secondary 
vorticity components w,, w, since the variation of total pressure near the 
downstream face of the grid might include large static pressure gradients. 

The method of grid design described in this paper has as its purpose 
the calculation of a resistance grading corresponding to a constant total 
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pressure gradient far downstream together with a predictable and, if possible, 
small static pressure gradient near the grid, and is applicable to weakly- 
sheared flows in the sense that the maximum change in velocity across the 
wind tunnel is small compared with the mean velocity. 

On the basis of the method, a particular grid was designed and tested 
in a wind tunnel with a working section 20 in. square at a mean speed, 
near the tunnel axis, of approximately 75 ft. sec"'. Static pressure and 
Pitot pressure traverses across the working section at a number of stations 
downstream of the grid revealed, over a large part of the area, an approxi- 
mately linear velocity distribution, u(y), which was almost invariant with 
both x and zs. Furthermore, the static pressure was found to be sensibly 
constant over planes parallel to the grid at distances from it of one tunnel 
width and greater; the region in the immediate vicinity of the grid could 
not be explored with a conventional static probe for reasons explained 


in §5, 
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Figure 1. Arrangement of the grid and coordinate system. 


One of the objects of producing a shear flow on a large scale is to examine 
the flow about bodies in a stream of uniform vorticity, and thereby to 
undertake the experimental counterpart of theoretical investigations made 
principally by Hawthorne & Martin (1955) and Lighthill (1956). It is 
hoped to report later on the results of these experiments. 


2. METHOD OF GRID DESIGN 
The wind tunnel may be treated as a long channel with walls y = 0, 
y =A, in which a grid is placed in the plane x= (0. At great distances 
upstream the velocity is uniform and has the value U; far downstream 
the velocity is again parallel to the walls and its magnitude is given by 
u= U+A(y— 3h). (1) 
We suppose the fluid is inviscid and that Ah/U is so small that the 


departure of any streamline from a straight line is correspondingly small. 
Accordingly, the stream function may be written 


us = Uy i us’, 7% 0, 


= Uy + 3X(y?—hy) +", x > 9, (2) 
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where x’, which is everywhere small compared with Uh and tends to zero 
as x—> + o, satisfies V*s’ = 0. Appropriate solutions satisfying the 
boundary conditions 





Cy'(x,0) — dyp'(x,h) _ 


7 = 0 
Ox Ox 
for all x are 
yp ~ 
— = ¥ A, eh sin(nzy/h), ex 
we ; 
Las 
= > B,e*™ sin(nzy/h), x > 0. (3) 


n=1 


The presence of the grid requires that (i) dy%/dy is continuous through 
x = 0, (ii) the velocity component 0:5/dx obeys a certain refraction condition 
across x = (0), and (iii) the difference between the total pressures on a given 
streamline for large positive and negative values of x must be equal to the 
resistance per unit area imposed by the grid. ‘To these may be added the 
condition that the static pressures far upstream and far downstream of the 
grid are independent of y. 

Condition (i) leads to 


aa aa ) l Al 
az) nA, cos(nzy/h)—7 S nB, cos(nzy/h) = (; - 3) a 
ai ? n= h = l 
tle Sy(y—h) = h? > C, sin(nzv/h), 
n=l 
A, —B, =C,,Ah[U. (4) 


The refraction condition (ii) requires the y-component of the velocity 
to change by a factor x in passing through the grid. « for a grid of uniform 
resistance is known to be 1-1(1+ A)-"?, where K is the resistance coefficient 
(Taylor & Batchelor 1949), and it is unlikely that this value will be seriously 
in error for a grid with slightly graded resistance provided that K is given 
its local value. It follows that, if the resistance grading is represented by 








rh 
K = K,{1+ (y)], such that € is small and | ¢« dy = 0, 
1-1 | Kye ] (<2) (5) 
e = ar —---- €"} ~ 
(1 Ky) 7 2(1 T Ko) 

With x given by (5), the refraction condition is 

> 0B, sin(nzy/h) + (a+ be) > nA, sin(nzy/h) = 0, 

n-1 n= 


which, noting that A, and B, are O(Ah/U), reduces to 
B,+aA, = 0+ O(Ahe/U), (6) 
where a = 1-1(1+ K,)-!?. 
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Finally, we have to consider the change in total pressure along a 
streamline. According to the assumption of an inviscid fluid, the total 
pressure remains constant along streamlines upstream and downstream 
of the grid but, at the grid itself, there is a decrease in total pressure 
amounting to K}p[u(0)]?. Hence, 

py + SpU2— py — 3p[U +N yy — MANE = K4p[u(O)P, (7) 
where py, and p, are respectively the static pressures far upstream and far 
downstream of the grid. The streamline which has the ordinate y, for 
large 1 meets the grid at y = yp, and the relation between y, and y, follows 
from (2); 

V1 = Yo + U'(0, ¥9)/U + O(AU'/U?). (8) 
Equation (8) implies that, to the first order in AX/U and wv’ Uh, the 
v-component of the velocity at the grid is given by 


u(Q) —-U + A( V9 — th) + or (Y) 


which, together with (7), leads to 


Al , 1 2Ky cw 0 0 
PoP - K{l+e)+2—(14 Ky) ;) See (10) 


LoU/* = l h r cy 
Since (P)—p,) is independent of v, (10) resolves into 
OED oe (11) 
spl . 
and : , Ah ho 1] x Ky cds’ (0, Vo) + 
Aye+27 (1 KF = 5) oe ay (), (id) 


The distribution e(y,) follows directly from (12) if we note from (3), (4) 
and (7) that 


eys'(0, Yo) > a < 


a y ] 
a ait-—— 8 aC codeeyii) = — = (3 0 - (13 
ev | (l+a),< : i i dia "its ra)\h 5] oe 


hence, ro ped l % 4 
clang) 1 rra(% 5): _ 


This result could have been arrived at without introducing the Fourier 
series, but the series representation (3) is useful for calculating streamlines 


and pressures away from the grid. 


PRESSURES ON THE FACES OF THE GRID 


Equations (9) and (13) demonstrate that, of the ultimate perturbation 
\(yv— 4h) to the velocity in the channel, a fraction 1/(1+ a) appears at the 
grid and is produced irrotationally in the flow upstream. Clearly, as a 
approaches zero, a condition which could be realized by making Ay very 
large, a progressively greater proportion of the velocity adjustment is 
achieved before the flow reaches the grid; consequently the static pressure 
on the downstream face of the grid must tend to a constant value consistent 
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with a uniform distribution of vorticity w.. (The static pressure on the 
upstream face will, of course, not be constant and we shall refer to its etfect 
presently). Whilst the implication of this argument is clear enough as a 
criterion for avoiding the development of secondary circulation in the flow 
downstream of the grid, it is not very helpful in practice because there is 
always a limit to the value of the resistance beyond which a useful per- 
formance cannot be obtained from the wind tunnel in which the grid is 


placed. 
To explore the possibility of a compromise between constancy of static 
pressure in the plane v = +0 and a tolerably small value of Ay—approxi- 


mately 1, say—note that the pressure on the downstream face of the 
grid is given by 


C.(+0,y,) = 2-2 a (J) a (; is 3) r (tz): 
1oU? U/(l+a)\hk 2 U? 
(15) 
Although (15) shows that the pressure coeflicient departs from a constant 
value only by O(AA/U) when Ay, and therefore a, are O(1), it cannot be 
inferred immediately that such a pressure variation is acceptable unless the 








associated secondary vorticity, w,, is of a still smaller order of magnitude ; 
for the primary (non-dimensional) vorticity w,h/U induced by the grid is 
itself only O(Ah/U). 

The development of secondary vorticity arises from the action of the 
pressure gradient op/cy on the fluid in the boundary layers on the tunnel 
walls ats = +h. Such fluid, on account of its reduced velocity in the 
y-direction, more readily acquires a motion in the y-direction in 
response to the pressure gradient than the faster-moving fluid in the main 
stream. Once a cross-flow is well-established in the boundary layers a 
compensating flow must appear in the main stream in order to preserve 
continuity, and it follows that, if the maximum cross-flow velocity in the 
boundary layers is v’, the secondary vorticity in the main stream will be 
O{v'd/h?), at most, where 6 is the boundary layer thickness. 

Now the maximum value of the pressure gradient dp/oy is, from (15), 
pUXa/(1 +a) on the downstream face of the grid, and it falls off with x like 

“zh On the assumption of laminar boundary layers on the walls of the 
tunnel (it is known from experiments with curved pipes that the secondary 
tlow in these circumstances is more powerful than when the flow is turbulent), 
a rough argument based on consideration of the balance between the flux 
of boundary-layer momentum in the y-direction and the combined forces 
due to the pressure gradient and skin friction indicates that the maximum 
value of v’ attained some small distance downstream of the plane of the grid 
is O(Ah). It follows that w, is O(AS/h) at most and, owing to the presence 
of the small factor 6h, is of an order of magnitude smaller than w.. Since 
the boundary layers downstream of a grid are certain to be turbulent in 
any apparatus of useful dimensions, the mean resistance may safely be 
confined to values near 1. 








526 P. R. Owen and H. K. Zienkiewicz 


Turning our attention to the flow upstream of the grid, we notice that 
any attempt to ameliorate the pressure gradient at x > 0 by increasing A, 
will result in an increased pressure gradient at x < 0, for the maximum 
values of op/éy, which occur on the faces of the grid, are related by 


~ : ee, 
\ oy t 0 a ey z +0. 


and so the danger of inducing a secondary circulation is merely transferred 
from x > 0 to the region x < 0. We have however already rejected the 
use of large Ay on the grounds of wind tunnel performance, and it remains 
to consider the suitability of values of A, near 1. The problem can be 
disposed of quickly by appealing to the arguments of the preceding 
paragraph. ‘The conclusions reached there hold equally well for the region 
x <0, with the reservation that the increased severity of the transverse 
pressure gradient in the flow approaching the grid (a < 1 for Ky = 1) 
will be offset by the smaller boundary layer thickness on the vertical walls 
and by the attenuating action of the grid on any secondary vorticity that 
happened to develop. 





4, RESISTANCE GRADING 

At large Reynolds numbers the resistance of a grid of rods in a uniform 
stream can be predicted with good accuracy if the drag coefficient based 
on the blocked area and the average interstitial velocity is taken to be 1 
(Wieghardt 1953). For our present purpose it will be assumed permissible 
to adopt the same procedure for each portion of the grid, with the possible 
exception of places where the spacing changes rapidly with y, since the 
scale of variation in the y-direction of the velocity through the plane of 
the grid, u(0), is much larger than the greatest spacing between adjacent 
rods. Accordingly, Ky(1+«) = &(1—€) 2, where € = d/s; dis the diameter 
of each rod and the spacing is s(y). With e« given by (14), 


é rr ee 
oo, é a) Al) a aoe _ eae 
(i- “Fi = - ky 2a (x ra) 3) } (6) 


It will be observed from (16) that the choice of Ay is not entirely arbitrary, 
since for each value of AK, there is a maximum value of Ah/U which can be 
reached with a physically-realizable grid. ‘This maximum, which is given 
in the following table, is governed by the condition that € > 0 at y=A. 





] 

| 

| 0-8 
| 
c | 

| 





| | 
| 0-2 10-4 | 0-6 
| | 
| | 


0:56 | 0-64 | 


(MW/U) max | 0-18 | | 0. 
| 





5. EXPERIMENT 
A grid composed of parallel rods of 0-125 in. diameter was constructed 
according to (16) with Ahk/U = 0-45 and A, = 1:15. A conservative choice 
of Ky, as compared with the values shown in the table above, was considered 
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desirable in order to avoid unduly rapid rates of change of the rod spacing 
with y near y/h= 1. The appropriate values of € are shown in figure 2, 
which also contains the values measured on the centre-line, s = 0, of the 
grid after it was built. 
























































OSI 0 Py 
A “eo. 
= a o 
ee 
# 
= 
oO” 
el 
| 
"O 02 4 6 8 4k 1 
Figure 2. Spacing of the rods; ———— design values of €, -- 0 --0-—-— measured 
values on the grid centre-line. 
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Figure 3. Calculated streamline pattern for Ah/U = 0-45 and K, = 1:15. 


The grid was placed at the entrance to the working section, 20 in. 
square, of the low-turbulence wind tunnel at the Mechanics of Fluids 
Laboratory. ‘l'raverses across the working section in planes parallel to 
the grid at a number of distances from it were made with Pitot and static 
tubes of 1mm outside diameter at a mean windspeed U of roughly 
1 


75 ft. sec 
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No variation in static pressure could be detected from the traverses, 
which were confined to distances greater than approximately one tunnel 
height from the grid. ‘The grid could not be approached much mor« 
closely with a conventional static probe owing to an appreciable inclination 
and curvature of the streamlines, the predicted pattern of which is shown 
in figure 3. 
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Figure 4. Velocity contours in the plane xh — 3; - measured, 


- - calculated. 


‘The contours of uw U as deduced from the Pitot and static pressure 
traverses in a plane parallel to the grid and 5 ft. (v/h = 3) downstream ot 
it are shown in figure 4+. ‘The absence of any systematic departure from 
lines parallel to the s-axis, together with the previously-mentioned fact 
that the static pressure is constant over the region, suggests that no secondary 
flow was present on a large scale. 

The profiles of u/U in the plane of symmetry of the tunnel, y = 0, 
at distances of 2, 5 and 7 ft. from the grid (xh = 1-2, 3-0, 4-2) are presented 
in figure 5. ‘The most striking feature of this figure is the close agreement 
between the measured velocity distribution and that predicted by the theory : 
an agreement which is remarkable in view of the fact that the grid was 
designed for a value of Ah, U of 0-45, which could hardly be accepted as 
small enough to come within the scope of a linear theory. ‘To be sure, the 
presumed effects of neglecting squares and higher powers of AW/U can be 
interpreted in the experimental velocity profile as a deficit in the average rate 


of shear together with a small departure from linearity, although the latter 
effect is equally likely to be due to an inadequacy in the method of relating 
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the spacing of the rods to the local resistance. (‘The curvature of the velocity 
profile in figure 5 is most marked in the interval 0-70 < y/h < 0-95, where 
according to figure 2, € changes rapidly with y/h.) 

‘The small-scale irregularities in the velocity profile, which were found 
to be repeatable in the experiments, can be attributed to inaccuracies in 
the manufacture of the grid, as displayed in figure 2, and a detailed 
examination of the measurements (not all of which, for the sake of clarity, 
are reproduced in figure 5) suggested that the small amplitude waviness 
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Figure 5. Velocity profiles measured in the plane of symmetry of the tunnel. 


of the velocity profile in the region y/h < 0-5 could be associated with the 
corresponding abrupt fluctuations in €. In particular, the pronounced 
dip in the velocity profile near y/A = 0-5 might be caused by the error in € 
at roughly the same value of yh; for it is easily shown that the change 
in u/U brought about by a small change in € in the interval 


y—sOy <y <yt+ddy 


- 1 (,, Ko\(_1+2Ky _,\Au - 
Janet. ole ) (i+ 4K, i 7 


0 


is given by’ 





provided that outside the interval € is altered by a constant amount so as 


} 


to preserve the condition | edy = 0. In the present case, with A, = 1-15, 
a = 0:75 and Ag = 0-015, (17) gives Au/U = —0-03 which is near to the 


This relation offers the possibility of an iterative method of designing the grid, 
by adjusting € to compensate for the differences between the required velocity 
profile and that observed in an initial experiment. Such complication was not 
considered necessary for the case considered here, although in choosing a design 
value of Ah/U as large as 0-45 we were anticipating an iterative approach. 
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observed amplitude of the dip (Au/U = -—0-02). Closer agreement could 
not be expected because the error Ag, being caused by a sag in some of the 
rods rather than a bodily displacement, is not independent of 3, and th« 
resulting distortion of the flow from that of nearly-uniform shear is essentially 
three-dimensional and must resemble the effect of a wake embedded in the 
stream. Indeed, a comparison between the observed velocity profiles at 














s/h = 0 and z/h = —O0-2 displayed in figure 6 for x/h = 3 shows that the 
dip had almost disappeared at 2/h = —0-2 (4 in. from the central plane 
of the tunnel). 
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Figure 6. Velocity profiles for z = 0 and 2/h 0-2 measured at x/h = 3. 


The above calculation provides a stricture on the remark made in the 
Introduction that it might not be difficult (in principle) to devise empirically 
the correct spacing between the rods, because the misalignment that gave 
rise to the error 0-015 in € amounted only to 0-02 in. ‘The chance of adjusting 
by experiment the positions of 62 rods, as used here, to an accuracy 
significantly greater than ()-02 in. seems very remote! 

The velocity profiles of figure 5 exhibit complete departure from linearity 
near the walls, a behaviour evidently to be associated with the presence 
of the boundary layers. In this respect, two points require comment. 
In the first place, near each wall the velocity reaches a distinct peak which 
tends to subside with increasing distance from the grid. Secondly, the 
peak near y=A occurs closer to the wall than that near y=0. ‘The 
appearance of these velocity peaks is due to the accelerating action of the 
pressure drop across the grid on the fluid in the boundary layer which, by 
virtue of its reduced velocity, suffers a smaller change in total pressure on 
passing through the grid than fluid in the main stream. Since the change 
in total pressure diminishes progressively as the wall is approached, whereas 
the static pressure is constant, the velocity in the boundary layer would, 
in the absence of shearing stresses, increase towards the wall. This 
tendency is however resisted by the shearing stresses and wall friction, 


with the result that a maximum velocity occurs in the outer part of the 
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boundary layer. With increasing distance from the grid, the shearing 
stresses smooth out the velocity peak, as can be seen in figure 5. The 
difference between the positions of the velocity maxima on the top and 
bottom walls is obviously a consequence of the difference in boundary 
layer thickness; for the boundary layer on the upper wall is subjected to 
a favourable pressure gradient whereas the pressure gradient on the lower 
wall is adverse, as is evident from the calculated wall pressure distribution 
shown in figure 7. The observed difference in boundary layer thickness 
on the upper and lower walls agreed quite well with a rough calculation 
based on 'Truckenbrodt’s method (Schlichting 1955) using the pressure 
distributions of figure 7. 
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Figure 7. Calculated pressure distributions on the top and bottom walls of the tunnel; 
Ah/U = 0-45, K, = 1°15. 








A final remark may be made on figure 5. Ostensibly there is no 
systematic change in the rate of shear between x/h = 1-2 and x/h = 4-2, 
and even the small irregularities in the velocity profile appear to suffer 
no appreciable decay. ‘This is not entirely surprising, because a uniform 
shear flow would be expected to persist almost indefinitely if, as is likely, 
the turbulence generated at the grid gives rise to an approximately constant 
eddy viscosity. All modification to the rate of shear would then be confined 
to regions excluded from the action of a uniform shearing stress, that is, 
to the boundary layers on the tunnel walls, which only gradually encroach 
on the main stream. 
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An exact solution for progressive capillary waves of 
arbitrary amplitude 


By G. D. CRAPPER 
Department of Mathematics, University of Manchestei 


(Received 4 March 1957) 


SUMMARY 

An exact solution is found in a fairly simple form for two- 
dimensional progressive waves of arbitrary amplitude on a fluid 
of unlimited depth, when only surface tension and not gravity 
is taken into account as the restoring force. ‘The calculated wave 
forms are exhibited graphically for various amplitudes, and the 
relation between wave velocity and amplitude is plotted. ‘The 
wave of greatest height occurs when the vertical distance between 
trough and crest is 0-730 wavelengths (compared with 0-142 for 
gravity waves). Higher waves are prevented from appearing by 
the enclosing of air bubbles in the troughs. 


1. INTRODUCTION 

Periodic gravity waves of finite amplitude have been investigated by 
Stokes (1847; 1880, pp. 197 & 314; see also Lamb 1932, § 250) and others, 
and a good approximate solution has been found. ‘The gravity wave of 
greatest height has also been determined (Stokes 1880, p. 225). It is 
therefore interesting to see if analogous solutions can be obtained for purely 
capillary waves, for which gravity is neglected; it is particularly interesting 
as a first step towards investigating the problem of waves of finite amplitude 
under the combined effects of gravity and surface tension. 

Capillary waves are found to be remarkable in that an exact solution 
exists for arbitrary amplitude. ‘The case where the fluid has infinite depth 
was considered first, and, following Stokes, a series approach was made. 
Inspection revealed regularities in the coefficients, however, and the series 
obtained by assuming that these regularities continued in the higher terms 
was found to have a sum in closed form. This sum was then shown to be 
an exact solution. From this original investigation the present method, 
leading directly to the solution, was evolved. ‘The present method indicates 
that there is also an exact solution if the fluid has finite uniform depth. 
‘The analysis is, however, rather complicated, involving elliptic functions, 
and the solution 1s not considered worth evaluating in detail. 

The solution for infinite depth shows that the wave velocity given by 
the theory of waves of infinitesimal amplitude (Lamb 1932, $265 & § 266) 
is accurate only for zero amplitude and that, if the wavelength is fixed, the 


velocity decreases as the amplitude increases. ‘This behaviour is exactly 
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opposite to that of gravity waves. The waves themselves are very rounded 
in section, and the wave of greatest height is reached when the surface 
bends back to touch itself, enclosing a bubble of air. 


2. EQUATIONS OF MOTION IN THE (4, y)-PLANE 

Our investigation will be confined to two-dimensional travelling waves 
on the surface of an ideal fluid which has infinite depth. It is convenient 
to choose Cartesian axes with x measured horizontally to the left and y 
vertically downwards. ‘To make the flow steady we bring the waves to rest 
by superimposing a uniform velocity on the system, and therefore we 
assume that the undisturbed fluid (that is, the fluid at great depths) is 
moving in the positive x-direction with velocity c, the wave velocity. 

If the motion is generated from rest, the flow must be irrotational, and 
so the equation of motion is 








= 0, (1) 


where ¢ is the velocity potential of the flow. The stream function % also 
satisfies this equation. 

In a steady flow the surface is a streamline, say % = 0, and Bernoulli’s 
equation holds, so we have 

pip +3q = 3c. (2) 

Here p represents the difference of pressure from its hydrostatic value ; 
the constant on the right is }c? because at great depths, where q = c, this 
difference is zero. Surface tension creates a pressure difference across the 
surface which is given by 


P-Po= T/R, (3) 
where p is the pressure in the fluid at the surface, py is the ‘atmospheric’ 
pressure, 7 is the surface tension, and R is the radius of curvature of the 
surface, counted positive when the centre of curvature lies inside the fluid. 
The use of p to represent the absoluie pressure at the surface as well as its 
ditference from the hydrostatic value is permissible since gravity is neglected 
as a restoring force. Now when R is positive, d*y/dx* is positive, and so 

] d*y/dx* 





iS 4 
R~ tit (dyjdee (4) 
The boundary condition is therefore 
T d*y/dx : 7 
—- ——~___, +ig=ice ond =0. (5) 
p {1+(dy/dx)?P*  * 
As the fluid is of infinite depth there are the further conditions 
ch ous 
—— = = P é. 
Ox oy (6) 
as y > 00. 
ch Cds te 
— sont seit: SF 0, 
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The problem as formulated is very complicated, but it becomes much 
simpler on taking new independent variables (¢,%) and new dependent 
variables (7,4), where + = log q and (qcos 4, qsin @) are the (x, y)-components 
of the velocity of the fluid. 

On any streamline 


—_=—, (7) 


where ds is an elementary arc of the streamline, so 
l (=) (3 
ip — tion ’ 5) 
R ag y= const 


9 


kg®=ic? ond = 0. (9) 


and we can write (5) as 


‘This can be simplified further if we take 7/pc? as unit of length and c as 
unit of velocity, and this gives 
eo t 
e ad tse" =} ony = 0. (10) 
Thus, the boundary condition in the (¢, %)-plane is 
06 


== —sinhr ong =0. tt) 
0d ! 
If w= 6+ and z = ¥+1y, 
dws 
_ * sie 
and loe( 2) — ee | (13) 
dz 


These must be regular functions, so we have, equivalent to the original 
equation of motion (1), 

A) A) 

ovr O°T 
and for the same reason we can use the Cauchy—Riemann equations to 
write (11) as 


(14) 


OT 


— = —sinhr onp=0. 15 

om) Ss ub (15) 

In the new units g -> 1 as y + ~%, so that the conditions (6) become 
7>Q and@>0 aso o~. (16) 


3. GENERAL SOLUTION FOR 7 
‘lo solve this non-linear problem we look for a solution of (14) which 
satisfies the equation 


. = —f()sinh'z (17) 


Oo 
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i 
i 








An exact solution for progressive capillary waves 535 





everywhere in the flow, for some function f(s). Such a solution will satisfy 
the boundary condition (15) provided that 


f(0) = 1. (18) 

We can integrate (17) directly: 
logtanh $+ = F (ys) + G(¢), (19) 
W here na F (ib) f (us) (20) 


and G(¢) is arbitrary, or 


(X(b) + ¥(d)) 


GM aaa ——eeeee 2 
lo8) XGh) ¥(4){ oF) 
where X(s) = e-F), (22) 
and Y(4) = e@), (23) 


Functions X(ys) and Y(¢) can now be determined so that Laplace’s 
equation (14) is satisfied. For this to be so we must have 





y” X"(ys 
25(X"(ah))2 + (Y'(b))2} + {X2(b) — ¥2()} { ia, = vat =(). (24) 


Differentiating this with respect to ¢ and then with respect to y gives 


wy (Fy) + (Gay) wy (25) 


or mana tie jin 
(.X°(ys)) I (a) + (¥?(d)) 07 0. (26) 


Here each term must be merely a constant, say 





' , X"()\' 
| 2(as ees 27 
(.Y°(s)) nas) (27) 
and 7 ; ¥"(d) 
24) = —Al- rs 
(¥?(d)) + y Tar) (28) 
Integrating, we have ts X" (ys) 
g g, 214), AW. 2 
Xs) =h Yh) (29) 
™ Y"(¢) : 
0 ’ +m, 3 
Y*() h (4) m (30) 
where /, m are arbitrary constants, and (24) reduces to 
7 , a NX" (ys) )) fe d) 
) "(ys))" ifs ys) 4 ———me =f ~) os 79) t ————— 
2(.N" (ub)? — N(b) NX (b) + m XUh) 2(Y'(4))?— Y(*) (4) ! Tg) 0. 
(31) 
cae X" f 
iia 2(.N'(ys))? — Xs) X"(ab) + m- ik: ) n, (32 
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Y"(¢) 





"YV'(db))2 — V(b) V(b) 4 ~_ 33 
2(Y'(¢))? — Y(4)Y"(d) +1 V4) n, (33) 
where, again, m 1s a constant. 
: 29 = X(u), .-, 7 
From (29) X"(p) = MO) x4) —0), (34) 
rom (32 as ‘ 
cee ee 5 (Xb) — I(X%h) — m). (35) 
But differentiating (35) gives 
X (uh , ‘ 
A “(ys) : Care 2(ys) —/+X *(us) m), (36) 
and therefore l=m, (37) 
panies Je : 
2( \ (ys))* nN -r Ace 2(ys) - /) (35) 
simil: , ea. 
iain, 2 ¥'(p))? = —m— 7 (V%$)— De. (39) 
Thus the general solutions for X(s) and Y(q) satisfy 
(X'(b))? = a, +a, X*(ys) + a, X*(eb), (40) 
(Y'(¢))? = —a,—a, Y*(¢)—a, Y(), (41) 


where @,, @, dy, are any constants. 

‘The quadratic case (a, = 0) has been found to satisfy our requirements 
for fluid of infinite depth, but inspection of the full quartic has shown that 
it will give a solution when the fluid depth is finite. It has already been 
pointed out, however, that the analysis needed is too elaborate to make this 
solution worthwhile. 


4. DETAILS OF SOLUTION 
If we put a, = 0 in (40) and (41) we have 
X (ib) = (a,/ao)"?sinh( + pal" +C), (42) 
¥() = 1(a,/a,)'*sin( + das/* + D), (43) 
where C, D are constants of integration. Clearly {| ¥(d)/X(uys)} + 0 as 
ys» x, so that the condition 7 — 0 (16) is satisfied. 
From (20) and (22), 


A 


f(b) = = log X(W) (44) 
= +kcoth(+k’b+C), (45) 
writing k for ai”, and hence the boundary condition f(0) = 1 determines C: 
& + Slog(+7/A)?, (46) 

there wa] 
where I - (47) 


k+1° 
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Now, by (21), (42) and (43), 





( sinh( + ks + C)+ sinh( + kd +7D) 
7 = logy — - =. : (48) 
| sinh( + kb + C)—sinh( + tk6 + 7D) 
tanh P 49 
og ———. , 
~ tanh O 7) 
where P=3(+kb+C+ikh+iD), | 
onl | (50) 
O = 1( tkb+CFikdb — iD). 
Hence, OT sech?P __. sech?O =: : 
— = tik — —, + | pee“ S —— oe (51) 
Cys “tanh P * tanhO = ad 
so that # = tlog(coth Pcoth Q) (52) 
and #6 + 0 as ¢ 0. Thus 
dw 
log — = logcoth?O, (53) 
E dz ~ 
‘nce lw eFikw+C-iD _ ]\2 
ind hence . nm ( satelite ) (54) 
z eee ee 
From (46), e& = +1/A_ or FiA, (55) 
and, as it only amounts to adding a constant to ¢, it is permissible to put 
D= 37. ‘Then we find that both the alternatives (55) lead to 
dz 1 F Ae*hr\2 6 
— = |—— } . 5 
dw ~ \T+Aere ar) 


Taking the alternative sign is equivalent to adding a constant z/k to w 
(i.c. to d), so the solution is unique: 

47 l 47 = 
— ——- +. 57 
k 1+ denn + F 07) 


z= w- 


The constant has been chosen to make s = w when A = 0. Returning to 
the original length and velocity units, we get 
w 41 l 47 59 
25 — — ~—__—__ + -. 5 
k 1+ Ae* k a 
If w is increased by 2z7c/k, the only effect is to increase z by 27/k, so we 
have, for the wavelength A, 


A = 2n/k 


{1-A*\ T 
pn] (EE ee nt § 
2n (j =) = (59) 


Hence . = 12/1 — A2\172 a 
z) (=z) (60) 


and 


\ 


21 1] 2 
= og a ‘a... mn Cae 1 = 
aw 1+ Ae? 


(61) 
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where x = w/c, (62) 
and the range 0 4. 1 is one wavelength. On the surface & = 0), 
x = (¢/cA) and 
\ 2 A sin 27% 
x Wes, pices 
A a | | ZA Cos Z7re 
r (6.5) 
x Zz Z 1+ A cos2z7zx | 
A nm 714+A?4+2A cos2r7x J 


From these relations it is easily seen that, if ais the amplitude of the wave, 
defined as the vertical height between trough and crest, then 





a 4A 64 
X~ m1— AB)’ i 
or 2A ( ma"\ 12 ) z 
y path 4 on 05 
A =4(1 at) If (65) 
RESULTS 


The full equation of the surface, in terms of the parameter z, is 


2 21 2 7a" 12 hic . 3 
~= @ a Poe -—, (66) 
A 7 7a J 7 


Preliminary calculations showed that for large enough values of a/A 
the wave surface crosses itself, and this: suggested that the wave of greatest 


height would occur at the value of a/A for which the surface was tangent 








to itself. 
For this to be possible we must have x = 0 for x #0; Le. 
4. I sin 2 TY. 1 
1+ 4A2+2A cos2z7a = . : (67) 
7 
or A? —2A f(27x)+1 = 0, (68) 
where . 2sint 
f(t) =- —cost (69) 
and we want the least value of a/A for which this condition is satisfied ; 
A = f(to)— {(f(to))?— 1)", (70) 


where f(f,) is the minimum value of f(t). This gives for the wave of greatest 


height a/A = 0-730. (71) 
(0-142 (Michell 


By contrast, the gravity wave of greatest height has aA 
1893). 

An interesting point to notice is that the surface y 
A, is the streamline ys = const., given by 


A,e~7" = A, (72) 


() for waves with 


amplitude corresponding to A= 


for waves with amplitude corresponding to 4 


fact simplifies the computation because the streamlines for the wave of 
greatest height, shown in figure 1, are themselves surface shapes for the 


A, (A, > A,). This 
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Figure 1. 
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Streamlines for the wave of greatest height. 
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shape for the stated value of a/A, and if any particular one is taken as surface 
those belowit are still streamlines. 
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for waves of length A. 
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stated values of (a/A), and if any particular line is taken as surface, those 
below it are still streamlines. ‘The almost circular shape of the crests in 
figure 1 is to be expected in what the streamlines indicate to be a region ot 
slow flow, since here surface tension is the dominating force. 

Finally we have, from (60) and (64), 


2a T\12 ma*\-14 me 
i (5 ) (13 4 (73) 


‘The function (1 +7%a?/4A?)"14, which shows how the wave velocity 
falls away from the value given by linear theory as the amplitude increases, 
is drawn in figure 2. 
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The flow past a rapidly rotating circular cylinder in 
a uniform stream 


By D. W. MOORE 


Jesus College, Cambridge 
(Received 5 April 1957) 


SUMMARY 

A circular cylinder is rotating in a uniform stream of viscous 
liquid, the circumferential velocity of the cylinder being much 
greater than that of the uniform stream. It is shown that the effect 
of the uniform stream can be regarded as a small perturbation of 
the rotary flow due to the cylinder, and that a uniformly valid 
first approximation to the flow field can be obtained in this way. 
‘This approximate solution has the interesting property that the 
circulation about any circular contour concentric with the cylinder 
is the same. ‘The lift on the cylinder is shown to be that given 
by the classical formula, and the drag force is zero to the order of 


the approximation considered. 


1. INTRODUCTION 
In inviscid fluid theory an important type of two-dimensional flow is 
that generated when an infinite cylinder is placed in a uniform stream, with 
its generators normal to the plane of flow. Since the space outside the 
cylinder is doubly connected the velocity potential is determined only if 


the circulation K = | u.dl is specified about some contour I" encircling 


- | 

the cylinder. ‘Vhe cylinder is then found to experience a force of magnitude 
KpU per unit length perpendicular to the direction of the uniform stream, 
where p is the fluid density, and U is the velocity of the uniform stream. 
This result is due to Kutta (1910). ‘Taylor (1925) showed that this same 
formula for the force on the cylinder is valid for a viscous fluid, provided 
that K is defined by a contour which is everywhere distant from the cylinder, 
and which crosses the cylinder’s wake normally. This result explained the 
experimental work of Bryant & Williams (1925), who had compared the 
measured value of the lift with the value obtained from the inviscid formula 
by inserting the measured value of K. 

To determine the value of K for given boundary conditions on the cylinder 
it is necessary to take the viscosity of the fluid into account. Sometimes 
this can be done without reterence to the detailed viscous flow pattern. 
In the case of a slender aerofoil terminating in a cusp, for instance, the 
well-known Kutta—Joukowski condition, that K must be chosen so that the 
velocity of the potential flow is finite at the cusp, is a very plausible version 
of the requirement that steady viscous flow should be possible in the 
boundary layers near the cusp. In general, however, it is not possible 
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to determine K without solving the Navier-Stokes equations for the flow, 
and it is with one such case that the present paper is concerned. 

The problem considered is that of a circular cylinder of radius a rotating 
with a circumferential velocity V in a uniform stream of velocity U. ‘The 
viscous stresses set up by the cylinder’s rotation induce a circulation in the 
fluid, and the cylinder experiences a transverse force in accordance with the 
general theory mentioned above. ‘lhe existence of this force was discovered 
by Magnus in 1852, and the phenomenon is now usually known as the 
‘Magnus Effect’. 

When the Reynolds number Va v is large the viscous forces may be 
expected to be small outside the boundary layer, and therefore the well- 
known inviscid theory will apply outside the boundary layer. ‘The flow 
pattern given by the inviscid theory depends critically on the ratio 27Ua/K. 
If this ratio exceeds } there are two stagnation points on the circumference 
of the cylinder. If 27la/K ;, on the other hand, there are no stagnation 
points on the cylinder, and there is a region of closed streamlines near the 
cylinder. In this case, closed boundary layers of a well-behaved type 
might be expected to form if the Reynolds number of the flow were large. 
If the circulation is produced solely by the cylinder’s rotation it is plausible 
to assume that A is O(al’), so that if UV is small closed boundary layers 
may be expected to be formed. ‘These are observed experimentally when 
U/V < { (Prandtl & Tietjens 1936, plate 7). Wood (1957) has considered 
this case from a boundary layer point of view, by expanding the velocity 
in the boundary layer in powers of UV and substituting this expansion 
into the boundary layer equations. 

In the case U/V = 0 the problem has a simple exact solution for all 
values of R, and it is reasonable to expect that when U/V < 1 a solution 
might be obtained, for all values of R, by regarding the uniform stream as 
producing a small perturbation of this exact solution. It is the purpose 
of the present paper to examine this approximation. Since the perturbation 
velocity becomes comparable with the velocity of the basic flow at a great 
distance from the cylinder, the perturbation method is obviously in danger 
of being invalid in this region of the flow. In $2 it is shown that this 
difficulty can be overcome by using two approximations, one appropriate 
to the region near the cylinder, and one appropriate to large distances from 
the cylinder. In particular a uniformly valid first approximation to the 
velocity distribution can be found for all values of the Reynolds number, 
and a uniformly valid second approximation to the velocity distribution 
is found when the Reynolds number is greater than v48. — Further 
approximations are not obtained explicitly, but their general nature can be 


inferred, and their relation to Wood’s expansion elucidated. 

The remainder of the paper is concerned mainly with the first approxima- 
tion; the nature of the flow pattern, and the forces on the cylinder, being 
considered in $3 and § 4, respectively. It appears that, to the first order 
in U/V the circulation round every circular contour concentric with the 
cylinder is the same. ‘Thus, to this order, either the circulation at the 
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cylinder’s circumference, or round a contour very distant from the cylinder 
can be inserted in the formula for the lift. ‘This is verified in § 4 by a direct 
calculation of the lift. ‘The drag force is also calculated, and is shown to be 
zero to this order in U/V. ‘lhe results provided by a second approximation 
are also considered briefly in $3 and §4. In particular, the circulation 
round a distant contour is determined, and, in the limit of infinite Reynolds 
numbers, is shown to be the same as that predicted by Wood from boundary 
layer considerations. 
2. ‘THE FIRST APPROXIMATION 
Using polar coordinates (ar,#) with the origin at the centre of the 
cylinder, and the line # = 0 in the direction of the uniform stream, the 
equation for the stream function Vay is 
Vig = R E oS ii. oe (vy) | (1) 
r cOcr r or od ‘ 
where, it will be remembered, R = Va/v is the Reynolds number based 
on the circumferential speed of the cylinder. ‘The boundary conditions 


to be satisfied by % are, 


ods 
wu = 0, ay = —], att = |Z (2) 
tf 
y~ersind asr> o, (3) 


where ¢ = U/V is the small expansion parameter. 
When ¢ = 0, the relevant exact solution of (1) is the irrotational purely 


rotary flow 


by = — log Yr. (4) 
Here, in accordance with the ideas discussed in § 1, an approximate solution 


of the form 


uy = by + ex, (5) 
is sought. ‘The function %, must satisfy the boundary conditions 
Cys, 
uy = 5 = 0 a7 = f- (6) 
yy ~rsind asr> ow. (7) 


If terms in e€? are neglected, substitution of (5) in the exact flow equation (1) 
shows that /, satisfies the equation 


Roc 19) 


eS r2 00 (V*4,) = 0. (9) 


The general solution of (8) can be obtained by assuming that it may be 
represented by a Fourier series 

t= > frlrie”; (9) 

this is suggested by the fact that 4, must be a single-valued function of 8. 


Fortunately, it is found that the boundary conditions (2) and (3) can be 
satisfied by taking only the terms forn = +1, so that convergence questions 
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do not arise. It is then easily found that 
p, = Ale’[ —1(r—1) — Fr} + Fre tio)}, (10) 
In this solution 
+(3-—p)i 
rae q ( p) (11) 


~ G+(3—py 
and p and q are real constants given by 
P, —9 = Biv/(1+ R?) + 1}? (12) 
A description of the flow pattern represented by (10) will be given in §3; 
the remainder of the present section is concerned with the validity of the 
approximations leading to this solution. 

The exact solution of (1), (2) and (3) is, of course, unknown, so that the 
validity of (10) must be considered by indirect methods. ‘Thus if (10) is 
substituted into (1), the value of the terms that are neglected when (8) is 
derived can be calculated; and the magnitude of these neglected terms 
compared with the magnitude of the retained terms determines whether 
or not the approximation is consistent. It must be emphasized that this 
procedure cannot rigorously establish the validity of the approximation ; 
for if the approximation is invalid it cannot be used to calculate the neglected 
terms, and the procedure fails. However, since the approximation is 
physically plausible, it seems reasonable to accept it as valid if the procedure 
should show it to be consistent. 

The exact equation for x, is 


R C “OC [s C : |; 
j (V?b,) = Rd me de lee 


Vid, -— 3 

r of 
where the second term on the left-hana side is the retained inertia term, 
and the right-hand side consists of the previously neglected inertia terms. 
When r = O(1) the retained inertia term is obviously dominant. At large 
values of r, however, one finds, on inserting (10) into this equation, that the 
retained inertia term is O(r~”-*) whilst the largest neglected inertia term is 
O(er”"). ‘Thus the approximation breaks down when r = O(1/e), and 
since the outer boundary condition (7) must be applied, the validity of the 
solution (10) requires further justification, even when r < IL/e. 

Now when r = O(1/e) the contributions to the velocity from both the 
basic rotary flow and the uniform stream are O(e), and it is reasonable to 
suppose that the exact solution has this same property. If, then, the 
equation (8) is solved subject to the exact boundary conditions on the cylinder 
and the order of magnitude conditions on the velocity when r = O(1/e) 
one finds that the solution must be of the form Cys, where C is a constant 
of order unity and y;, is given by (10). In order to proceed any further the 
nature of the flow in the region O(1) < er < o& must be examined. Now 
the rotational part of the velocity distribution corresponding to this solution 
is O(e”) when r = O(1/e) where, from (12), p is always greater than unity 
and increases with R. Hence, with the mild assumption that the order of 
magnitude of the vorticity does not increase with r in the outer region, the 
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flow in this region is, to a first approximation, irrotational. ‘The problem 
in the outer region is therefore that of finding an irrotational stream function 
which satisfies the boundary conditions 


y~psind asp=er—> o, (14) 


yf ~ C(loge—logp+psin@+o(p)) asp. (15) 


;quation (15) is the condition that the inner and outer solutions should 
match. Clearly a solution is possible only if C = 1, and the unique outer 
solution is then 

w= —logr+ersin#é. (16) 


Thus, subject to the assumed conditions on the orders of magnitude of the 
velocity and vorticity distributions, it has been established that (10) provides 
a first approximation to the velocity distribution throughout the whole flow 
field. It should be noted, however, that this solution provides an approxi- 
mation to the vorticity distribution only when r < O(1/e). 

It is clear that, in general, one cannot obtain a better approximation to & 
without taking account of viscous effects in the outer region. A method of 
obtaining outer expansions which included the effects of viscosity was 
considered recently by Proudman & Pearson (1957), and it is reasonable 
to suppose that their technique could be applied to this case*. However, 
when p > 2, the discussion of the rotational term given above shows that 
it is O(e*) in the outer region, so that the flow in the outer region is 


irrotational to a second approximation. In this case one seeks a second 
approximation in the inner region by writing 
b = ot ey +eby, (17) 


and neglecting terms in e? when (17) is substituted into (1). One finds 
that 4, must satisfy the equation 


I 1 1 oa 
Vha— Fag (VMs) = Rl SES )- StS 3y(PD | (18) 


"1 
iY 
The general solution of (18) can be effected by expanding y, as a complex 
Fourier series. One then imposes the exact boundary conditions 
Cys 
tf = = = 0, atr=1, (19) 


or 





and the order of magnitude requirement that the vorticity corresponding 
to this solution should be O(e*) when r = O(1/e). These conditions show 
that the most general form ys, can have is 


wb. = Ap, —f'(1)logr —f(1) + f(r) + Fe" [ar? + Br + yr’? +.g(r)]}. (20) 
In (20), f(r) and g(r) are known functions arising from the inhomogeneous 


* Filon (1928) considered some general properties of the flow at great distances 
for a cylinder when circulation is present, but did not consider the possibility 
of matching his solutions to an inner expansion. 
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part of (18), and z, 8, A, and y are unknown constants satisfying the equations 
4+B+y+g(1) = 0, 

(21) 

28 +dA+g'(1) = 0. 

The constant 6 is the value of 2+(++27R)!* whose real part is negative. 

The most general form that the expansion in the irrotational outer region 


can assume 1s 


b = log «+(—log p+p sin 6) +eA#[A(log € —log p) + Be'*/p + Ce™/p? +...) + 
+ «2.4 A'(log e —log p) + B’e'’/p + C'e**"/p” ccm (22) 


If one puts y = p/e in (20) and imposes the condition that the expansions (17) 
and (22) should match when p — 0, one finds that 


Ga) g=U, 
A=B=C-=...=0,} (23) 
1 FC): B’=C'= = () 
Thus ys, is uniquely determined, and the outer expansion is now 


wb = loge +(—logp+psin@) + ef'(1)(log « — log p). (24) 


The calculation of the functions f(r) and g(r) is straightforward, but rather 
lengthy, and details will not be given. 

The solution just obtained ts valid only when p > 2, or from (12), 
R > v.48. ‘The determination of the second approximation would, as 
has been remarked, require the consideration of a viscous outer expansion 
if this condition were not satistied, but in view of the greater practical interest 
in higher Reynolds numbers this is not attempted. 

As R is still further increased the vorticity in the outer region decreases 
in order of magnitude, and one can find progressively more terms of an 
inner expansion in powers of e. Eventually in the limit of infinite Reynolds 
number, one would arrive at the situation considered by Wood (1957) in 
which one expands the velocity in powers of € and determines the coefficients 
by substitution in the boundary layer equations. It is not clear that in the 
case of finite Reynolds numbers the inner expansion will be a power series 
in €, as it may well contain terms of the type e”(loge)”. 


3. ‘THE NATURE OF THE FLOW PATTERN 


‘The flow field surrounding the rotating cylinder presents some interesting 
features, and it will be discussed in this section on the basis of the approxima- 
tions obtained in § 2. 

‘The first step in elucidating the flow pattern is to calculate the circulation 
about a circular contour of radius r, concentric with the cylinder. ‘This is 
given by the formula 


K(r) Var | = (y+) dé, (25) 
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and on inserting the expressions given in $2 for the functions ys, and uJ, 
one finds that %, makes no contribution to the integral, so that 
K(r) = 27Va. (26) 
Thus the circulation about all circular contours concentric with the cylinder 
is the same, and is equal to the circulation at the circumference of the 
cylinder. By Stokes’s theorem, this implies that the net amount of vorticity 
contained between any two such circles is zero. It is interesting to recall 
that ‘l'aylor (1925) made the hypothesis that the net amount of vorticity 
in the wake contained between lines perpendicular to the wake was zero. 
If p > 2, the second approximation derived in § 2 can be used to calculate 
the circulation. It is found that in the inner region 


K(r) = 27Va{1+ef'(1)—-7f'(r)] (27) 
and that in the outer region 
K(r) = 27Va{1+ef'(1)]. (28) 


The function f(r) is given by 
Vert 6 +1 l/ 5 

: 20 uv ir? *"(g—1) 

= A) ee Toe - 2+z =i? | a eS 

Hr) (o+a)(o+o -2)(o+1) (o+1) (a - 1)(o- 3)(o 4 1) 

(29) 

where o = 2—p+zq, and the overbar denotes the complex conjugate. ‘Thus 
the circulation at very large distances is constant, and is given by 


—o £9 e?q(3p —5) ; 
K = 2zl| al 1 Ape ea | (30) 








It may be shown from (12) that 
Pp, -q~ (SR)? asR> ow, (31) 
and so one finds that 
K — 27Va(1—3e*) asR> &. (32) 
The result (32) was obtained by Wood (1957) from boundary layer 
considerations. 
A first approximation to the distribution of the vorticity ¢ can be 
calculated from (10): 
C= —eVa Air’ rite” F[(2 — p +iq)* — 1}. (33) 
The result assumes a simpler form when R is large compared with unity. 
It is then found that 
C ~ Ua R!?r-+ 4 cos( 32 + 6 — (4 R)"? log 7). (34) 
This expression shows that the vorticity suffers a rapid oscillatory decay as 
the wall of the cylinder is left, and that it is effectively confined to a boundary 
layer of thickness O(aR~!?). Wood (1957) has shown that this oscillatory 
character is an inherent property of closed boundary layers. 
A closer examination of (10) throws further light on the nature of the 
flow pattern. ‘The first term is simply the potential flow around the cylinder 
induced by the uniform stream. The second and third terms arise from the 
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viscosity of the fluid. At large Reynolds numbers, the third term decays 
rapidly as r increases, and corresponds to a thin, closed, boundary layer 
attached to the cylinder. ‘The second term, on the other hand, decays only 
as fast as r-', and corresponds to an irrotational secondary flow induced by 
the normal velocities at the outer edge of the boundary layer. ‘The character 
of this secondary flow is made clear by writing down the stream function 
from which it can be derived. When R is large compared with unity this 
stream function is 


ys’ ~ UarR-! *cos(6 + fz). (35) 
This secondary flow, like all such flows, is irrotational, and decreases in 
magnitude as R > » The inflow into the boundary layer is shown 
diagrammaticaily in figure 1. 
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Figure 1. The distribution of inflow velocity into the boundary laye1 


‘The mechanism producing the secondary flow is of the same general 
type as is found in other rotary flows. ‘The viscous forces modify the velocity 
of a given fluid particle, so that the centrifugal force due to the curvature 
of the streamlines no longer balances the normal pressure gradient. ‘There 
seems to be no obvious physical reason for the angular distribution of the 
flow, particularly the symmetry about the line 6 = }7. 

It appears, then, that at large Reynolds numbers the flow pattern outside 
the boundary layer is very nearly that predicted by the inviscid theory. 
More precisely, at any point outside the boundary layer, 


(viscous stream function) = (inviscid stream function) + O(R-!?). (36) 


In particular, there is no wake; a fact which will be found to have 


important consequences for the drag force experienced by the cylinder. 











The flow past a rapidly rotating cylinder 549 


The formation of a boundary layer as R + 0, with the corresponding 
occurrence of large velocity gradients, may make some of the terms neglected 
in the derivation of (8) of comparable order to those retained. Inspection 
of the exact equation (13) for ys, reveals that the terms retained are O(:be/5*) 
whilst the largest neglected terms are O(:/s?«2R/8*), where 4 is the dimension- 
less boundary layer thickness. In the boundary layer + = O(R~'?) and 
5 = O(R'!?), so that the retained terms are O(c R®*) and the largest neglected 
terms are O(e?R??). Thus the neglected terms are still smaller than those 
retained by a factor e, and the validity of the approximation is unaffected by 
the presence of a boundary layer. 


4. 'THE FORCES EXPERIENCED BY THE CYLINDER 
The solution obtained in $2 can be used to calculate the forces exerted 
on the cylinder by the fluid. 
It may be shown that on the surface of the cylinder 





a ot Yr oF 
and P,, = —p(i, 9), (38) 
where p(1, 4) satisfies the equation 
pt, @) _ _ wt ( ee ( 
“— Sea +1) : (39) 


If the solution ys, +s, found in §2 is used to calculate these stresses, 
the force (X, Y) acting on the cylinder can be determined by integration, 
and one finds that 

X = 0, (40) 
Y= = 2raV pl 2 (41) 

Thus the drag force is zero to the order of the approximation considered. 
This is perhaps a surprising result, but it follows from the character of the 
flow at large distances from the cylinder. ‘There is no wake in this flow, 
since the boundary layers are closed, and consequently the flow at great 
distances is (to this order) everywhere irrotational. It follows from 
momentum considerations (Taylor 1925) that the drag force is zero. 

The result (41) for the lifting force is also in agreement with Taylor’s 
conclusions (discussed in $1), since, as was shown in § 3, the circulation is 
2zaV to the first order. 

The above results have been derived from the first approximation 
ys) +x, which, as was shown in § 2, provides a valid approximation to the 
flow pattern for R > 0. If attention is restricted to the case p > 2, the 
second approximation to the flow field can be inserted in the above 
expressions for the stresses, and the total force on the cylinder recalcuated. 
One finds that the second approximation makes no contribution to either XY 
or Y, so that 

X = O(e*), (42) 


Y = —2zaVpU{1+ O(e)}. (43) 
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The result (42) is not surprising in view of the argument given above, and 
(43) is in agreement with the value obtained when the circulation at large 
distances calculated in §3 is inserted in Taylor’s formula for the lift. 
Unfortunately the measurements that have been made of the forces 
experienced by a rotating circular cylinder (Goldstein 1938, pp. 545-548) 
have been made at high Reynolds numbers where the stability of the flow 
is in doubt, so that no comparison of the above theory with experiment 
seems to be possible. 

The viscous dissipation in the boundary layer absorbs the work which 
must be done against the torque to maintain the cylinder’s rotation. ‘The 
torque G can be calculated from (37), and, using the second approximation 
one finds that 


G(R) = —tmuav | 1+ 5 | (44) 

+p(p*— 1) 
For large values of R, p may be replaced by its asymptotic expression, and 
G ~ —4npaV(14+ 4v2eR??), (45) 


The author is indebted to Dr I. Proudman for his encouragement and 
advice, to Mr J. R. A. Pearson for some stimulating discussions on the 
nature of the approximations, and to Mr M. B. Glauert for some valuable 
comments on an earlier draft of this paper. During the period of preparation 
of this paper the author was in receipt of a maintenance allowance from the 
Department of Scientific and Industrial Research. 
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Experiments on the onset of wave formation on a film 
of water flowing down a vertical plane 
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SUMMARY 

The supply of water to the outside of a long vertical tube was 
adjusted until the water film was seen to be only just disturbed 
by a train of travelling waves. Under these conditions the 
Reynolds number, defined as the discharge per unit width 
divided by the kinematic viscosity, was 4+. ‘The wave train was 
slightly irregular, and average values of the length and velocity 
of the waves were about 0-45 in. and 53 in./sec, the temperature 
being 19 C. 


1. INTRODUCTION 

When a thin sheet of water passes over a smooth inclined plane, such as 
a slab of rock or a tarred road, its surface is sometimes disturbed by 
travelling waves even in the absence of wind. ‘These roll waves, as they are 
often termed, occur with their crests normal to the direction of streaming, 
and photographs of those produced in highly turbulent form in an embanked 
torrent of uniform slope and cross-section were published by Cornish 
(1934, plates 21-23). Under laboratory conditions in still air numerous 
experiments have been made, mostly on vertical surfaces, and the observa- 
tions can be recorded in terms of the Reynolds number R defined by 
R = O/v, where O is the discharge per unit width of the stream and vr is 
the kinematic viscosity of the liquid. If the slope is vertical and the liquid 
is water, the film to be examined is very thin and the experimental difficulties 
are considerable. In vertical flow the motion is not obviously turbulent 
below about R = 250, and over most of this range the surface of the liquid 
is agitated by irregular ripples, single photographs of which were published 
by Grimley (1945). Dukler & Bergelin (1952), after giving a detailed 
survey of the experimental literature, reported their measurements of ripple 
profiles at Reynolds numbers down to about 196 which showed steep-fronted 
slopes with heights that exceeded the mean film thickness. But most of the 
investigations have been concerned with finding the minimum Reynolds 
number at which waves can be observed. ‘The results show a marked 
scatter. Friedman & Miller (1941) gave R = 6-25 as the limit, and Grimley’s 
value (1945) was 6:2, but his photograph obtained at R = 6:65 showed 
waves that were apparently rather large as well as irregular. On the other 


hand Kirkbride’s result (1934) was 2. 
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In the next section an account is given of an attempt to determine not 
only the smallest Reynolds number at which waves could be seen under 
good conditions of illumination but also the velocity and length of these 
waves. ‘The liquid used was water, which was caused to flow down the 
outside of a long vertical glass tube. ‘The thickness of the film was so small 
in comparison with the radius of the tube that the motion could be taken 
to be the same as that down a plane surface. 


2. DESCRIPTION OF EXPERIMENTS 

The apparatus, which is shown in diagrammatic elevation in figure 1 
(plate 1), was built round a fixed 5 ft. length of T'ruebore glass pipe of 
outside diameter 1-045in. Distilled water passed from an elevated 
reservoir through a control valve and two alternative Rotameters into the 
upper Perspex chamber. A hole in the thick base of this chamber was 
slightly larger than the glass tube, and it could be centred by four adjusting 
screws threaded in the wall and pressing against the glass tube. In this way 
the water was delivered to the outside of the tube in as uniform a stream 
as possible. At the bottom the water was caught in another chamber and 
fell into a lower reservoir whence it was returned to the upper by a pump. 
The two chambers were joined by a Perspex pipe, held in position by rubber 
tubing. ‘The pipe served as a dust cover to the working section, and in it 
a large window was cut which was removed only when photography was 
actually in progress. Four sprays were added at the extreme top of the 
working section so that, when required, the surface of the working section 
could be thoroughly wetted. ‘hey were supplied from a branch on the 
delivery side of the pump. For observation by eye, illumination at about 
30° incidence at the middle of the working section was provided by bulbs 
and reflectors. A closely-fitting cylinder of glossy white paper was pushed 
down the inside of the glass tube; and when the waves were small, only 
their shadows on this paper were visible. ‘The wave system extended all 
down the tube, and the motion appeared to be independent of the distance 
from the top. In spite of the care taken with cleaning and with precautions 
against the entry of dust, much experience with the apparatus was required 
before it would operate consistently at very low velocities. 

The whole of the photographic work was very kindly carried out by 
Dr R. H. J. Brown, who used the repetitive flash apparatus that he has 
already described in detail (1952). Here it is sufficient to state that the 
flash duration of the lamp was about 10 microseconds, and that up to 
95 pictures per second could be taken; thus the range of his apparatus 
was amply sufficient to meet the requirements of the present problem. 
The focused beam from the lamp was directed upwards at about 20 
incidence, and it also illuminated a clock and a vertical scale placed close 
to the working section, which was photographed horizontally. 

The procedure, when the final photographs were taken, was as follows. 
After the apparatus had been running at a moderate velocity for half an hour, 
the pump was stopped and the discharge slowly reduced while the tube 
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Figure 1. Elevation of apparatus. Figure 2. Shadow of limiting train of 
Waves. 
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was repeatedly re-centred to produce waves of apparently the same amplitude 
all round the periphery. At the lowest velocity tried, waves could barely 
be seen on the front of the tube and over some but not all of the remainder 
of the periphery. ‘The discharge was then very slightly increased, causing 
waves all round; and thus it was checked that the waveless zone had not 
dried off. When the velocity was again reduced to its previous value, it 
was seen that the wave pattern was the same as before. After this process 
of raising and lowering the flow had been carried out several times, it was 
thought that the minimum condition had been attained. ‘The eye could 
just perceive a continuous stream of waves with horizontal crests, all moving 
with the same velocity but having slightly different lengths. Occasionally 
a larger and faster wave appeared that passed through the procession 
without permanently disturbing it. ‘Iwo bursts of exposures were then 
taken and the flow determined by diverting and weighing the discharge 
from the base of the apparatus. ‘The discharge per inch width of the tube 
circumference was 6:9x 10-%in.3/sec at a temperature 19°C, thus 
R = 6-9/1-58 = 4-4. 

Eighty-three pictures were obtained at intervals of 0-047 sec, and 
they were examined with the aid of a projector. ‘The specimen, reproduced 
in figure 2 (plate 1), shows one of the most regular wave trains that were 
photographed. Measurements were confined to the top 3 in. of the field 
of view, about 12 in. below the start of the working section, because there 
the shadows were most sharply defined and the distortion due to lack of 
parallelism in the light beam was least. ‘The progress of typical single 
waves could be traced without difficulty on successive pictures because 
of the slight irregularities in the wave spacing, and recognition was helped 
by small density variations in some of the shadows, which suggested 
departure from the simple sinusoidal profile. Five crests followed in this 
way were found to travel 1-57, 1-53, 1-62 in. in 6 frames and 1-83, 1-87 in. 
in 7 frames, giving velocities 5-6, 5-4, 5-7, 5-6, 5-7 in./sec. In one picture 
all the waves had lengths between 0-45 and 0-5 in., except for one that was 
only 0-31 in. from its follower. A train of four equally spaced crests, 
occupying a distance 1-30 in., moved 1-60 in. in 6 frames, thus the wave- 
length was 0:43 in. and the velocity 5-7 in./sec. Average values of the 
length and velocity were therefore about 0-45 in. and 55 in./sec. 


In addition to Dr R. H. J. Brown, whose help is described above, I have 
also to thank Professor 'T’. R. C. Fox for the loan of the hydraulic apparatus 
and Professor J. M. Kay for assistance in the early stages of the work. 
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SUMMARY 

‘This paper deals theoretically with a problem of hydrodynamic 
stability characterized by small values of the Reynolds number R. 
‘lhe primary flow whose stability is examined consists of a uniform 
laminar stream of viscous liquid running down an inclined plane 
under the action of gravity, being bounded on one side by a free 
surface influenced by surface tension. ‘The problem thus has a 
direct bearing on the properties of thin liquid films such as have 
important uses in chemical engineering. 

Numerous experiments in the past have shown that in flow 
down a wall the stream is noticeably agitated by waves except 
when R is quite small; ona vertical water film, for instance, waves 
may be observed until R is reduced to some value rather less 
than 10. ‘The present treatment is accordingly based on methods 
of approximation suited to fairly low values of R, and thereby avoids 
the severe mathematical difficulties usual in stability problems at 
high R. ‘The formulation of the problem resembles that given 
by Yih (1954); but the method of solution differs from his, and 
the respective results are in conflict. In particular, there is dis- 
agreement over the matter of the stability of a strictly vertical 
stream at very small R. In contrast with the previous conclusions, 
it is shown here that the flow is always unstable: that is, a class 
of undamped waves exists for all finite values of R. However, the 
rates of amplification of unstable waves are shown to become very 
small when R is made fairly small, and their wavelengths to become 
very large; this provides a satisfactory explanation for the apparent 
absence of waves in some experimental observations, and also for 
the wide scatter among existing estimates of the * quasi-critical’ 
value of R below which waves are undetectable. In view of the 
controversial nature of these results, emphasis is given to various 
points of agreement between the present work and the established 
theory of roll waves; the latter theory gives a clear picture of the 
physical mechanism of wave formation on gravitational flows, and 
in its light the results obtained here appear entirely reasonable. 

‘The conditions governing neutral stability are worked out to 
the third order in a parameter which is shown to be small; but a 
less accurate approximation is then justified as an adequate basis for 
an easily workable theory providing a ready check with experiment, 
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This theory is used to predict the value of R at which 
observable waves should first develop on a vertical water film, 
and also the length and velocity of the waves. These three 
predictions are compared with the experimental results found 
by Binnie (1957), and are substantially confirmed. 


1. INTRODUCTION 


The flow of liquid in a thin film down a solid wall may often be observed 
in everyday life, as when rain water runs in a sheet down a window pane, 
or when paint, or mayonnaise, drains from some solid object which has been 
dipped in the liquid. It is also a subject of practical importance in chemical 
engineering, and has been studied by many experimental workers concerned 
with that field (e.g. Kirkbride 1934; Friedman & Miller 1941; Grimley 
1945; Dukler & Bergelin 1952). ‘The character of the flow has been shown 
to depend largely on the Reynolds number R = O/v, where Q is the rate 
of volume flow per unit span of the stream and v the kinematic viscosity 

= viscosity j/density p), although the influence of surface tension clearly 
may be important in most of the cases studied. In flow down a vertical 
plane, for example, the motion is apparently turbulent when R is greater 
than about 300 (Jeffreys 1925). When R is less than this, the mean flow 
is evidently governed by a law of laminar friction, and the mean depth is 
found to be approximately that given by the simple theory due to Nusselt 
(1916) and Jeffreys (1925) which assumes a uniform flow. Nevertheless, 
waves are to be observed on the free surface throughout almost the entire 
range of laminar flow, the possible exception being when R is reduced below 
a certain small value which may be about 4 if the liquid is water. 

The apparent absence of waves on very thin films has led several 
investigators to assume that, for the flow down a vertical plane, there exists 
a critical value of R below which uniform laminar flow is entirely stable: 
that is, the flow is in a condition where small disturbances of every kind are 
suppressed. ‘The well-proven significance of the Reynolds number in its 
more familiar contexts has undoubtedly encouraged this belief. Indeed, 
numerous estimates of the supposed critical Reynolds number have been 
stated (see Binnie (1957) for a review of some of them). On the theoretical 
side the work of Kapitza (1948) led to an estimate of 5-8 for the critical 
value, and that of Yih (1954) gave about 1:5. 

In the present paper it will be argued that a critical Reynolds number 
in the usual sense does vot exist for the particular case of uniform flow down 
a vertical plane. In other words, for all finite Reynolds numbers there is 
a class of wave-like disturbances which undergo unbounded amplification 
according to a linearized theory. ‘The presence or absence of surface 
tension does not alter this general conclusion. ‘The theory is, however, 
far from being incompatible with the various experiments seeming to 
indicate a critical Reynolds number, since an alternative explanation is 
readily forthcoming for the absence of observable waves. In fact the new 
explanation demonstrates a unity among various experimental results where 
2P2 
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the respective values of the supposed critical Reynolds number are con- 
flicting. On the other hand, it is necessary to dispute previous theoretical 
vork on this subject. 

‘he problem to be considered could be adapted to a simplified treatment 
akin to the usual theory of ‘roll waves’ on turbulent streams. In the latter 
only the ‘ overall’ features of the flow are considered: that is, the frictional 
force against the stream is taken to vary in some prescribed way with the 
local depth and mean velocity and to be independent of the finer details 
of the motion, while the action of gravity (whose component perpendicular 
to the flow tends to keep the free surface flat) is expressed in a momentum 
equation based on the mean velocity. In his theory of roll-wave formation, 
Dressler (1952) derived a stability condition for a generalized law of friction, 
which extends to the example of laminar flow in question here. ‘Thus, as 
the terms representing the restoring forces can easily be modified to 
encompass surface tension as well as gravity, immediate use may be made 
of his result to estimate when waves should appear. In $5 of this paper 
the result following from Dressler’s theory is shown to correspond to a 
limiting case of the present theory; but it is appropriate to note here that 
an approach to the problem directly by way of Dressler’s work effectively 
demonstrates the instability inherent at all Reynolds numbers when the 
plane is vertical. ‘his is particularly encouraging since roll-wave theory 
keeps a much closer connection with physical reasoning than is possible 
in the present treatment. 

Since the primary (waveless) flow is laminar and completely known, the 
stability problem can be formulated mathematically in a more precise way 
than that of roll-wave theory. ‘The formulation is made in the well-known 
manner which has been extensively used for stability problems concerned 
with laminar tlows between fixed boundaries (cf. Lin 1955). It is implicitly 
assumed that the Fourier components of an arbitrary small disturbance 
are dynamically independent; thus the solution for a wave of arbitrary 
yeriod and velocity constitutes a complete solution. ‘The equations of 
motion for a wave of small amplitude are then set up, together with appro- 
priate boundary conditions; this defines a characteristic-value problem 
leading to a relationship among the essential parameters of the primary 
How and the wave. ‘There is particular interest in the conditions under 
which the wave is propagated unchanged (i.e. its velocity is a real constant), 
since these refer to the dividing line between stability and instability. 
l‘hese few sentences merely outline a very well-known type of problem, 
and we may refer to the book by Lin (1955) for an extensive account of the 
subject. The novelties of the present application are, first, the introduction 
of boundary conditions appropriate to a free surface under the action of 
gravity and surface tension, and, second, the use of a method of approximate 


solution suited to fairly small Reynolds numbers. ‘The latter avoids the 
serious difficulties which beset many of the well-known problems of 
hydrodynamic stability, in which instability is expected only for large 
Reynolds numbers (Lin 1955, p. 7). Furthermore, in a later part of this 
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paper justification is found for throwing off most of the weight of an analysis 
along the usual lines of stability theory, and hence proceeding with an 
easily workable approximation which is able to give a satisfactory account 
of all the main features of wave formation in practice. A rather similar 
formulation of this problem has already been given by Yih (1954), except 
that no account was taken there of surface tension; however, since the 
present method of solution is different and the results entirely so, there 
scems to be sufhcient reason for setting out the argument completely from 
the beginning. 

The analysis will be restricted to two-dimensional wave disturbances 
in a vertical plane (x,). ‘This simplification may be justified by appeal 
to the well-known argument first given by Squire (1933; restated by Lin 
1955, §3.1) in connection with the problem of disturbed flow between 
parallel planes. He pointed out that any three-dimensional disturbance 
is governed by the same equations as a certain two-dimensional disturbance 
in a similar flow at lower Reynolds number. ‘This means that two- 
dimensional waves have a greater tendency to instability than three- 
dimensional ones; and it follows that a two-dimensional analysis is 
completely adequate if only the stability of a flow is in question. ‘The 
suitability of a two-dimensional analysis is further demonstrated by Binnic’s 
experiments, noted immediately below, which showed that waves at fairly 
small Reynolds numbers are very approximately uniform along the horizontal 
line of their crests. 

The experimental results reported by Binnie (1957) have been found 
to provide confirmation of the theory in several important respects. ‘lhese 
results were kindly made known to me during the course of my theoretical 
work, and I was encouraged by them to explore the long-wave approxi- 
mation described in §5. ‘The experiments were made with vertical wate 
films, and their results bear out the theoretical predictions on the following 
three counts: (1) the Reynolds number at which large amplifications of 
unstable disturbances first occur, so that observable waves develop; (ii) the 
velocity of the wave of maximum instability; (iii) the corresponding wave- 
length. ‘lhe theoretical argument also accounts for the irregularity of the 


observed waves. 


2. PROPERTIES OF THE PRIMARY FLOW 

The laminar flow whose stability is to be examined is a uniform two- 
dimensional stream bounded on one side by a fixed wall and on the othei 
by a free surface; the fluid has constant density and viscosity. ‘The velocity 
U is everywhere parallel to x (see figure 1), and the graph of U ws y is 
parabola which has its vertex in the free surface since, as the fluid is viscous, 
the rate of shearing dU /dy must vanish there. In this paper the coordinates 
(x,y) are made non-dimensional by taking the stream depth / as the unit 
of length, so that y = 0 and y = 1 define the free surface and the fixed 
boundary respectively. Also, velocities are made non-dimensional by takin; 
the undisturbed velocity at y = 0, denoted by up, as the unit of velocity. 
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Thus we have U=1-y*. (2.1) 
By integration this immediately shows that the dimensionless mean velocity 
is u,, = % (i.e. in dimensional units u,,= O/h = Zug). ‘he motion is 


steady, so that the shear force —(ju,)/h)dU/dy on the wall is balanced by 
the total gravity force on the stream in the direction of flow; thus 


2uuy/h = phg sin 6, 
where @ is the slope (figure 1). ‘This gives directly 
uy = Sh®gsinG/v, ( 


R = Ov = Ruyh/v = sh®gsin 6/r*. (2.3) 








Figure 1. Diagram of the undisturbed flow, showing the velocity profile. 


The existence of an explicit relation (2.3) between flow rate and depth 
suggests that any gradual or ‘long-wave’ disturbance from the uniform 
state would be propagated downstream with a (dimensional) velocity 
C = dO/dh (cf. the account of ‘kinematic waves’ given by Lighthill & 
Whitham (1955)). In fact (2.3) leads to 


C = h®gsin@/v = 2up. (2.4) 
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‘This simple result provides an interesting check on the subsequent analysis, 
which shows by a totally different method that very long periodic waves 
are indeed propagated with a phase velocity equal to 2u. 


3. "THE STABILITY PROBLEM 
Suppose that the free surface is given a perturbation represented by 
the equation 
Y = (x, t) = dee, (3.1) 
where x is the dimensionless wave number (= 2zh/wavelength), and 
t is time made non-dimensional by taking h/u, as the unit. Either the real 
or the imaginary part of (3.1) may, of course, be taken to describe a wave 
on the physical surface. In the usual manner of stability analyses, the 
wave amplitude 6 is taken to be a small quantity whose square is negligible. 
‘The corresponding perturbation of the stream function may be written as 
— 8 f(y)ee—, 
which implies that the velocity components satisfying the continuity 
condition for incompressible flow are 


u = 5f’(y)e! u(c—cl) 1 r r ] 
4 


wr 
bo 
— 


v= — 148 f (ye ct), | 
The elimination of the pressure from the linearized equations of motion 
leads to an equation for f(y): 
(U-—c)(f" —#f)-U"f = — (f* —2ef"+2'f), ci 
Mf" —atf)— Uf = = (f"—2atf” +24f) 
which is commonly known as the Orr-Sommerteld equation (Lin 1955, 
$1.3). After substituting for U from (2.1), it is convenient for our purpose 
to rearrange (3.3) in the form 


f® = (n—cn + 2a? —ny*)f" + (2n—o2n + co2n— 04+ 2ny*)f, (3.4) 
and then f* =(pt+ gy" sf" +(r+sy")f, (3.5) 
where we have first put m = 1au,h/v = 3ia2R/2, and then 
p =n-—cn+22?; g=—-N; 
r (3.6) 
r = 2n—2?n+co2n—a!; f= ata: | 


This greatly simplifies the writing-out of the subsequent heavy algebra. 
The problem entails five boundary conditions as follows. First, the 
kinematical surface condition is 


=. Os 


where U, and vy are the values of U and v at y=0; hence, because of 
(2.1), (3.1) and (3.2), 
f(0) =e-1. (3.7) 
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This equation simply shows the connection between the assumed perturba- 
tion of the free surface and the corresponding perturbation of the stream 
function; it bears out the intuitively reasonable expectation that, provided 
c is not excessively large, the two perturbations are of the same order of 
magnitude. Since both velocity components must vanish at the fixed 
boundary, we also have 


f(Q1) = 9, (3.8) 
f'(1) = 0. (3.9) 


The remaining boundary conditions relate to the continuity of stresses 
across the free surface. It is assumed that the only stress acting on the 
outside of this surface is a uniform pressure py; therefore, since the fluid 
is viscous, the rate of shearing in the stream cu/ey+cv/ex must vanish 
at y = 7, which requires that 


f"(0) = 2—a*(c— 1). (3.10) 


The equation expressing the continuity of normal stress across the 
free surface may be written, after division by the density p, 
>) a. 
LVu,, [C7 = 
Oe oe = & +01 ~g cos 4 n, (3.11) 
p h OV] y=0 p 
where p, is the (dimensional) pressure just inside the stream, o 1s the principal 
curvature of the surface, and I’, which may be called the ‘ kinematic surface 
tension’, is equal to the usual surface tension coefficient divided by density. 


To the first order in 6, the curvature is 


1 0% 470 
/ : ix(r—e 2 45 
o=7>5- —e 3.12) 
h ox? h ( 
Also, the pressure p satisfies the Navier-Stokes equation 
eu ou ou l cp vy feu cu) a ete ie 
Se AAG aes ARE on — = t+ ea t+ ap thesind/us. (3.13) 
ot Ox oy purcx hits \ x ove f 


Hence, with the use of (3.11) and (3.12), separation of the variable part 

of (3.13) (i.e. the part with 6 as a common factor) at vy = 7 leads to 
(n —cn +302) f'(0) = f'"(0) 4 — (231 + xh®e cos 8). (3.14) 

0 

The first four of the boundary conditions, (3.7) to (3.10), will be used 
to determine the four constants arising in the solution of the fourth-order 
equation (3.4). ‘he boundary condition (3.14) will then be used to find « 
as a function of « and the parameters of the mean flow. ‘This procedure 
resembles the usual one in stability problems (cf. Lin 1955). Instability 
is revealed for a particular wave mode when the respective ¢ turns out to 
have a positive imaginary part, for then the mode has an exponentially- 
increasing time factor. As usual it is convenient to express the results of 
the calculations by setting down the conditions which make ¢ wholly real ; 
cis then simply the phase velocity of a neutral disturbance. ‘These conditions 


represent a division between stability and instability. 
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‘he present problem is somewhat unusual in that two stress conditions, 
(3.10) and (3.14), replace the two additional kinematic conditions, like (3.8) 
and (3.9), which arise in problems concerning flow between fixed boundaries. 
It may, incidentally, be asked why five boundary conditions have been 
written down whereas most treatments of stability problems (cf. Lin 1955) 
make do with four? However, the limited meaning of the first boundary 
condition (3.7) has already been noted above; it contributes nothing to the 
essential characteristic-value problem, since an assumed perturbation of 
either the free surface or the stream function may equally well be taken as 
a starting point, and one simply follows from the other according to (3.7). 

The only important difference between this and the better-known 
problems of hydrodynamic stability lies in the method of solution of the 
Orr—Sommerteld equation set out in the next section. In view of existing 
experimental evidence, attention may confidently be restricted to the case 
of fairly small Reynolds numbers and wave numbers, and indeed the end 
results of the calculation confirm the adequacy of this approach. 


4. POWER-SERIES APPROXIMATION TO THE STREAM FUNCTION 

It was Kelvin (1887) who first noted that in linearized problems like the 
present one where the function f(y) satisfies the fourth-order differential 
equation (3.4), the solution may be expressed as an ascending power series 
in y. He demonstrated that this series is convergent for all physically 
realizable conditions, but observed that the rate of convergence is very slow 
when R is not small. ‘The approach suggested by Kelvin is of little use 
at fairly large Reynolds numbers, with which the greater part of the existing 
work on hydrodynamic stability is concerned. For the present problem, 
however, the method is definitely useful, since the values of R and of « are 
expected to be fairly small in the range of physical interest. 

If one formally puts 


fo) = > Avy, (4.1) 


N=vU 
then this series is seen to constitute a solution of (3.5) when the coefficients 
A,, are made to satisfy the recurrence relation 
N(N—1)(N—2)(N—-3)Ay = (N-—2)(N—3)pAy_2+ 
-{r+(N—4)(N—5)q}Ay_,t+s*Ay_¢ (4.2) 


for N > 3. ‘This relation follows directly from the differential equation. 
The first four coefficients A, to A, are simply the four constants arising 
essentially in the solution of a fourth-order equation, and (4.2) gives every 
other A, in terms of these four and the parameters p, q, 7, s. If these 
parameters are now regarded as small quantities of the same order of 
magnitude, the four coefficients A, to 4; may be said to be O(p), and each 
successive set of four coefficients is of decreasing order of magnitude with 
reference to powers of p. From (3.4), (3.5) and (3.6) it is seen that 
provided « is small, then p, q, r, s are of the order of m: that is, the order 
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of «R. ‘Thus, a limited number of terms of the series would amount to 
an approximate solution in terms of ascending powers of «R (and, 
incidentally, of x”); this may be contrasted to the commonly-used asymp- 
totic solutions in descending powers of «R, suitable when R is large (see, 
for instance, Lin 1955, chapter 3). 

We choose to approximate f(y) in this way as far as O(p*) = O(n’), 
which requires sixteen terms of the expansion (4.1). In this procedure 
there is an implicit assumption that ¢ is not large, which is confirmed by 
the final results. In order simply to keep symmetry in the approximation, 
terms which are O(z?), for instance, are given equal status with terms which 
are O(n), although it may be anticipated that a < mn (i.e. « < R) for 
the physically interesting cases discussed later. 

Although the various manipulations made with the approximate solution 
for f(y) are very laborious, they are quite straightforward. Accordingly, 
only an outline of the method of calculation need be given, together with 
the results of the more important stages. 

The third-order approximation to the series solution (4.1), as calculated 
by means of (4.2), is most conveniently arranged as follows: 


f(y) = Ay[1 + (1/4 !)rvt + (1/6!)( pr + s?)y§ + (1/8 !)(12@r + 1? + p*r + 2ps?)v> + 
+ (1/10!)(2pr? + 42pqr + 60qs? + 32rs*)y™ 4 
+ (1/12 !)(672q?r + 68qr? + ny m4 + 
+ Afy+(1/5!)ry? + (1/7 !)(pr + 6s?) 
(1/9!)(20gr +r? + p?r 4 oy 
+ (1/11 !)(62pgr + 2pr? + 252gs? + 48rs*)y"! 
+ (1/13 !)(1440g"r + 92qr? + r°)y"9] 4 
A,| y? + (2/4!)py? + (2/6!)(2q¢+7r + p?)y® 4 
+ (2/8!)(14pq + 2pr + 12s? + p*)y* 
+ (2/10!)(60g? + 32gr + 1? + 44-p?q + 3p?r + 42ps*)y!? + 
+ (2/12 !)(560pq? + 10pgr + 672qs" + 12rs*)y!?] + 
+ A,[ v3 + (6/5 !)py? + (6/7 !)\(6g +r + py" 4 
+ (6/9!)(26pq + 2pr + 20s? + p)y? + 
+ (6/11 !)(2529? + 48qr + r? + 68p?q + 3pr 4 ballast 
+ (6/13 !)(2124pq? + 218pqr + 3pr? + 1872gs? + 92rs?)y + 
+ (6/15 !)(2772093 + 5532g?r + 158qr? + r)y??]. (4.3) 
To use this result, the boundary conditions (3.7), (3.10) and (3.14) are first 
expressed in terms of the coefficients 4, ; these equations give respectively 


Ay =¢-1, (4.4) 
2A, = 2—2a%(c—1), (4.5) 
(n—cn+3*)A, = 6A,4 ae (a3 + ah*g cos 8). (4.6) 


vuy 





si Sl AA Oca 
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Thus, 4, and A, are related very simply to the main parameters of the 
problem. ‘To find A, and A, for use in (4.6), the boundary conditions (3.8) 
and (3.9) are employed. By setting equal to zero f(y) and f’(y) as given 
by (4.3), we have a pair of simultaneous equations for 4, and Az, in terms 
of A, and A,. ‘The solutions of these equations are then expressed in 
ascending powers of Pp, q, r, s by a straightforward, though rather lengthy, 
inversion process. By this stage of the calculation, some of the numerical 
coefficients would become very cumbersome if left as fractions, and so all 
but the simple coefficients of the leading terms are expressed as decimals. 
Note that to satisfy (4.6) as far as O(n*), A, is required to O(n*), but A, is 
required only to O(n?). Accordingly, although A, was calculated to O(n*) 
to provide certain checks on the result for A3, the third-order terms in A, 
need not be reproduced here. 

After some condensation made possible by the fact that r = —2g in 
the third-order terms, the results may be expressed as follows: 


A, = Af a. ( + Bf r) + (—0-0013 095p? — 


40°* TH0%* TOs 
— 0-0006 052pq — 0-0003 373 pr — 0-0000 7734? — 


bo] Ww 


— 00000 167gr + 0-0000 3577? + 0001091382) | + 
_— 


] . . ! ; 0-Q007 73S8p2 
tAel — 5 +) egPt ga02- SO") aia i 
— 00-0001 389pq +(0-Q001 190pr — 0-000 054g" + 


+ 0-0000 123gr — 0-0000 0497? — 0-0001 ost) | (4.7) 


| 1 3 11 
Pe Be Se. ae Be anne 8 (OT es 
Ag = Ad 5 t ( x P 730 7 380 r) (0-0032 143p? 4 
+ 0-0012 83 1pq + 0-0010 450pr + 0-0001 3499? + 
+0-0000 320qr — 0-0000 8917? — 0:0022 487s") + 


+ (—0-0001 759p — 0-0000 315p2r + 0-0000 130pr? + 


+ 0-0000 0307? + 0-0001 040ps? + 0-0000 507%) | 4 


| 3 1 l 
+ Ad — = + — PR wee ot aa 7) + (08-0023 413 9"- 
uf 2 ( 40? ~ 3154 515") ( —" 
+ 0-0003 538pg — 0-0003 LO8pr + 0-0000 1399? - 
— ():0000 276gr + 0-0000 1197? + 0-0003 638s?) + 


+ (—0-0001 0683 — 0-0000 758p2r — 0-0000 192pr? — 


— 0-0000 020r? — 0-0000 445ps? — 0-0000 1221s | (4.8) 








564 T. Brooke Benjamin 


When p, g, 7, s are expressed in terms of n and z? by means of (3.6), thes 


results lead to 


am, | ee 
‘1, f1) = m+ gag "+ yt) 


+ (—0-0013 095c2n? + 0-0026 S85cn? — 0-0012 798n* 4 


+ 0-0147 619cx2n 0.0138 09522 — 00147 6192!) | : 


l l 11 : * 
+ Aa 3 +( a 340 "+ 39%) 


+ ( —0-0007 738c2n? + 0-0011 706cn? — 0-0004 463n7 5 


+ 0-0019 048cx2n — 0-0013 49222n - o-0019 0482") |, (4.9) 


33 Ee dae re 
As Ay] 5 + (5H 8 — 56” wo“) (0-0032 143c2n? — 0-0072 354¢n? + 


+ 00037 357n® — 0-0521 429ca*n + 0-05.15 0792Fn + 
- ()-0521 4292!) + (0-0001 759c%n? — 0-0004 963c2n? + 
+ (0004 519¢n? — 0-0001 34523 — 0-0032 712c?a2n? 4 


t ):0054 363¢72n? — 0:0022 48322? + 0-0083 413cx4!n — 


0-O068 2542!n — 0-0065 926z! ] + 


l 3 11 3 
1 "ey el per ee (ANI? 3c2n2 
A,| 5 (F< 168" x) (00023 413c"n 


0-0037 O70cn? + 0-0014 8$26n? — 0-0061 905cx2n 4 
+ 00046 03222n + 0-0061 9054) + (0-0001 068c3n? 
)-QO002 447c2n? + 0-0001 88 1en? — 0-0000 48573 


0-0003 300c2z2n2 + 0-0004 767¢c22n? — 0-0001 729%2n? + 
+ )-Q003 492¢x4n — 0-0002 103 24n — 0-O0002 3282) | (4.10) 


When these expressions are substituted into (4.6) and as much condensation 
as possible is made, the result is 





7 : 6 OS 24 a 
— (x31'+ ah? cos @) + 3c—6- nf ~¢(?+ ==C =) + 0?(5-4c — 48) + 
VUy ) Es Bb 
+ n*( — 0-0057 143c? + 0-0215 873c? — 0-0254 628¢ + 0-0102 907) 4 
+ y?n( — 0-1628 57 1c? + 0-3029 095c — 0-1397 619) 4 
+ #4( — 0-1371 429c + 0-0742 857) + 
+ m3( — 0-0002 540c? + 0-0011 412c8 — 0-0018 O07 ce? + 
+ (-0012 414¢ — 0-0003 178) + 
+ #2n*(0-0073 968c3 — 0-0240 506c? + 0-0257 OOS8e — 0-0090 556) 4 
+ a'n( — 0-0506 667c? + 0-0966 906c — 0-0468 574) 4 
~ 2°(0-0547 302c — 0-0504127) = 0. (4.11) 
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This result can be regarded as an equation for c in terms of m and z, 
hence of R and z, and of the ‘restoring force’ due to gravity and surface 
tension as expressed by the first group of terms in the equation. In 
particular, the equation can be used to find the relation between R and « 
which makes ¢ wholly real, i.e. which gives the case of neutral stability. 
To do this, we first separate the real part of (4.11), noting that ” is purely 
imaginary, and hence quite easily obtain an explicit expression for ¢ by a 


process of successive approximation in terms of x? and n? = — 9a?R?/4. 
We finally get, for ¢ real, 
c = 2(1—2? + 112!/6 + 00-0077 58122R? — 3-3555 5568), (4.12) 


where the values of R and zx must be consistent with the imaginary part 
of (4.11) being satistied. ‘The error in (4.12) is known to be considerably 
less, due to numerical factors, than O(a'R*). Note that there is no term 
in («R)! in (4.12); this result occurs by the cancellation of terms with fairly 
complicated numerical factors. (This may seem perhaps to suggest some 
special significance; but I can see no other way of arriving at such a result 
beyond a direct calculation like the present one.) 

The imaginary part of (4.11) is now considered. On eliminating c 
by means of (4.12), and using the following relations found easily from (2.2) 
and (2.3): 

uh = (3°3/4)(gsin A)! yt 3R3, 
uh = iosin OR, 


ts 
we finally obtain the approximation 
i + = - ee. 
7° 3R 5 
— 0-0000 639%? R? — 14-6352952%! = 0, (4.13) 
where € is defined by 
¢ = P(gsiné)-**y-**, (4.14) 
Equation (4.13) expresses the required approximate relation between R 
and « for neutral stability. ‘he equation is easily solved, since for a specified 
value of R it becomes a quadratic in z?. 

In passing it may be noted that the formula for real values of c, i.e. (4.12) 
subject to (4.11), appears to differ from the result obtained by Yih (1954), 
which was presented only in graphical form. ‘The discrepancy is not 
understood. For x 0 (i.e. for very long waves), equation (4.12) gives 
simply c = 2, which agrees with the result (2.4) obtained by a very simple 
argument. ‘This also agrees with the result found by Kapitza (1948), 
who used a different form of argument applicable only to very long waves. 

Let us now consider, for example, the case of a vertical wall (@ = 90 ). 
The term involving cot 4 in (4.13) now vanishes, and we have ¢ = lg"! 3y-43, 
(The parameter ¢ is about 3300 for water at 19 C.) Figure 2 shows the 
relation between R and x according to (4.13) for various values of € and a 
range of R up to 20. ‘The curves shown in the figure are curves of neutral 
stability; and it can be inferred that, for a particular ¢, the region lying 
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above the respective curve represents stability, while that below it represents 
instability. ‘Che following features of figure 2 may be specially noted. 


(i) Every curve for ¢ > 0 passes through the origin, and has infinite 
slope there. (The case ¢ = 0 is discussed below.) ‘Thus, for all finite 
values of R there is a finite range of unstable x, and therefore the flow is 
never completely stable. On the other hand, equation (4.13) shows that 
stability occurs for sufficiently small R when 0 < 4 < 90. 


Pm $:0 
me 300 
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10,000 
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Figure 2. Curves of neutral stability for laminar flow down a vertical wall with various 
values of the parameter = T g~'/3 p-4/3, 

(ii) ‘he stabilizing tendency of surface tension is demonstrated by the 
fact that the curves bend nearer the R-axis for increasing ¢. ‘Thus, as might 
be expected intuitively, the range of unstable x for a given R is reduced 
as I’, and hence ¢, is increased. Note that surface tension cannot induce 
complete stability, although it may stabilize waves whose length is less than 
a certain limit. 

(11) ‘he curve according to (4.13) for ¢ = 0 is included in the figure, 
being shown as a broken line, although the values of x in this case are rather 
too large for confidence in the accuracy of the result. (We recall that the 
theory has been developed on the assumption that x is fairly small.) 
However, the trend of the other curves in the figure shows definitely that 
the unstable region of the (R,«)-chart is largest when €=0. ‘This fact 
is useful, since it seems to rule out the possibility that the neglect of surface 
tension in Yih’s treatment of the problem may be responsible for the 
discrepancy between his and the present general conclusions regarding 


stability. 





sd, ‘eae 











Wave formation in laminar flow down an inclined plane 567 


In connection with the general role of surface tension in the mathematical 
problem, the following matter is worth noting. Equation (4.13) shows 
that « > O(R°®/¢!) tor R + 0 with 6= 90. ‘Thus, for small R the curve 
4 vs R°® is a parabola whose vertex touches the z-axis at the origin. ‘The 
curvature at the vertex is directly proportional to ¢, and so in the limit 
¢ 0 the curvature is zero and the curve coincides with the z-axis. It 
may hence be inferred that the curve « vs R also coincides with the x-axis 
for > 0, although for this curve the curvature at the origin is infinite for 
for all finite values of ¢. Clearly, the same result should occur if surface 
tension were ignored (i.e. € = 0) from the outset of the analysis, and so 
one might question why no such result is forthcoming if we put directly 
¢ = Oand @ = 90 in (4.13). However, in proceeding from the dimensional 
form of the dynamical equation (3.13) to as far as (4.13), a division of all 
the terms by u5h, which is proportional to R®*, has been entailed. Hence, 
to preserve the appropriate ‘ physical’ status of surface tension with respect 
to the limiting process R — 0), equation (4.13) should be multiplied by R®*; 
when this is done there is no ambiguity about the fact that R= 0 is a 
solution regardless of « when ¢= 0. ‘This result that the flow is unstable, 
or at least neutral, for all finite R and « seems correct intuitively, since for 
a vertical wall and zero surface tension there is no restoring force to act 
on the disturbed flow. 

(iv) Owing to the particular method of approximation used for the theory, 
the most accurate part of the figure is the vicinity of the origin. ‘The upper 
half, at least, of the figure is somewhat uncertain; but the figure as a whole 
is quite adequate as a qualitative picture of the neutral-stability conditions. 
Moreover, the line of argument to be developed in § 5 shows that it is sufficient 
to consider a range of R and x very near the origin to account satisfactorily 
for the relevant experimental facts. 


5. SIMPLIFIED THEORY FOR LONG WAVES 

In this section some formulae are derived which are suitable for a 
straightforward check with experiment. In contrast with the results of 
the last section, they are compact and easily manipulated; and although 
these advantages are gained by sacrificing the higher-order terms calculated 
in $4, the loss of precision entailed is found not to be serious. Let us 
therefore proceed on the assumption, to be justified later, that x is quite 
small. Our knowledge of the second- and third-order terms in (4.11) 
shows conclusively that a valid first approximation to this equation is 
simply 

ans 

— \2°1'+ ah? cos 6} +3c—6+4 

VUyg 


6 68 24) 
sek ge A ae a. 0. (5.1) 
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Suppose now that c = c,+ic,, where c, and c; are wholly real, and that 
c, > ¢; The latter assumption will also be justified later. On separation 
of real and imaginary parts, equation (5.1) gives, as a first approximation, 


Cy = 2, (5.2) 
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and hence 


|. td thg cos 6) — 
C; 5 R < = &@ ae. tJ.) 
- hu uz 
‘Therefore, since unstable waves have c, 0, the condition for instability 
1S 
8 271° hgcos# e 
rd : 9 T 9 (5.4) 
5° hue u2 


Some special cases following from (5.4) deserve to be noted. If one puts 
x » Qin (5.4) and substitutes (2.2) and (2.3) in the second term on the right- 
hand side, this instability condition becomes 


R> 2 cotd. (5.5) 
6 

Thus, when (5.5) is just satisfied, very long waves (z > 0) become unstable. 
Moreover, if surface tension is zero, all fairly long waves are unstable 
when (5.5) is satisfied, which rather suggests that in practice a bore (in 
other words, a discontinuous ‘roll wave’) would be formed. It is notable, 
though not at all unexpected, that except for the numerical factor the 
instability condition (5.5) is identical with the condition found by Jeffreys 
(1925) for the formation of roll waves on a turbulent stream, assumed to be 
subject to the Chézy law of friction. Indeed, the present result can be 
deduced directly, with only a very small discrepancy in the numerical 
factor, from the roll-wave theory developed by Dressler (1952) for a 
generalized law of friction. ‘The details of this comparison need not be 
given here; it is sufficient to note that the result follows in a straightforward 
way when the primary-flow relationships for laminar flow (as in §2) are 
introduced into Dressler’s generalized instability condition. 

‘The action of gravity may be expected to encourage instability if 6 > 90 , 
that is, if the fixed wall slopes backwards and the stream runs down its 
under side. ‘The ability of liquid streams to adhere beneath solid surfaces 
is popularly known as the ‘tea-pot effect’, and will be familiar to any one 
who has poured liquid from a vessel with an ill-designed spout. The 
destabilizing influence of gravity in this case is properly expressed by the 
gravitational term in (5.4), which becomes negative when 4 90. As 
an extreme case, let 6 = 180° and hence uw, = 0 in accordance with (2.2). 
‘The liquid now forms a stationary film on the under side of a horizontal 
plane, and the instability condition (5.4) becomes 

gh® > al’, 
or 


A > 2n(P/g)!2, (5.6) 


5 


if the wavelength A = 2zh/x is introduced. This is seen to check with 
the result obtained by Bellman & Pennington (1954, equation (3.7)), who 
calculated the effect of surface tension on the ‘Taylor instability of the 
boundary between two fluids accelerated perpendicularly to its plane 





‘ 
{ 
| 











Wave formation in laminar flow down an inclined plane 569 


(Taylor 1950). (It should be noted that the effect of gravity is here equivalent 
to an upward acceleration.) 

‘The final special case which we shall discuss arises when # = 90°, that 
is, when the fixed boundary is vertical. ‘This is the case for which a 
comparison between theory and experiment will be made in §6, and also 
the one in which there is marked disagreement between the present work 
and that of Yih (1954). The rest of this section will be devoted to it. 

For 48 = 90°, (5.4) becomes 

8 


Ss Ta 
5 hw 


Hence, the flow is always unstable, since sufficiently small values of « can 
always be found to satisfy this instability condition; indeed, this fact is 
already clear from the discussion at the end of § 4 and also from (5.5), which 
becomes simply R > 0 when @ = 90. As pointed out in an early part of 
this paper, this result is remarkable in view of existing experimental work 
on vertical liquid films, much of which seems to indicate an absence of waves 
at sufficiently low Reynolds numbers. However, the following straight- 
forward argument appears to bring the present theory into a completely 
satisfactory agreement with experiment, in particular the work mentioned 
in $6. 

It is reasonable to suppose that the wave most prominent in an 
experimental observation will be the one whose rate of growth according 
to linearized theory is largest. ‘The use of such an assumption has many 
precedents; for instance, Rayleigh (1894, $87) discussed its theoretical 
basis, and used it in several of his investigations. ‘There is various experi- 
mental evidence to suggest that often the growth of a ‘wave of maximum 
instability’ eventually, through non-linear effects, tends to suppress other 
waves which are unstable according to linearized theory but have smaller 
rates of growth; thus, the well-defined pattern of the particular wave mode 
of greatest rate of amplification may be observed (see Benjamin & Ursell 
1954 for a discussion of this point in a rather similar context). In the 
present problem unstable waves grow in amplitude according to a time 
factor exp(zc;t). ‘lo tind « for the most unstable wave, we have therefore 
only to maximize xc;, which, by virtue of (5.3), amounts to maximizing the 


is 


function 


Hence, the optimum value of x (say %,,) is given by 


i thu? — 
“e, = 3 (5.5) 


m 51° 





When the formulae of §2 for the primary flow are used (with @ = 90°) to 
eliminate u, and / from this expression, it leads directly to 


Lm, = LAZLr2BgV6P 12) RIE, (5.9) 


m 
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where the numerical factor is an approximation to (3/5)!73¥8, This together 
with (5.3) shows that the maximum value of ¢; is 

(c,),, = 0-224{ 2% 3gl6] V2) RU, (5.10) 
lhe quantity within the braces in (5.9) and (5.10) depends only on the 
particular liquid considered. For example, taking the values appropriate 


to water at 19 C: v=0-0103 cm*,sec and VT = 72:9 cm?/sec?, and also 
g = 981 cm/sec”, we get from (5.9) and (5.10) 

ty, = 00195 R°®, (5.11) 

(c;), = 0-00391 RY *, (5.12) 


Provided R is not much greater than unity, these results justify the initial 
assumptions that x and ¢; are small. 

Now, to illustrate the likelihood of waves being observed experimentally, 
we shall calculate the amplification experienced by the wave of maximum 
instability on a vertical water film as it travels a distance of 10cm. Since 
the wave velocity is approximately 2u,, the amplification factor is approxi- 
mately ./ = exp{10z,,(c;)),,/24}, where h# is in cm. Hence, using (5.11) 
and (5.12), and eliminating 4 by means of the primary-flow relations of § 2, 
we finally get 

J = exp{0-0543R°3}. (5.13) 
Values of 7 given by (5.13) are plotted as a function of R in figure 3. The 
tigure shows the amplification to be large only for values of R greater than 
about +; above this the rate of increase of .c7 is enormous. ‘Thus, one might 
expect observable waves to develop reasonably near the beginning of the 
stream as soon as the Reynolds number is brought up to a value near 4. 
Of course, the linearized theory ceases to apply when large amplifications 
occur; but the present result does seem to offer a convincing explanation 
of why, in experiments where the Reynolds number is gradually increased 
from small values, waves should suddenly appear on a previously untroubled 
stream. However, it has been shown here that a critical Reynolds number 
in the usual sense does not exist for a vertical stream; for, even when the 
Reynolds number is extremely small, there are unstable waves which, 
presumably, would grow to appreciable size if the stream were long enough. 
lhe present result suggests that experimental measurements of the * quasi- 
critical’ Reynolds number (R,,. say, which seems from the theory to be 
about + for water) may depend to a large extent on the method and precision 
of the observations (even perhaps on the eyesight of the observer); this 
may explain the wide scatter among existing estimates. Note that the value 
of R,. depends in a fairly complicated way on the kinematic surface tensioi 
and viscosity; according to the present theory, it is definitely not the sam« 
for different liquids. 

In conclusion a formula for the dimensional wavelength A, corre- 
sponding to the most unstable wave number z,, is worth noting, being mort 
convenient than (3.9) for comparison with experiment. Putting A,, = 27h z 


ind using (2.3) and (5.9), we find that 
A, » Sir Rs. (5.14 
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For water at 19° C, this gives 


A,, = 2:21R-™* cm = 0-871 R-? in. (5.15) 


m 
This is suggested as the most likely wavelength of observed disturbances, 
at least when instability first becomes noticeable, and a wave structure can 
be distinguished. 


60} 
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Figure 3. The amplification experienced by the most unstable wave on a vertical 
water film as the wave travels 10 cm. 


Nevertheless, one can scarcely expect waves to appear with a strictly 
uniform and distinct periodicity ; because under all conditions infinitesimal 
waves with a wide range of wavelengths are unstable, and the wave with 
length A,, comes into prominence only through a rather uncritical selection 
process depending on differences in the rates of amplification at different 
wavelengths. ‘he ultimate state of the amplified waves is, of course, 
determined largely by non-linear effects which remain unknown. 


6. COMPARISON WITH THE EXPERIMENTS BY BINNIE 


‘The reader may refer to the paper by Binnie (1957) for an account oi 


these experiments, which were carried out with water at 19 C. ‘The waves 


2Q2 
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described by him were somewhat irregular in length even when the flow 
was delicately adjusted until the waves were just perceptible (see the 
photograph in his paper). ‘This feature is understandable from the con- 
siderations made in the final paragraph of § 5, and suggests the unlikelihood 
of a very precise agreement between experiment and any linearized theory. 

‘The three main points of comparison are as follows: 

(i) ‘The Reynolds number at which observable waves first occur was 
estimated experimentally as 4-4. Although the theory does not provide a 
correspondingly precise estimate, this value appears quite reasonable in the 
light of equation (5.13) and figure 3. 

(ii) ‘The records of wave velocity and length were taken when the 
measured value of O was 6-9 x 10°? in.?/sec per inch span. On the substi- 
tution of this value and the value of v quoted in §5 for water at 19 C, 
equations (2.2) and (2.3) (with 6= 90) lead to h = 44x 10-%in. and 
hence 2u, = 4:7 in./sec. According to (5.2) the latter figure is the first 
approximation to the wave velocity. The mean experimental value for the 
wave velocity was 5-5 in./sec, i.e. 17°, in excess of the theoretical estimate. 
The order of magnitude of the extra terms involved in an approximation 
toc, better than (5.2) is found to be too small to account for this discrepancy, 
which seems more likely to be due to the finite size of the measured waves. 
(‘The velocities of other kinds of surface wave are known to increase with the 
wave amplitude; e.g. irrotational waves in an inviscid fluid (Lamb 1932, 
$250).) However, the comparison seems quite good in view of the obvious 
limitations of linearized theory in this type of investigation; in any case, 
it appears definitely to support the present theoretical results in favour of 
those found by Yih (1954) which indicate that the wave velocity should be 
many times greater than the experimental values. 

(iii) With R = 4-4 as measured, equation (5.15) gives A,, = 0-42 in. for 
the most unstable wavelength. ‘This is in surprisingly close agreement 
with the value 0-45 in, obtained as the mean of the experimental observations. 

‘The above three items seem in large measure to confirm the usefulness 
of the simplified method developed in $5. However, further experiments 
would be necessary for a complete check; for instance, measurements 
of RX, for different liquids would provide a particularly severe test. Suitable 
experiments on films with a finite slope would be more difficult than those 
on vertical films, since the convenience of a cylindrical wall is lost and 
difficulties with edge effects are bound to arise, but such experiments would 
be very desirable. 


| wish to acknowledge the encouragement given throughout this work by 
Mr A. M. Binnie, and his invaluable help in checking the lengthy algebraic 
calculations of §4. I should like also to express my thanks to Professor 
T. R. C. Fox, who first stimulated my interest in this problem, and to 
Dr G. K. Batchelor, who noticed that a powerful criticism based on energy 
considerations could be brought against my results for a vertical stream ; 
this proved indecisive in the end, but led to a helpful clarification of some 


aspects of the problem. 





i 
} 
: 











ws 


Wave formation in laminar flow down an inclined plane 


REFERENCES 

BELLMAN, R. & PENNINGTON, R. H. 1954 Quart. Appl. Math. 12, 151. 

BENJAMIN, T’. B. & Urseti, F. 1954 Proc. Roy. Soc. A, 225, 505. 

BInNIE, A. M. 1957 ¥. Fluid Mech. 2, 551. 

Dresster, R. F. 1952 Gravity Waves, Nat. Bur. Stan., Wash., Circular no. 521, 
p. 257. 

DuKter, A. E. & BeRGELIN, O. P. 1952 Chem. Engng Progress 48, 557 

FRIEDMAN, S. J. & Mitver, C. O. 1941 Industr. Engng Chem. 33, 885 

GRIMLEY, 8. S. 1945 Trans. Inst. Chem. Engrs 23, 228. 

JerFreys, H. 1925 Phil. Mag. (6), 49, 793. 

Kapitza, P. L. 1948 }. Exp. Theor. Phys., U.S.S.R., 18, 3. 

KELVIN, Lorp 1887 Mathematical and Physical Papers, Vol. 4, 321. Cambridge 
University Press. 

KIRKBRIDE, C. G. 1934 Industr. Engng Chem. 26, 425. Trans. Amer. Inst. Chem. 
Engrs 30, 170. 

Lams, H. 1932 Hydrodynamics, 6th Ed. Cambridge University Press. 

LIGHTHILL, M. J. & WuiTHaM, G. B. 1955 Proc. Roy. Soc. A, 229, 281. 

Lin, C. C. 1955 The Theory of Hydrodynamic Stability. Cambridge University 
Press. 

NussE_T, W. 1916 Z. ver. Dtsch. Ing. 60, 541. 

RAYLEIGH, Lorp 1894 The Theory of Sound, 2nd Ed., Vol. 1. London ; Macmillan 

Seuire, H. B. 1933 Proc. Roy. Soc. A, 142, 621. 

"TayLor, G. I. 1950 Proc. Roy. Soc. A, 201, 192. 

Yiu, C.-S. 1954 Proc. 2nd U.S. Congr. Appl. Mech., Amer. Soc. Mech. Engrs, 623. 


NOTE ADDED AT PROOF STAGE 


In the course of a recent correspondence with Prof. C.-S. Yih, light has 
been thrown on the possible reasons for the discrepancy between the present 
results and those published by him in 1954, to which frequent reference 
has been made in the text above. He has also kindly communicated some 
valuable comments on general aspects of the stability problem. It seems 
desirable that the main points to have emerged from the correspondence 
should be put on record here, particularly since this may obviate the con- 
fusing position which the differences between our two papers may other- 
wise have created. ‘The points listed for convenience as follows are all duc 
essentially to Prof. Yih, although the wording is chiefly mine. 

(1) The direct numerical method used by Yih (1954) in calculating the 
neutral stability conditions for a vertical film was based on only two terms 
of a certain series expansion, and was rather sensitive to small computa- 
tional errors. ‘Thus, although his results were an effective demonstration 
of the low Reynolds numbers at which instability can be expected, they 
cannot be relied upon as more than a rough indication of the range of 
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unstable conditions ; in particular, they cannot be held to confute the 
present precise results concerning instability near the origin of the (R, ~)- 
diagram. 


(2) In unpublished work, Yih has obtained a proof that, for a vertical 
film, disturbances with sufficiently large wave number are stable irrespec- 
tively of surface tension or Reynolds number. ‘This conclusion is perfectly 
in accord with the present work, although beyond its scope since the nature 
of the approximations introduced here preclude an accurate account of 
very short waves. 


(3) The case of zero wave number (i.e. indefinitely long waves) is some- 
what ambiguous from the mathematical point of view. For the long-wave 
approximation used in this paper, the wave number is taken to be small 
enough for its square to be negligible but remains a finite quantity through- 
out the calculations. This approach, which seems the most reasonable 
physically, reveals instability in a vertical film even for vanishingly small 
wave numbers. On the other hand, the long-wave case might be ap- 
proached by setting the wave number equal to zero from the start. ‘This 
step corresponds to a strictly uni-directional disturbance, with the free 
surface entirely undisturbed and undisplaced. In fact, the disturbance 
need not be infinitesimal, since the Navier-Stokes equations can be exactly 
solved. Yih has observed that the latter approach leads to a result indicat- 
ing stability. ‘This result appears to indicate that there is a higher mode 
which is always damped at zero wave number. 


(+) The distinctive features of this problem, notably the instability of a 
vertical film at all Reynolds numbers, all stem from the fact that the energy 
supply to the flow is from gravity potential—in contrast to, say, the energy 
supply from the free stream in the case of boundary layers. If gravity were 
absent, any motion of the film would eventually be damped out ; in this 
sense the film would be stable, although, of course, any form of static dis- 
turbance is neutral when surface tension is also absent. 


(5) From the mathematical point of view, a distinctive feature of the 
problem is that the complex wave velocity ¢ appears not only in the differ- 
ential equations but also in the boundary conditions. ‘This fact, in addition 
to those presented in the last three items, points to the possibility that for 
each pair of values of R and « there may be two values of c—one for the 
primary mode and one for the higher mode. ‘There are no numerical 
results for this possible higher mode, which in any case is not as 


significant physically as the primary mode. 











A note on the application of the supersonic area rule to 
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SUMMARY 
A proof is given that, in spite of certain singularities in the 
oblique area distributions, the supersonic area rule can be used to 
determine correctly the wave drag, according to linearized theory, 
of thin symmetrical rectangular wings at zero incidence. ‘The proof 
could be extended to cover the case of swept wings with supersonic 
edges, for which similar singularities also exist. 


1. INTRODUCTION 

It has recently been suggested that the ‘supersonic area rule’ (Jones 
1953), when used to determine the wave drag at zero lift of wings with 
straight supersonic edges, gives results which disagree with those of con- 
ventional thin wing theory, due to certain singularities which appear in 
the oblique area distributions. It is the purpose of the present note to 
show that this is not the case, provided that these singularities are handled 
with care. A proof is given of the equivalence of the wave drag as calculated 
by the two methods for unswept wings of rectangular planform; the proof 
could easily be extended to cover the case of untapered swept wings with 


supersonic leading and trailing edges. 


2. ‘THE SUPERSONIC AREA RULE 
We consider a rectangular wing of unit chord and span 6 = 2s, and use 
standard rectangular axes with origin O at the centre of the leading edge, 
Ox in the direction of the free stream, Oy along the leading edge and Os 
vertically upward; we shall also use cylindrical coordinates r = \/(v? + s*) 
and 6=tan's/y. ‘The wing section, which must be symmetrical with 
sharp leading and trailing edges, is taken to be 


s=+f(x) (0<x <1), (1) 


and we define fl v) to be zero outside this range. 

The area rule for calculating the wave drag, according to the linearized 
theory of supersonic flow, of any body which can be represented solely by 
a distribution of sources, may be deduced from considerations of momentum 
flux across the large cylinder r = R; use is made of Hayes’s theorem on 
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equivalent source positions (see, for example, Heaslet, Lomax & Spreiter 
1952), which states that, for closed source distributions at large values 
of y and for a given value of 6, the flow due to a unit source is independent 
of its position on the plane 


v = x) +A(ycosé+ssin8§), (2) 


where xy is constant and 8 = \/(M*—1). When applied to plane wings, 
for which the source strength is proportional to the local slope of the surface 
in the x-direction, the theorem shows that the planar source distribution 
may be concentrated onto the x-axis so that, again for a given value of 4, 
the effect at large distances is the same as that of an ‘equivalent body of 
revolution’ whose cross-sectional area at x = € is S(€,@), defined as the area 
of the projection on a plane normal to the axis Ow of the section of the wing 
by the plane. 

v = €+ Bycos#. (3) 


In this way the following well-known formula for the wave drag D may 
be obtained: 


D a } ‘| dé  (E.ONS” } 
int 5 2 a dt i dédy log E—y S'(E, 4)S"(y, 8), (4) 
where py and U are the density and velocity of the tree stream, S(&, 4) ts 
defined above, and primes denote differentiation with respect to € (or 7). 
For symmetrical unyawed wings this reduces to 
D [ 7+ »(( oe n me MYON. B a 
Tp t “Wee 3 |, dt | | dédy log ¢ 1) S (AL yS (7, 8). (5) 
In order that equations (4) and (5) shall be valid, it is necessary that 
S'(€,@) shall be a continuous function of €, though S”’(é,6) may have a 
finite number of discontinuities. But it is evident that, in the case of a 
rectangular wing, this requirement is not satisfied when 6 = Sa; S'(é, 37 
is in fact proportional to f’(€), with discontinuities at € = 0 and 1. ‘lo put 
this in another way, we see that when the source distribution representing 
the wing is concentrated onto the axis along lines x = € (cf. (3) with 6 = 57) 
parallel to the leading edge, the resulting axial distribution is of finite strength 
at the leading and trailing edges; and such a distribution is inadmissible 
since it leads to infinite velocities on the Mach cones from these points. 
Physically, the situation is due to the fact that (according to unmodified 
linearized theory), the flow in the plane 6 = }7 is always partly two- 
dimensional in character, with leading and trailing edge shocks of first 
order strength, however large R may be; and such a flow can never be 


represented by a continuous distribution of purely axial sources. ‘The 
effect is however concentrated, so far as the control surface r = R is 
concerned, in the immediate neighbourhood of the points (8R,0, + R) 
and (8R+1,0, +R) (when R is sufficiently large); and such isolated 
singularities do not invalidate the use of equations (4) and (5) for the overall 
wave drag, provided that care is taken near 6 = + 57. 


: 











i ati et 


Application of the supersonic area rule to rectangular wings 577 


3. DETERMINATION OF OBLIQUE AREA DISTRIBUTIONS 
It is convenient to write 
o(@) = Bscosé. 
‘There are two principal cases to be considered, according as o 2 }. 
(a)o > J or secé 28s (see figure 1a). Since the sections in question 
are simply oblique sections of a cylinder, it is easily seen that 











a 2secO pete | : 
S(€,0) = B [ f(x) dx (-o <é < 1-v) (6a) 
2sec6 /! 
== | f(x) dx (-o <é<®@) (7a) 
b J0 
2secé@ /} : , 
=— | f(w)dx (0 <é <1 +0). (Sa) 
2) oa 
(6) 0 — «6 — sor secé 2fs (see figure 16). In this case the sections 
for o ~ € < 1—o cut both wing tips, and the results have to be 
modified as follows: 
/ 2sec6 site | ; 
S(E,0) =—— |_ f(x) da (-0 <é <o) (6b) 
1) te 
2secO e+ 
-——— | f(x)dx (a © E < 1-0) (7b) 
2sec fl : 
—— |  f(x)dx  (l-o CE <1 +0). (8b) 
2 























£ 
a) 7>72 b) o< 3 


Figure 1. Determination of oblique area distributions. S is the projection of the 
shaded area on a plane perpendicular to Ox. 


Ditterentiating (6) to (8) with respect to €, we find that in both cases 


S'(E, 0) = 2B! sec 0 f(E+0)—f(E-c)! | (9) 


SE, 0) = ZB eee OF (E +oa)—fi'(€ —@)} ; (10) 
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remembering that f(x) = f(x) = 0 for x 0 and a 1. When 6 = !n, 
it is evident that 
S’(€, dar) = 2s f'(E). (11) 
‘The function S’(€,6) is sketched for some typical cases of a biconvex 
section in figure 2. It is clear that, except when 6 = $7, S(€,@) is 
everywhere continuous as required in (5), though S"(&,@) has four 
discontinuities. 


Sie ' > to cats 
‘i S (26) =a255 (2+c)-f(€-o) 
© | 3 Ws Vi if 
5 (2) 
| 
- + | ley & iro | +o 
-C D < Ss 
oe b) o <2 

Figure 2. The function S’(&, @). 


+. EVALUATION OF THE WAVE DRAG 


Substituting from (10) in (5), we obtain 


D 4 i : 
5 | sec*4 di | | dédy log E-—y if (E+e0)f'(y+o)+4 


t f'(E—a)f'(y-—0)-f'(E+oa)f'(n-—0)—-f'(E-—o2)f'(y4+a)}. (12) 


Ikach of the four double integrals with respect to € and 7 may be simplified 
by linear transformations of the tvpe 


leading to the result 


D 4 ; | ms ee: | 
re =sps. | sec’ dt | | dxd\ f (v)/ (y)log |1 
Po , J ) | 


l 6°b* cos?é | 


hs ») 


(x—y)° 


When the function f(x) is given, the double integral with respect to v and \ 
can be evaluated directly, and it is found that the resulting function of @ is 
integrable over the range 0 to }7 and vields the correct result for the wave 
drag D. But in order to obtain a general result it is clearly necessary to 
perform the integration with respect to 4 first; this must however be don 
with care, since the function 


87h? cos?4 


F(0, x, v) sec"4 log! 1 —_——_— 

: (x—y)* 
although regular as 6 —- iz if a yv and having only logarithmic singu- 
larities if @ 7, has an essential singularity when 6 = }7, x=. The 


triple integral (13) for D cannot therefore be expressed as a volume integral, 
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and the order of integration cannot be inverted in the usual way; but this 
difficulty may be overcome as follows. 
An alternative expression for D is 














D S pr. ee ee Bb cos?@ 
tp, U2 = 7B |. sec?@ dé [, } (x) dx I. f G )log 1- “(w—y)> dy, (14) 
and this may be written in the form 
D a . 
Gp ~ Hip 0+ as 
where 
“\m -| rc 6 2h2 . 524 
I,(«) = | sec" dd | f'(x) dx f'(y)log}1 - chal dy, (16) 
/0 Jo. Jo (x—y)? 
and 
“ln am | re } | B2h2 . 524 | 
1,(€) | sec?0d0 | f(x) dx f'(y)log}1 ili | dy. (17) 
me gh OP degre (wy |“ 





‘Vhe singularity at 6 = $7, x = y, has been isolated from the region of 
integration in /,, so that the integration with respect to @ may be performed 
first. 

Using the result (see Appendix) 


¥ F(6,x, y)d6 = -—7 (e v-y Bb), 


} 

| 

J0 } 
f (18) 


we obtain 


| ve 
L(e) T | f(x) dx | f'(y)dy + 
J0 “0 


2 F “Bh Jf B22 ) ses 
ree g f (x) dx I, f(y) 4 re yt ® : 


the second integral occurs only if Bb < 1. 
Now 
“| rr—e vy 
| f(x)dv | f’'O)dy = | f'Qdf(v-@) dv = Of), 
~ 0 “0 “0 
since | f (x)f(x) dx = 0. 
“0 
Hence 
-} . “1-Bb ‘ B2h2 
I(e) = T | f(x) dx | rv, | — ——,} dy+ Ole). (19) 
J Bb /0 ; =F} * 
In evaluating /, we must retain the original order of integration. | Now 
x 7 B2h2 . 524 . “7 | Behe . $74 | 
| f’(y)log|l — =" | dy = f'(x){1 + O(e)} | — log/1- Earl? 








(x—y)? 
1 + (8b/e)cos 4| | 
{1 —(fb €)cos 8 | r 


Bo 
l —ao cos*@ t 
e2 





Bb cos 6 log 


f'(x){1+ O(e)} 4 €log 














4 
; 
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1 


Ife) = (1+ O(e)j | sec?0db | Uf'(x)}* x 











~0 “0 
3? b> | 1 +(Bb/e)cos 4| ) 
; d¢ log|1 — cos*4| + Bb cos loz aie a | ax 
{ €- ] (Bb/e)cos 4| f 
1 
{1+ Ole)! | tf '(x)}* dx x 
~ O 
“hm ( | Bb? + (Bh e)cos 6 | 
x sec?4 < « log|1 — —— cos*| + Bb cos 6 log! = —| } do. 
e | € |\1 —(8b/e)cos 6 | 
Using the results proved in the Appendix 
rim Q2h2 | 
| sec3d log! 1 5 cos*$/d6 = —7 (e pb), 
“0 | 
"ja + (Pb/e)cos b| 
sec 4 log|— | ¢ Lar? € 3b), 
J, ' 1 (2b €)cos 6] ( is 
. “1 . 
we find that I,(e) = \2®Bb |  f'(x) 2 dv {1+ O(e)}. (20) 
0 


Substituting from equations (19) and (20) in (15), and taking the limit 








as « » 0, we obtain finally 
D 4 /l 
Se — ae \ f! 2 0 1: 
D UB Ju (w)}* da (Bb = 1), (21a) 
“] , 
B if (x)? dx 
2J0 


8 geet pon Fo  ela Oe ee 
Bap, | Rh aah Po. F "fs yet? (£6 ). (21 b) 

I;quation (21a) is simply the well-known result that the drag coefficient 
of a rectangular wing is equal to the corresponding two-dimensional value, 
being unaffected by the presence of the wing tips, provided that the aspect 
ratio is greater than 8 -'; equation (21 b) contains the additional term that 
has to be added when the aspect ratio is less than £~', and it has been verified 
that for a biconvex section this gives results in agreement with those of 
Harmon (1947). 

5. CONCLUSIONS 

‘The argument given above shows that the area rule can be applied 
successfully to determine the wave drag of a rectangular wing, in spite of 
the singular area distribution which occurs when 6 = !z, the meaning of 
which has been discussed in $2. In fact the complications introduced by 
this singularity in the anlvsis of §+ are due to the necessity of inverting the 
order of integration in (13) in order to obtain a general proof; if the 
integrations with respect to » and v had been performed first no such 


difficulty would have been experienced, though it should be mentioned 
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that, in view of the simplicity of the final result, the labour involved is quite 
considerable, even for a simple biconvex section. 

It is evident that the method of proof given in the present note could 
easily be extended to cover the case of untapered swept wings. Again 
singular area distribution would occur when @ = cos-'(8-' tan A), where A 
is the angle of sweep, but again they would not lead to serious difficulties. 


APPENDIX. EVALUATION OF TWO DEFINITE INTEGRALS 
rin 


(1) Ha) = sec*4 log 1—acos?6| d@ (a > 0). 





“0 
Write tan@ = ¢: then 
* 2 


J} = log 


~ 0 


t??7+1-—a 








(a) If a < 1, write 1—a=c?; then 
oes f2 i Cc 
} = |, log (5) dt 
?? 4 c t oc 
= | tlog(— + 2ctan-!- —2tan—'t 
t?+ 1 c 








t=0 
=(e~1)n = —a{1- (1 -a)} 
(6) If a > 1, write a—1 = d?; then 
; £ Pa 
J - |. log 7a? | a 





t? — d? lt+d| o 
= ee , | ey ‘ 1 
| toe (+1) dlogi, “al 2 tan | ; 





1+acos@ 


: oe cos U 
(ii) (a) | sec logl; ppeer 


dé (a > 1) 








= | sec #log'1+acos 6) dé. 
0 


We shall first show that ](a) is independent of a, if a > 1. For 


] = 0 
Ka) = { -—— 


» (acos#— 1) 


da 


where the Cauchy principal value must be taken. Writing tan }@ = 7, we 
tind that 
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‘naa 1—cosé 

Hence (a) I L) = | sec @log ion) au 
15 1+cosé 


| logt 
—=—— dr where 7 = tan 36, 





r—>1 


‘The work described above was carried out in the Aerodynamics Division 
of the National Physical Laboratory and is published by permission of the 
Director of the Laboratory. 
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buoyant fluid 
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SUMMARY 
Isolated masses of buoyant fluid were released in a water 
tank. ‘Their width, 2r, and the distance travelled, z, were measured 
as functions of time and were found to follow roughly the laws 


== NB, w = C(gBr)!2, 


where w is the vertical velocity, B the mean buoyancy, and ” and C 
are constants. ‘These equations are predicted by dimensional 
analysis, assuming viscosity to be negligible, and the constants 
appear to be independent of the Reynolds number. It is found 
that C = 1-2, and n is in the neighbourhood of 4. 

Since the Froude number relating the buoyancy and inertia 
forces is the same as for isolated masses of buoyant air in the 
atmosphere, it is concluded that the constants will have the same 
value in this latter case. This is confirmed roughly by observation 
of cumulus cloud towers. 

Some of the characteristics of the motion observed in the 
experiments are described and comparison is made with vortex 
rings. 


1. ‘THE NATURE OF THE PROBLEM 

It was desired to study motion produced in a fluid by buoyancy forces 
alone far from solid boundaries, in order to simulate buoyant convection 
in the atmosphere. From the behaviour of cumulus clouds and_ the 
experience of glider pilots over many years, it appears that convection, 
by which in this paper we mean convection produced by buoyancy forces, 
consists largely of more or less isolated masses of buoyant air rising into 
and mixing with their surroundings from intermittent sources on the ground 
or within clouds (see, for example, Ludlam & Scorer 1953; Welch, Welch & 
Irving 1955). 

For simplicity in the experiments, masses of heavy solutions were 
released from rest at the top of a water tank, and their growth downwards 
was observed photographically. If the volume at release is sufficiently 
small compared with its subsequent volume and no impulse is applied, 
the subsequent growth can be expected to be largely independent of the 
initial configuration of the substance released. In such case the factors 
determining the motion are so few that a very simple dimensional analysis 
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can be applied. ‘(he method has been discussed by Batchelor (1954), 
and experiments of this kind have been carried out by Morton, ‘Taylor & 
‘Turner (1956), and Scorer & Ronne (1956). ‘lhe motion is so complicated 
that no way was seen whereby the simplest dimensional analysis could be 


improved upon, and so the numerical coefficients relating a bulk dimension 
of the cloud of buoyant fluid to its buoyancy and vertical velocity were 
measured, with a view to applying them to the atmosphere. 


2. ‘THE CHARACTERISTICS OF ‘ THERMALS’ 

Following the usage of glider pilots, the buoyant masses of fluid are 
called ‘thermals’ herein. Over the advancing front of the thermal, mixing 
with the fluid ahead takes place. Inside and at the rear the fluid circulates 
rather as in a vortex ring, the motion being relatively smooth. As a result 
of the mixing the volume continuously grows, and, although the vertical 


MOTION RELATIVE | MOTION RELATIVE 
TO GROWING THERMAL | TO A STATIONARY 
tt OBSERVER 


Figure 4. Approximate distribution of velocity in a typical thermal. The black area 
illustrates successive positions of a portion of the thermal, into which mixing 
has just occurred, beginning with the square. 


momentum steadily increases through action of the buoyancy force, the 
velocity decreases on account of the incorporation of exterior fluid. 
Examples are seen in figures 1, 2 and 3 (plates 1 and 2). The pattern of 
motion is represented in figure 4. An idea of the exterior motion is gained 
from figure 3 (plate 2). 

The thermals were released from a thin hemispherical copper cup 
(seen in figure 3, plate 2) pivoted about a horizontal axis, the level of the 
fluid inside being made the same as outside before release so as to avoid any 
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Figure 1. Series of pictures taken from cine film showing the growth of thermals in 
a water tank, made visible by white precipitate. (Experimental numbers 
106, 107, 111). No. 107 is unusually well formed. 


Figure 2. Pictures illustrating how a shape assumed early in the life of a thermal may persist while the 


volume increases several times. This mav be contrasted with no. 111 of figure 1 in which 
the protuberances are shed. Note the * cauliflower’ nature of the surface. No shearing 


1 


iDerances 


motion is visible in these prot 





Figure 


. Scorer, Experime on convection of isolated masses of buoyant fluid, Plate 2 





3. In each picture there are two exposures about sec apart. ‘The first picture shows the free 


rate of sink of some lowly sinking particles In the second and _ third pictures an idea of 

velocity field outside the thermal can be gained. ‘The particles immediately in front of the thern 
are displaced downwards as much as the faster sinking ones lower down, while those at the side ar 
in an upcurrent. ‘Those that have disappeared have entered the advancing *‘ front’ of the thern 


not the rear. ‘Che markings on the vertical rod are at 10 cm intervals. 


Figure 9. View of a thermal from above showing the hollowed-out rear. 
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initial impulse. ‘The cup was turned over quickly by hand, and a negligible 
amount of motion was produced thereby. A dense white precipitate was 
used to make the buoyant fluid visible. Immediately after release the front 
surface of the thermal became covered with protuberances and the volume 
began to increase. After the main mass escaped from the surface there 
was often a stem left behind. But this was very dilute and had negligible 
velocity in it, and was therefore assumed to contain a negligible fraction 
of the original material. ‘This stem is thought to be a kind of ‘splash’ 
produced as the original material becomes rearranged into the configuration 
which it ultimately assumes. 

The motion consists initially of an acceleration. After the thermal 
has travelled about 1-5 diameters deceleration begins, and the measurements 
were concerned with the subsequent motion. 

In order that simple dimensional analysis should be applicable it is 
necessary to suppose that the thermal behaves as it if originated at a virtual 
point origin. ‘Therefore, in order to allow location of this origin as closely 
as possible, the same cup (radius 3 in.) was used in nearly all the experiments 
and was filled to approximately the same depth each time. ‘Thus at least 
the geometry of the initial conditions was approximately the same, and the 
origin was assumed to be in the same position in all these cases. Some other 
cases were included when they were susceptible of easy measurement, 
that is, when the thermals grew with reasonable symmetry. 

The density differences were kept small, and were initially less than 
15°,. If this is not done two difficulties arise. First, the motion may be 
different according to whether the heavier fluid is inside or outside (downward 
or upward motion), and, second, it may vary according to the magnitude 
of the density differences. During the measured part of the thermals’ 
lives the density differences were mostly much less than 5°, and were 
reduced to around 0-1°,, towards the end of each experiment. 

Some thermals grew in a grossly asymmetrical manner. ‘This is thought 
to be due to internal motions of the thermal, or occasionally in the tank, 
at the moment of release. ‘These cases were rejected. 

Occasionally odd characteristics of shape seen in the early stages were 
retained almost throughout the observed life of a thermal (see figure 2, 
plate 1). In particular the angle subtended at the point of origin tended 
to remain the same, though it varied somewhat from thermal to thermal. 


3. SIMPLE DIMENSIONAL ANALYSIS 
‘The only factors which can determine the velocities are the buoyancy 
forces and the size. ‘The local buoyancy force on unit mass of fluid is 


fe = gB, (1) 


5 


Po 
p; and py being the interior and exterior densities. ‘The density ratio B 
varies within the thermal because the dilution varies, and so the mean value 
B will be used, the distribution of B being at present unknown. If r is 


F.M. 2R 
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the radius of the largest horizontal section, the vertical velocity of the front 
of the thermal must then be of the form 


w* = C2gBr (2) 


from dimensional considerations. C? is a Froude number, being a ratio 
between inertia forces and buoyancy forces. 

The same formula can likewise be derived for the rise of a bubble of 
buoyant fluid which will not mix the surroundings. One example of this 
is the formula for the rise of bubbles of air through water, derived by 
Davies & ‘Taylor (1950). There are some interesting differences. The 
buoyancy B in that case is virtually unity, but it is also constant, and for 
this reason it is not necessary to assume that B is small in the case of immiscible 
fluids. ‘The drag is here due to the formation of a wake, so that the buoyant 
energy is converted into energy of turbulent motion of the surroundings ; 
in our case there is no wake, but the ‘ drag’ is due to the continuous incorpora- 
tion of outside fluid. 

If z is the height of the foremost part of the thermal the two velocities 
dz/dt and dr/dt are proportional to all other velocities, and so 


z= nmr (3) 


if the origin of z is suitably chosen. ‘This, the virtual origin, is at the apex 
of the cone swept out by the largest horizontal section. mn has to be 
determined by experiment, and, like C, must be the same for all thermals, 
according to the assumptions made. Since only one representative velocity 
exists (which is equivalent to saying that there is no Froude number other 
than C?), the motion is similar at all stages in the life of a thermal and all 
thermals are similar. 

The volume J’ of the thermal, which is seen in the experiments to possess 
a clearly defined outline, is given by 


V = mr, (+) 


where m is another number to be determined by experiment. The total 
buoyancy of the thermal is constant, so that if suffix 0 denotes the values at 
the moment of release, when p; is uniform, 


gBmr3 = 9B, V9. (>) 


In the experiments B, and IV’, were measured before the release of the thermal. 
Writing w = dz dt and integrating (2), we obtain 


kz? =t, (6) 
where k = m'2/2nC (2B, V,)*2, (7) 


the origin of ¢ being suitably chosen. 

The objective in the experiments was first to confirm the relation (6) 
for each thermal separately, and then (7) for a variety of values of gB, TV, 
the total weight (deficiency or excess) of the thermal. 





: 
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4, "THE EFFECT OF VISCOSITY 

In previous work on experiments and analysis of this kind (Batchelor 
1954; Morton, Taylor & Turner 1956), the motion has been described as 
‘fully turbulent’, and it has been assumed that the viscous forces are 
negligible in comparison with the inertia forces (including the eddy stresses) 
and the buoyancy forces. In the case of a conical plume from a maintained 
source of buoyancy, the fact that the plume is conical confirms the correctness 
of this assumption because the Reynolds number varies up the plume. 
Furthermore, the motion is quite evidently turbulent throughout the volume 
of the plume. However, when small drops of coloured fluid are let fall in 
clear fluid they do not develop turbulent motion but retain a smooth outline. 
In anindividual thermal from a source which is not maintained, the Reynolds 
number wr/v is equal to 1/2kv, which is the same throughout the life of a 
thermal; so that it would in any case be expected to grow along a cone and 
retain a constant shape even if viscous forces were important. But the 
Reynolds number varies from one thermal to another, and so if (7) is 
confirmed we can assume that variations in Reynolds number make no 
difference and that viscous forces are negligible. 

In the case of thermals in cumulus clouds (see §7 below) the Reynolds 
number is about one thousand times greater than in these experiments. 
The constant C is roughly the same, and so the assumption that viscosity 
can be neglected is further justified. 

5, EXPERIMENTAL RESULTS 

‘The thermals were released into a tank with horizontal section 2 ft. x 4 ft. 
The thermals traversed the depth of 35 ft. in between 5 and 30 seconds 
according to the weight excess. 

The virtual origin was determined roughly by plotting r against sz, 
extrapolating backwards, and taking z= 0 where r=0. Variations in 
the position of the virtual origin by 3-4 cm had a negligible effect on the 
value of k obtained when s* was plotted against ¢, and so the origin was 
taken to have the same position relative to the bottom point of the cup in 
all cases in which Vy was nearly thesame. Ina few cases Vy was substantially 
different and for those the origin was determined individually. 

The crudeness with which the measurements have to be made must 
be mentioned. Successive outlines of four typical thermals are shown in 
figure 5. With them are the plots of 2? against ¢ and the straight lines drawn 
to determine k. As expected, the points deviate from the line in the early 
part of the life while the thermal grows into its ultimate shape. Some of 
the thermals grew somewhat asymmetrically or irregularly but nevertheless 
yielded fairly definite values of k. ‘The results are given in table 1. 

No qualitative difference in the motion was evident when the motion 
was upwards. A fluorescent dye was used to make the thermal of methy! 
alcohol visible, but the outline was less sharp than with the white precipitate. 
Only one upward moving case (no. 201) was measured, and this agrees 
closely with the others. 


2R2 
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Figure 5. Successive outlines of thermals traced from photographs. Beneath each 


is a graph of =? against t. (a) and (6) show fairly typical thermals; 106 (c) 
is an example of a very large value of n; 113 (d) shows a very asymmetrical 


thermal. 


In figure 6 we have plotted (V5 By)!? against k-', and the straight line 
was drawn by eye through the origin to fit the points. ‘There is no evident 
systematic deviation from this line and so viscosity can reasonably be 
assumed to be negligible. The straight line represents the relation 
k-1 = 180(V,, B,)"?, i.e. (see (7)) 

2nC'[ne* = 1G0e"* = 5-7. (8) 
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It was observed that ” and m were not the same in all thermals, rough 
mean values being 
m= 3, we 4. (9) 
‘The variations in m did not appear to be related to other properties of the 
thermals, but some systematic variation in m was detected (see below). 
With these values we find 























C = 1-2. (10) 
stiisian.| we es | 
Experiment | (I 0 Bo)" | “f (Bo) 1/2p-1| Z (Vp Bo)-!2k-n?2 
no. cm cm* sec 
] | 

25 6-63 1140 170 3°8 87 
110 6:53 1390 210 4-8 96 
107 6-31 1090 170 3°7 89 
112 5:81 1320 220 4:8 100 
104 5-49 1010 200 4-1 99 
102 4:97 1000 200 4-1 99 
109 4-54 880 190 4-2 93 
24 4-45 850 190 3-6 100 
72 3-86 625 160 3:1 91 
111 3°45 568 160 3-0 93 
30 2°82 490 170 3-3 93 
106 — 477 — | 50 | — 
108 2-44 406 170 3-5 | 91 
118 1-94 357 | 190 3-6 100 
201 147 | 250 | 170 | 48 | 78 
117 1:37 250 180 | 4-0 90 
113 1-26 200 =| ~~ 150 2-9 | 88 
105 1-01 172 170 29 | 100 








Table 1. Measurements of total mass (deficiency or excess), k (giving the velocity 
as a function of height), and rate of widening 7. 


6. VARIATIONS IN THE ANGLE OF THE ENVELOPING CONE 
We may attempt to take account of the variations of n in the following 
way. Ifthe volume is expressed in terms of rz, and suffix 1 is used to denote 
values at a given time, then 


Bite co 

Be®s = By 9; 2, (11) 
In the place of (2), we have 
wo = Cn“! gB2z)!? = Cr'*(2B, 7; z,)!2="1. (12) 
Accordingly instead of having C as a constant we might have Cn'? as a 
constant, for now gB, 172, represents the weight (deficiency or excess) of the 
thermal. In that case 

4on'2( 0B, Vo) *k 1 n-? = Cn"? = const. (13) 


To test this we have plotted (V, By)~!?k-! against n in figure 7, and the 
curve drawn is part of the parabola represented by 
n = 0-00012(V, By) 1k. (14) 
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If the points lay on this curve it would imply that Cn!2 was in fact constant. 
This appears to take some account of the variations in , because over this 
series of experiments Cn! ? varies by only half as much as C itself (see table 1) ; 
we have no indication of the mechanism, except that if the thermal grows 
along a wider cone it travels more slowly at a given distance from the virtual 
origin, presumably because it has been more diluted. If m = 3, equations 
(13) and (14) give 











Ce? = 25 (15) 
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Figure 6. Figure 7. 


Figure 6. (V9 Bo)!” plotted against k~! for the thermals enumerated in table 1. The 
straight line represents k~! 180 (V> By)}. 
Figure 7. The relationship n = 0-00012/V» By k?, indicated by the curve drawn, is an 
improvement upon the crude assumption that » is the same for all experiments. 


The tendency of a thermal to keep the same value of m throughout its 
life may be compared with the behaviour of the air bubbles of Davies and 
‘Taylor whose angular ‘aperture’ varied from one bubble to another. 

In general, the wider thermals (smaller n) tended to leave more ‘ debris’ 
behind and this may possibly account for their slower motion. The 
unusually slow one (113, shown in figure 5 (d)) developed very asymmetrically 
and only part of it was in effect measured, which would imply that its 
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total buoyancy was overestimated. ‘The point derived from it lies below 


the line in figure 6. 


7. COMPARISON WITH ATMOSPHERIC THERMALS 


The only published measurements of thermals in the atmosphere 
available for comparison are those given by Malkus & Scorer (1955). ‘They 
measured the rate of rise of isolated cumulus towers, and found that early 
in their life on emerging from the parent cloud the relationship 

w? = -9Br (16) 
was fairly well obeyed. ‘The cloud towers ceased to rise after ascending 
about 1:5 diameters, mainly because the buoyancy decreased to zero, or 
even changed sign, as the cloud droplets evaporated when the thermal was 
mixed with the surrounding dry air. In 1-5 diameters the whole of the 
interior of the thermal would have been exposed to mixing. Inside the 
parent cloud there is no such evaporation, so that soon after emerging the 
motion would have been little affected by evaporation in the exterior shell 
of the thermal. 

In the light of the present experiments and the nature of the mixing 
process it is probable that their estimates of buoyancy were too high; 
hence the value C = 89 obtained for clouds by comparing (2) and (16) 
would be expected to be a little below the true value. ‘This agrees satis- 
factorily with the value of 1-2 obtained from the present experiments, and 
is certainly within experimental error. Of course not all differences should 
be attributed to such errors. It should be pointed out that the cloud 
thermals were ascending into stably stratified surroundings and were not 
increasing in size as they ascended, so that the terminal velocity may have 
differed on this account; but the difference cannot have been very great 
because the thermals were decelerating. ‘The measurement of r in the case 
of the clouds was made soon after emergence from the parent cloud. 

Recently Ludlam & Saunders (1956) have made measurements of the 
vertical velocity of cumulus towers just as they began to emerge as thermals 
from the parent cloud. ‘They did not attempt to measure the diameter but, 
having (12) in mind, expressed the size simply as distance above the ground. 
It cannot be much in error to assume that the virtual source is at the ground. 
Consequently it was possible to make an estimate of the ‘constant’ Cn~}? 
in (12), and the mean value they obtained was roughly 0-5. ‘This was based 
on the behaviour of the most rapidly rising thermals at any moment and 
not on the motion of many slower ones which could have originated at 
higher levels with unknown z. The closeness of agreement with the value 
of 0-6 obtained in our experiments is certainly fortuitous in view of the 
roughness of the various estimates that are necessary, but it encourages the 
belief that no new dynamical process is involved even though the linear 
scale is much greater. ‘Their calculations are not given in the paper referred 
to, other aspects of the observations being discussed there, but they will 
be submitted for publication in due course. 
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An attempt was made during the summer of 1956 by Miss Betsy 
Woodward to make the appropriate measurements in clear air thermals 
below cloud base in a Skylark II glider. The humidity and temperature 
inside and outside of a thermal were measured using a wet and dry thermistor 
out-of-balance resistance thermometer, the readings being spoken into_ 
a wire recorder during flight. Spot readings had to be used to determine 
mean buoyancy of the thermal. ‘The size of the thermal was estimated very 
roughly from the size of the turning circle which she could use during the 
ascent. ‘The vertical velocity is not easily measured because it varies from 
one part of the thermal to another. It can be properly measured only if 
the glider has risen through the thermal to near the top, where, because of 
the decrease of vertical velocity and the outward radial component (which 
is equivalent to a downward velocity for a circling glider) the glider can rise 
only as fast as the whole thermal. 

Using the observed buoyancies, etc. (Woodward 1956), a vertical velocity 
of 2 m/sec was deduced using a value of 1-2 for C and r = 160 m, B = = 
but since the uncertainties are so many it is neither surprising nor gratifying 
that the observed rate of rise of the glider when near the top of the thermal 
was only 1 m/sec. It is not considered profitable to attempt to explain 
away this discrepancy because there are too many different ways in which 
to do it easily. ‘lhe observation is quoted because it shows that the model 
experiments again lead to an estimate of w of the correct order of magnitude. 


8. ‘THE INTERIOR MOTION OF A THERMAL 

Miss Woodward has made the interior motion visible by releasing a 
transparent thermal containing white pellets about 5mm in diameter 
balanced with pieces of fine wire attached to them so that they had a very 
small terminal velocity in water. The thermal was then illuminated by 
a flat beam of light about 3 cm thick, and the motion in the central vertical 
section of the thermal was photographed on 16 mm cine film. 

‘The positions of particles within the thermal in successive pictures were 
drawn by steadily reducing the picture of the thermal to a constant size 
on a screen as the motion proceeded. ‘The motion during the whole life 
of the thermal was thus concentrated into one picture, on the assumption 
that the motion was similar at all stages. Isopleths of horizontal and vertical 
component of velocity obtained in this way are shown in figure 8. Speeds 
are expressed as multiples of the forward velocity of the front of the thermal. 
This figure is based on observation of only one thermal. 

It is seen that the velocity is a maximum in the centre, whereas in a vortex 
in which the vorticity is concentrated in a ring the velocity is a maximum 
at the inner surface of the ring and decreases towards the axis. 


9, COMPARISON WITH BUOYANT VORTEX RINGS 


In a recent paper ‘Turner (1957) has discussed the behaviour of buoyant 


fluid ejected in the form of a ring vortex. Several years ago some glider 
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pilots were accustomed to think of the thermals in which they soared as 
possessing a configuration like a ring vortex. We now see that there is 
much justification for this. One important difference is that a thermal 
is not possessed of any circulation, kinetic energy, or forward momentum 
at the moment of release. ‘lhe circulation is generated gradually in the 
early stages of its existence. As Dr Batchelor noted to me in correspondence, 
the circulation round a thermal when it has achieved the final configuration 
is proportional to ws and is therefore constant (see (2), (3) and (5)). In 





Figure 8. The distribution of velocity obtained by observing the motion of particles 
inside a liquid thermal. The values of the horizontal (left hand diagram) 
and vertical (right) velocities are expressed as multiples of the vertical velocity 
of the front of the thermal. 


the early stages there is a hydrostatic buoyancy force acting on a column of 
buoyant fluid up the centre and increasing the circulation. Later on much 
of this column is replaced by exterior fluid, as is seen in figure 9 (plate 2). 
If the rate of change of circulation is zero, then 
$<? -() (17) 
J P 
for acircuit passing through the centre of the thermal and round the outside, 
p’ being the departure of the pressure from the hydrostatic value in the 
exterior fluid. If the fluid is buoyant along any part of the axis of the thermal 
the tendency to produce circulation must be counteracted by departure from 
the hydrostatic pressure due to the motion. Since the motion has been 
produced by the buoyancy forces we should expect the configuration of 
buoyancy to change to a distribution which does not continue to produce 
circulation; this means reducing the depth of buoyant fluid up the axis of 
the thermal. 
In Turner’s vortex rings the circulation remained constant because the 
buoyant fluid did not extend to the centre. The initial impulse was of the 
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order of 1500 gmcm sec"!. The buoyancy force was of the order of 
300 dynes and the impulse of a ring was increased by a factor of between 
2 and 10 during the period of observation. One would suspect that ulti- 
mately its motion would tend towards that of an isolated thermal as the initial 
impulse tended towards a small fraction of the total, but the rings did not 
reach this stage in his experiments. 

‘The only constant characteristics of a thermal are its total weight 


(deficiency or excess), gBV, and the fluid density p. Dimensional 
considerations show that its circulation A must therefore be given by 
K = const(gBV /p)'?. (18) 


It seems that if the circulation is less than this value the buoyancy forces 
create circulation. But if the circulation is much in excess of this value 
mixing towards the axis must be inhibited. ‘The vortex ring could only 
become a thermal if the circulation is reduced to the value appropriate to 
its total buoyancy according to (18). This might happen ultimately on 
account of viscosity, but then the viscous forces would interfere with the 
motion of the thermal. In ‘Turner’s case the total buoyancy of the rings 
was about 0-3 gm wt., which is much less than in the present experiments, 
in Which it ranged from 1 to 40 gm wt., so that they are not immediately 
comparable. 

Among those who have participated in this work Miss Betsy Woodward 
deserves special thanks for her assistance in performing the experiments 
and measuring the photographs. She has also done much pioneer work 
in the investigation of thermals in gliders. Her work is at present supported 
by the Munitalp Foundation. 
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ERRATUM 


The first equation in the Summary should read as s=nr. 
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SUMMARY 

An expression is obtained for the jump in the vorticity across 
a gasdynamic discontinuity in an inviscid flow. ‘This result 
generalizes results of Truesdell (1952) and Lighthill (1957) for the 
vorticity behind a steady curved shock in a uniform flow and that of 
Emmons (1957) for the vorticity jump across a steady flame. ‘The 
derivation is a dynamical one, and no assumptions on the com- 
position or thermodynamic properties of the fluid are made. ‘The 
jump in vorticity in the steady flow case is found to depend upon 
the jump in density and upon gradients along the discontinuity 
surface of the tangential velocity component and of the normal 
mass flow. An analogous result is found with unsteady flow. 


INTRODUCTION 

It is a well-known fact that the vorticity in a steady compressible flow 
is generally discontinuous across a shock wave. ‘Truesdell (1952) first 
obtained a general expression for the vorticity behind a curved two- 
dimensional steady shock in a uniform flow. ‘Truesdell states his result 
thus: the magnitude of the vorticity generated by a shock of given strength 
and curvature depends only on the magnitude of the tangential component of 
velocity and is independent of the form of the equation of state. Here the 
strength of a shock is defined as the density jump across the shock divided 
by the density in front. Other investigators, including this writer and 
Lighthill, unaware of ‘Truesdell’s work, later rederived the result. However, 
Lighthill (1957, pp. 14, 15) also provided a significant generalization in 
showing that the result was valid, when expressed in terms of the axes of 
principal curvature, for a steady shock wave of general shape in a uniform 
flow. Lighthill specified the density ratio across the shock to be the 
limiting value for a very strong shock, but this restriction is not used in 
his analysis and is unnecessary. 

The vorticity relations just discussed depend upon a number of 
assumptions, of which the most restrictive is that the flow in front of the 
shock is uniform. The methods used all require the application of Crocco’s 
vorticity law relating the vorticity with the entropy gradient, and the 
assumptions which underlie this law must be made. ‘These assumptions 
are that the flow is steady, is isocompositional, and is isoenergetic (i.e. 
has constant total enthalpy). 
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A similar question arises in the theory of the propagation of a laminar 
flame in a combustible gas mixture, as to the magnitude of the vorticity 
engendered by the flame. In this case the disturbance caused by the flame 
is felt in the fluid region in front of the flame, and the assumption of uniform 
flow in this region cannot be made. Emmons (1957) has obtained an 
expression for the change in vorticity across a flame in terms of the density 
ratio across the discontinuity and the jump in tangential entropy gradient. 
Emmons uses a modified form of Crocco’s vorticity law in his derivation, 
and his results are also restricted to steady and isocompositional flows. 
Emmons’s results and the results of this paper were obtained independently 
and are not directly comparable. 

The simplicity of ‘Truesdell’s relation and its independence of the 
thermodynamic state of the fluid suggests that the relation is a purely 
dynamical one. A purely dynamical relation must be derivable without 
recourse to a thermodynamic law such as Crocco’s vorticity law. In 
presenting a dynamical derivation for the vorticity jump we shall consider 
first the case of steady flow, in which the principal features of the derivation 
are not obscured by the complexity of an arbitrarily moving discontinuity 
Our method involves a vector decomposition normal to and tangential to 
the discontinuity surface. ‘The discontinuities considered are gasdynamic 
ones, including shock waves, deflagrations, and detonations, but excluding 
contact discontinuities or slip streams. 


‘THE VORTICITY JUMP IN STEADY FLOW 

‘The discontinuity surface is assumed to be sufficiently smooth and 
continuous, without discontinuities in slope. ‘The orientation of the 
surface is specified at any point by the unit normal vector n, which is assumed 
to be differentiable along the surface. A natural local coordinate system 
is imagined, based on any intrinsic two-dimensional coordinate system on 
the surface. ‘The additional coordinate n is distance along a straight normal 
from the discontinuity surface, with the corresponding coordinate surfaces 
geodesic parallels to the discontinuity surface. 

The unit vector n is the gradient of this normal coordinate, and we 
obtain immediately 

Vxan=0. (1) 
The negative of the gradient of n is the two-dimensional curvature tensor 
in the discontinuity surface 
K = —Vn, (2) 

which trom (1) must be symmetric. Because n is a unit vector Vn can 
have no component normal to the surface. The curvature is defined so 
as to be positive if the surface is concave on the side from which n is 
directed. 

In general, all vectors will be considered in a form decomposed, with 
respect to the surface, into a normal component and a tangential 
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two-dimensional vector. ‘Thus, with the velocity vector q for the flow 
field, we have 

q=n9, +4 (3) 
where q, =m.q, q, = nx(q xn). (4) 
The vorticity is the curl of q, and is expressed as 


C= Vx(nq,)+V xq, 


=n(V x q),, +(V x q)),—n x V4q,,; (5) 
we note that t, = (V x@q),. (6) 
Vector expansion of the identity V(m.q,) = 0 gives us the result 

(V xq) xn =n.Vq,+q,. Yn, (7) 
and from (5) we obtain an expression for 

G = nx (3 -a-&-Vian), (5) 


in which the subscript ¢ on the nabla operator V indicates that only the 
tangential part of the derivative is included. 

We need also an expression for the tangential part of the directional 
derivative of the velocity, q.Vq. <A straightforward analysis starting 
with (3) gives 


A 


ie » bes 
(4-Va)= mn (3 - 4) 4-4 (9) 
We now turn to the steady discontinuity relations, with the symbol 6 
indicating the jump of a quantity across the discontinuity. The continuity 
condition is 


5(pq,,) = 9, (10) 
while the tangential momentum condition is 
dq, = 0. (11) 


It is at this point, in establishing (11), that we assume non-zero mass flow 
through the discontinuity and exclude contact discontinuities from our 
consideration. ‘The normal momentum condition is 


—dp = PIn Gn. (12) 
The energy condition across the discontinuity is not used. Since relations 
(10), (11) and (12) hold everywhere on the discontinuity surface we may 
differentiate them along the surface. Also, the operators 6 and V, are 
commutable. ‘The increment in the vorticity is 


86, = 0, (13) 


ot, -nx9( Vig), (14) 


from (6) and (8). Equation (13) is a statement of the obvious fact that 
the normal component of the vorticity is continuous across a discontinuity. 
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We next write the tangential part of the momentum equation for the 
fluid flow. Since the flow is assumed steady, we obtain immediately with 
the aid of (9) 


VP = paa( st aR) + pq Veah (15) 
Our method now is to equate the tangential gradient of (12) with the jump 
in (15); 
— OVD = (09, )V 1 Gu + VlPIn)°In 


= PIn 5( =) 4 -G Vy Gop. (16) 
Combining (10), Ss and (16) gives us our desired expression for 5f,; 
6, = nx [V(pq,, )5(p*) — (pg, Qe» Ve, 8(p)]. (17) 


In this expression the only jumps which appear are jumps in the density, 
and these are multiplied by quantities which are continuous across the 
discontinuity. Because of (13), the subscript t on € in (17) is not needed. 

Finally, we obtain ‘Truesdell’s expression with Lighthill’s generalization, 
taking the flow in front of the discontinuity to be uniform with density p, 
and velocity U. We evaluate the quantities appearing in (17), 


pq, = po U,,, (18a) 
V(pq,,) = pp U. Vn = —poU,. B, (18 b) 
q,-V.q,; =nx(Ux U,.Vn) = U, U,.&, (18c) 
and obtain 
ot, =n x U,. K(p; '— pp ‘dp, (19) 
or, since the flow in front of the discontinuity is irrotational, 
(1 —e)? 


=- 


Here p, is the density behind the discontinuity, and 





nx U,.&. (20) 


€ = Po/P1 (21) 
is the density ratio across the discontinuity. It should be noted that this 
result is not limited to shock waves, but would hold also for a condensation 
shock or detonation. 

"THE VORTICITY JUMP IN UNSTEADY FLOW 

In dealing with the unsteady flow case we must take into account not 
only the unsteady nature of the flow field but also an arbitrary motion of 
the discontinuity surface. Many of the relations obtained for the steady 
flow case are still valid, principally the kinematic relations (6), (8) and (9), 
and the _ relations for the tangential velocity and the vorticity (11), 
(13) and (14). The primary modification which must be made is to the 
ae we a normal mass flow across the discontinuity. 

Of the available ways in which the shock motion may be specified, the 
most convenient is in terms of its normal velocity ng,. We use the notation 
d/dt to denote the time derivative taken at a point which always lies on the 
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discontinuity and which moves in a direction normal to the discontinuity 
when the discontinuity moves. If, is the angular velocity of the surface 
measured at such a point and is defined as a tangential vector, we may show 
that 


m= w, xn = — Vi qe (22) 
We designate the relative normal velocity component q, as 
Ir = In — Ios (23) 
and obtain the discontinuity relations 
5(pq,) = 0 (24) 
in place of (10), and 
— db = (pq,)0q, = (p9,)84,, (25) 


in place of (12). We use the symbol D to designate the classical total or 
material derivative with respect to time; and we introduce the symbol D, 
to designate a tangential total time derivative, the time derivative with 
respect to an observer who moves along the shock with a total velocity 
ng,+q,. If such a derivative is a vector, we shall include only its tangential 
part. Note that for such an observer the discontinuity is a normal one. 
The distinction between the derivatives D and D, is immaterial when they 
are applied to the vector n, because the space derivative of n in the normal 
direction is zero. ‘Thus we have 





1 
Dn = D,n = = +q,. Vn. (26) 
¢ 
For the tangential velocity q, we have 
dq, © oq, 77 
dt ~ ot + % Gn? or 
and e dq; ; 
’ D,q, = (Da): 4 a a (FB) +4 Vi Q- (28) 


In place of (15) we may now write 
—V, p = p(Dngq,,), + p(Day), 


cq, 
= PG, +9s)D, 0+ pq,= + pD,G- (29) 
Applying the same procedure as was used in the case of steady flow, we 
obtain the desired expression for the vorticity jump as 
80, = n x [V(pq,)8(p*) — (pq, '(D, Qr + 9s Dy)0(p)] (30) 
in place of (17). The quantity D,n is given in (26), or, with (2) and (22), 
may be re-expressed as 
Dn = —V,q,—q- B. (31) 
The quantity D,q, is given by (28) with (27). It may be checked that (30) 
is invariant under a velocity transformation to another unaccelerated frame 
of reference. Again, the subscript ¢ on § in (30) is not needed. 
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As an example of our result (30) for unsteady flow we apply it to the case 
in which the flow in front of the discontinuity is at a uniform density p, 
and uniform velocity U. We evaluate the quantities appearing in (30) to 

obtain, in analogy with (18), 
PIr = PAU, — 4s) ( ) 
Vi(pq,) 2 po D, n, (32 b) 
D,q, = nx (Ux Dn) = —U, Dyn. (94C) 

Ihe vorticity jump is then obtained from (30), 
st, = nx (U,.K+V, qi '—py Wp, (33) 

in place of (19), or 


(1-<)? 
C=- : os x (U,.K + V,q,) (34) 





in place of (20), with e givenas before by(21). This result is the generalization 
to unsteady flow of the ‘Truesdell-Lighthill vorticity expression. It may 
be noted that the fundamental quantity appearing in (20) or (34) is the 
angular velocity of the discontinuity surface with respect to the observer 
for whom a derivative is D,. We may rewrite (20) or (34) in the form 


ie (35) 

where the angular velocity w is 
w =nxD,n (36a) 
=w.+U,.K xn, (36b) 


with w, defined in (22). 
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SUMMARY 


It is shown that the interpenetration of two ionized streams is 
arrested, as a rule, not because of individual collisions between 
particles belonging to opposite streams, but because the whole 
system of charged particles is unstable. ‘The smallest wavelength 
of an unstable oscillation is Amin, where 


(7m, U? U2\-3/4 
Amin = J Tir )(1 = “4 ; 


Here + U are the velocities of the undisturbed streams, and N is 
the density of electrons in each. 

A further calculation for the non-relativistic case deals with 
the amplification of the plasma oscillations present in two 
colliding streams. It is shown that these grow rapidly and that 
Torit\V (mp U?/zNe*) is the distance of interpenetration achieved 
before the counterstreaming of the electrons is brought to a 
halt. The value of z¢rit depends only insensitively on the ratio 
of the internal plasma energy densities 7'p; to the kinetic energy 
densities 7in in the streams. For example, terit = 9-0 when 
Tpi: Txin = 1:10, and tert = 19-0 when Tp): Tein = 1: 10°. 


1. INTRODUCTION 

The interpenetration of two fully ionized streams of gas is considered 
in this paper. It is shown that, as a rule, the counterstreaming is stopped 
because of a collective instability among the electrons present, and that 
close encounters between charged particles belonging to the two streams 
are relatively unimportant in destroying the systematic motion. 

Advance information of some of the results in this paper was given by 
the author in an earlier paper (Kahn 1955) where it was suggested that 
this instability effect leads to the conversion of the kinetic energy of the two 
streams into the energy of irregular plasma oscillations. A_ particular 
example of this might be found in the collision of two galaxies—as in the 
radio source Cygnus A. One would therefore expect random electric 
fields to be present in the plasma left after a collision. A fast charged particle 
passing through such a medium would then be expected to radiate electro- 
magnetic waves. 

The instability effect may also be important in the laboratory. It may, 
for instance, explain the smallness of the distance to which two counter- 
streaming masses of plasma can interpenetrate. ‘That the distance is small 
is shown, for instance, by Bostick’s (1956) experiments. 


F.M. 258 
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In the following two sections we derive the condition that an oscillation 
of a particular wavelength may be amplitied. ‘The treatment given is 
non-relativistic in $2, and relativistic in $3. In $4 we consider the 
amplification of the oscillations which are present in the two streams before 
they have collided, and find a closer estimate of the time that elapses before 
the counterstreaming is stopped. 

Some calculations concerning this effect have been given by Lampert 


(1956), who uses the dispersion relations for a plasma. 


2. A NON-RELATIVISTIC TREATMENT 

Consider two interpenetrating streams, each of infinite extent, and each 
consisting of N protons and N electrons per cm*. Let stream 1 have a 
velocity U in the direction of increasing x, and stream 2 a velocity U in 
the direction of decreasing x. Let e« be the charge on the electron and m, 
its rest mass. Let the temperature 7 ineach stream be such that kT < m, U?. 
The thermal motion among the particles may then be neglected. 

An electron passing within a distance r of a charged particle from the 
other stream experiences a change e*/r in potential energy. Its kinetic 
energy relative to this particle is }m(2U)? = 2m,U*. An appreciable 
deflection in the direction of the electron’s motion occurs if 


e/r = O(2m, U*), 
that is, r= Ofe?|Zm, U2). (1) 


The etfective cross-sectional area for such a close collision is thus of the 
order of S = we?/4m-U*. ‘These encounters would eventually destroy the 
systematic motion. ‘The distance of interpenetration will at most be of the 


order of 


Wo = tr: (2) 


But there is a collective instability among the electrons which will stop 
the counterstreaming very much sooner. Assume, for the moment, that 
the protons, because of their much greater mass, continue their uniform 
motion, and provide a uniform background of positive charge. Let the 
densities in the electron streams be \(1+5,) and \(1+5s,), respectively, 
and let the velocities be U+u, and —U+uy; 8), s3, u, and wy are here 
regarded as functions of position and time. We shall confine ourselves to 
one-dimensional disturbances, so that the position enters through the 
v-coordinate only. The resulting electrostatic field will be parallel to Ox; 
let E be its magnitude. 

Poisson’s equation then gives 


cE r 
a (3) 


ox 
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and the linearized equations of motion for the two electron streams are 








ou Ou et 

— +U— =—_E, (4) 
ot Cx m 

OUs Oo  € 5 
oo Oe | gon es (5) 
ct Ox m 


The equations of continuity for the two electron streams are 


( 2 ( a . Ou 

E + ( l 2 n) = |S (1 = 51) +T N (1 + 5) a = 0, 
ra e ral . Ps OUs 
E —(U —Us = [va +5)+N +5) = 0, 





Is u 
+ sah, (6) 

t Ox Ox 

OSs ~ OSs Us 
+a + oe «8 (7) 

ct Ox Ox 

The combination of (3) and (4) leads to 
0 0 Ou 

~=t+ U = bt = ()2(s5, + So), 8 
( t - -) Ox ( 1 2) ( ) 


where (92? = 47Ne?/m,. With the aid of (6), 
(5 +U =) 4 = — 0%(s, +59). (9) 


cama d ( asl =) $= —2(5, +59). (10) 


Now look for solutions of (9) and (10) in the form 
5 =a¢°-™, Ramen, (11) 
If a solution can be found, with a given real k, for which w = » +i0, where 
uw and o are real, and o > 0, then the counterstreaming is unstable, for 
the amplitude of the corresponding oscillation can grow indefinitely. The 
maximum value kmax for which this is possible gives the minimum value 
Amin of the wavelengths at which there is instability, through the relation 
Amin = 27/kmax. If such counterstreaming were ever to occur it would 
be stopped within a distance of the order of Amin. 
The evaluation of Rmax is straightforward. Substitution from (11) 
into (9) leads to 
(w+kU)Pa, = 27(a, +44), 
or [(w+RU)P— Q? Ja, — Qa, = 0. (12) 
Similarly, from (10), 
~— Qa, + [(w-RU)P- Q?]a, = 0. (13) 
Elimination of a, and a, from (12) and (13) gives 
[(w+ RUPP — Q7)[(w — RU)? — Q?] — 24 = 0, 
or wt — 2w*(R?U? + OQ?) + R2U%(R2U? — 20?) = 0. (14) 
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Since (R27U02+Q 
for all real values of k, U and 
for w?. However, if kR?U2—2 


2)2 > k?U%(kR?U2 — 20?) 

Q, it follows that (14) always has real roots 
(2? is negative, one possible value of w? is 
negative, say —o?, and this leads to 


with instability as a consequence. 

The critical wave-number is given by k?U?—2Q? = 0, so that 

| ae = V2Q/U. 
The corresponding critical wavelength is 
Amin = 277/Rmax = v27U/Q. 

Any disturbance with a wavelength exceeding Amin will be amplified. 
The counterstreaming thus becomes unstable within a distance of the 
order of Amin. In terms of N, €, m, and U, 


am, U* 
Amin = JJ Ther) (15) 


This is much shorter than the limit set to the counterstreaming by close 


2m;U* | am, U? 
aNeé © (2Ne /? 


m ee ; 8 3 U6 
that is, if Wat cae (16) 


collisions if 


To give a numerical value, when U = 10° cm/s the condition states 
that need only be much smaller than 2 x 10? particles per cm’. This 
critical density is very large. ‘The instability effect is therefore usually 
the one that arrests the counterstreaming of the electrons. 

It is harder to say exactly what happens to the protons. ‘They have a 
larger value of Amin, but in their case the greater effectiveness of the 
instability, as opposed to that of close collisions, is even more marked 
(as is shown by substitution of a larger value for m, in (15)). However 
the presence of the oscillating electron gas cannot be left out of account 
in any calculation, and the theory is rather harder. It seems nonetheless 
clear that the protons also contribute their kinetic energy to the collective 
oscillation, rather than to thermal motion. 


3. "THE EXTENSION TO THE RELATIVISTIC CASE 
Let U+u, be the velocity of an electron in the first stream. Its 
momentum is 





7 m,(U +u,) - my U My Uy U2/c? 
Pi 1 ( [ +u,)) Ce M12 (1—U? @)! 2 (1-— U? 2 1—_Uzaj* 
my l My Uy 
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to the first order. The linearized equation of motion is 


0 sa 3 
(;; sis x)bs ~e 
or 3 af ; 
’ (5 + U =) = e£/M, 
it Ox 


where M = m,(1— U?*/c?) **. ‘There is a similar change in the equation 
for the second stream. ‘Thus, the only alteration is that the rest mass mp 
is replaced by a virtual mass M = m,(1— U?/c?) 3°. 

The linearized equations of continuity do not change. 

Finally there is no change in the equation relating s, and s, to F, unless 
w/k=c. ‘To see this we note that the vector potential is parallel to Ox. 
Let A be its magnitude. ‘Then 

a (4a Ne(1+5,)(U +u,)+47Ne(1+5,)(—U+u,)}. (17) 
— Sy 1 7 r Sy + Ug) 5. 


=~ ~,9 


> > 
ox or 4 





For the scalar potential ¢, 








— = —4nNe(s, +59). (18) 


We operate on (17) with —¢/cct and on (18) with —o/ex. ‘Then since, 


134 a¢ 


we find, on adding and linearizing, that 





1 2 
Ox? c* ot* c? Ox c? ot Ox 


CE 1 AE a @ : Os 1 o Os 
E = Sa) age toNef a (a + l S,)+ salle! milter ee (uy — Us) +4 —- | (19) 


But, according to the linearized equations of continuity (6) and (7), 


(oh 1 o* o? 1 0? [ 
ee ee eee E — 47Ne a. a $e Sa dx. 20 
(3-3 a) si aap) | (its) (20) 
Equation (20) is equivalent to equation (3), unless the wave motion is such 
that 


and fe , 
\3 3 5 


that is, unless the disturbance propagates with the speed of light. 
Since the only change made in the relativistic treatment is to replace my 
by m,(1— U?/c?) 3? we find that the critical wavelength now becomes 


am, U* U?\-34 
Amin = Je )( = =) . (21) 
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4. ‘THE AMPLIFICATION OF OSCILLATIONS 

We now return to the non-relativistic case. Suppose that two uniform, 
ionized, semi-infinite streams collide with one another. Once again let 
there be N protons and JN electrons per unit volume, and let + U and —U 
be the undisturbed velocities in the two streams. Let the first impact occur 
at x = O'at time t = 0: 

We assume further that the temperature in each stream is initially zero, 
but that there are present in it random electrostatic oscillations, whose 
energy per unit mass is small compared with }U?. We shall try to find out 
how long it is likely to be before these oscillations are amplified sufficiently 
to make the linear approximation break down. When this happens the 
electrostatic energy will be comparable with the energy of counterstreaming, 
and the latter may then be expected to be brought toa halt. The calculation 
should give a closer estimate of the possible distance of interpenetration 
of the electrons in the two streams. 

The neglect of the initial temperature 7, is justified provided 
U? > kT,/m,. For, as Bohm and Gross have shown (see, for instance, 
Spitzer (1956)), the thermal motions do not appreciably attect plasma 
oscillations whose wavelengths are much in excess of y (k7\,/ Ne?), whereas 
the oscillations that are amplified most readily in the present case are those 
with wavelengths of the order of \/(m, U?/Ne?). 

The introduction of dimensionless variables will simplify the further 
working. We set 


T= 


r= 042, A = Gel VW, 6 = €E/m, QU, v, = u,/U, v, = u,/U. 


(‘The meanings of E, Q, €, m,, u, and uw, are the same as in§ 2.) ‘he equations 


¢ ‘ o 1 a & ( 7 0 i | 
(st ax} 4 (s- ax) = 4 
o o\ Ov - 0 0 CVy | 
(5 + ax) + ay = Q, (~ — “y)it OX 0, f (22) 


aby = $1; + So, 


of §2 become 


where s, and s, also have the same meaning as in §2. The equations (22) 
combine to give 


| (23) 


or 





a a \87 23 a 2 a a \2 - o \2 
(= * x) (= - x) ‘12 ba (= - x) S12 ” (= ai sx) as a 0 (24) 


(the symbol s,. stands for either s, or 59). 
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A set of coordinates alternative to X and 7 is given by 


é=7-X, n= Tt+X, 
and leads to 


( ( ( ( ¢ tf 
——-~— =2—, —+=— =2—-. 
CF ( ix ( é OT aa X ( 7) 
Equations (23) then become 
0s, 
4a 1 +5, = —5_ (25) 
( a 
Os, 
+ AED Sy = —Sj, (26) 


and (24) becomes 





Now the first stream is moving in the direction of increasing X; at 
time 7 its front surface is at X = 7, while the front surface of the second 
stream is at -7. Counterstreaming therefore occurs only in the 
part of the half-plane + > 0 which lies between the lines X =7 and 


X —r, or, alternatively, between € = 0 and 7 = 0 (region A in figure 1). 


A 
} 
STREAM |2 ONLY 
an 


®) 


(ep). COUNTER 
STREAMING 


ay) id 
» is 


i 


STREAM |! ONLY 
| « 
| 


Figure 1. The nature of the motion in various regions. 


Regions Band D in figure 1 are occupied only by streams 2 and 1, respectively. 
Region C isempty. Our equations therefore apply only in region A ; to find 
the appropriate solutions adequate boundary conditions must be known 
on the positive parts of O€ and On. These boundary conditions are 
determined by the oscillations present in the single streams occupying 
regions B and D. 

In figure 2 the diagram is re-drawn with the €-axis horizontal and the 


y-axis vertical. Suppose we seek the solution of the equations within a 
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region to the right of and above a curve PO whose gradient never takes a 
finite positive value. ‘The solution of equation (25) may be written 


. OS p : , 
5,(€, 7) = 5,(€, no)cos $4 +2 — (E, yo)sindn—2 | 5(€,u)sin }(y-—u) du, (28) 


o7) 


7 


and this shows that the value of s, at a typical point R is determined solely 
by the values of s, and 0s,/on at the point Sy, directly below R on PQ, and 
by the value of s, at all points of RS). 

Repetition of the argument shows that the value of s, at R’, say, depends 
only on the values of s, and ¢s,/e& at 7), which lies on PO to the left of R’, 
and on all the values taken by s, on TR’. Continuing in this way we find 
that the values of s, and s, at R are determined completely by the values 
of s, and 0s,/0n, sy and os,/e€ on the stretch 7, S, of PO, where Ty lies to 
the left of R and on PO. In fact these values define the function within 
the area RS, 7). It follows that the solution within the region 0 < € < &), 
0 < » < 7 1s fully determined when s, and 0s,/éy are known on the 
stretch of the €-axis from 0 to €), and when s, and os,/0€ are known on the 


7-axis, from 0 to mp. 





* 
4 
} 
| 
x To | R . 
* | 4 
\ Ve 
* if 
To\s Ry 
\ | / 
\ 
\ } J 
\ Z 
\ 4 
\ 
—— - , 7 — ¥ 
pe 
Figure 2. The domain of dependence for a point R. 


But the y-axis has the equation \ = 7 in the other set of coordinates. 
Its positive part therefore represents the front surface of stream 1 from 
time 7 =0 onwards. ‘The required boundary conditions are therefore 
equivalent to a knowledge of s, and ds,/e€ at this front surface, and of s, 
and 0s,/Cy at the other front surface, that is, to a knowledge of the conditions 
at various times in the part of stream 2 which is just about to collide with 
stream 1, and vice versa. 

The equations are linear and so a general solution may be built up by 
the addition of various basic solutions corresponding to simple boundary 
conditions. Finally the solution for the first quadrant of the (€,7)-plane 
due to two semi-infinite streams will be the same as that due to two infinite 
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streams of the same density provided only the same boundary values on 
the positive halves of the €- and 7-axes apply in the two cases. 

We shall now find the solution that corresponds to thé following simple 
boundary conditions : 

5 = 3(E-&), 35/2 = 0, 
on the positive half of the é-axis; 
So = 0, OSy do = 0) 

on the positive half of the y-axis. Here 6()stands for the Dirac delta function. 

Let P be the point (&,0). ‘The solutions we shall find will be valid 
only to the right of and above P (see figure 3). In this region, the solution 
with the above boundary conditions is the same as that given by the 
boundary conditions s, = d(€—& ), 0s,/en = 0, s, = 0, es,/eE = 0, on the 
line AB through P, provided the slope of this line is negative, say —1. 
This alternative set of conditions leads to the same values on Pé and Pn’, 
since the values of s, and 0s,/0€ are zero at all points to the right of AP 
and to the left of Py’, and therefore by continuity on Py’. Similarly the 
values of s, and 0s,/dy are zero at all points above PB and below P€, with 
the exception of the point P. 


’ 


1 An 


Wy 
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Figure 3. Region of validity of the basic solution. 


The nature of the solution corresponding to these boundary conditions 
is not affected by the exact value of &; the solutions for the points & and &;, 
are identical except for a relative displacement é,—€, parallel to Of We 
therefore choose € = 0 for the sake of simplicity. The line AB then has 
the equation +7 = 0, or t = 0 in (X,7) coordinates. 

The equivalent boundary conditions are now 


; 4 es, 1 es, os, 
ee a) te ae ee 
Sy 3( XY), On (= *t) 


0 Qs, _ 1fes, ez) _ 9 
ial a = 23. ~ ax) = | 
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or 


— | ed 


or 


With the aid of the first of the equations (23) the third and fourth 


conditions (29) may be written 


>on 7 = 0. (29a) 


A possible solution of the differential equation (27) for s, 1s 
§, = ell KX uke] 
provided y(A) satisfies 
(u? — K%)?—2(u + K4) = 0, 
that is, pt —2u?(K2+ 1)+(K4-2K?) = 0. (30) 
(This equation is the dimensionless equivalent of (14).) Corresponding 
to any given positive value of AK there are four values of x, of which two 


may be imaginary. ‘The most general solution of (27) is, in the complex 
notation, 
Sy > | a,(K)ett** AK] dK. (31) 
K 


The four z, coefficients for each K have now to be determined from the 
boundary conditions for 7 = 0. In the complex notation, 


8(X) = =) e'KX dK = s,(X,0) 


0 1 x 4 
So that sy lk, 0) “= = | Ke'4 A dK. 
Hence, from (29) and (31), 
tor , 7 l — 
> | «,(K)e** dK = -| e# dK. 


Comparison of coefficients then shows that 


a ; l 

> a(K) = wa 
Setting B. = am,, (32) 
we have > B, = 1. (33) 


The other conditions give that 


> Bu, = —K, (34) 
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> B{K+u,) = 1, (35) 
r= 
pa br B,[K + p,]° = -K. (36) 
Now the operation (35) —2K x (34) — K? x (33) gives 
4 
> Bw; = K? +1, (35a) 


r=1 


and the operation (36) —2K x (35 a) — K? x (34) gives 
4 
> Bp? = —K°-3K. (36a) 
r=1 


But the roots of (30) occur in pairs, such that 
Pic ees sy = ~My wit ps = 2(K*+ 1). 
After some manipulation the solution for 8, is found to be 
1 K (p2—K?+1 


A 4 daa KF-1f" i 


We are especially interested in the disturbances that are amplified 
most rapidly, that is, in those for which p, takes its numerically largest 
imaginary value, which is readily seen to be + 42, and occurs when K = }\3. 
In this case 8, = }. It is seen that 1, = — 47 is the appropriate coefficient 
for an amplified wave. 

In general, when p, is imaginary, the second term on the right-hand side 
of (37) gives the imaginary part of f,. This is of the first order in 
«x = K—}v3, and so |B| = 4, to the second order in « 


The solution of (30) in terms of x is 
BT = ($v3 +x)? +1—-— vV{4(3v3 +«)*4+ 1} 
= —}(1—3x?+...) 
to the second order in x, and so 
by = —4i(1— 3x? +...) (38) 


also to the second order. 
Finally, the imaginary part of f, is given by 


IV3 +f —4—-({V3 +e)? +1 
AB) = FHP} 


= — je 


| 


to the first order in «, so that 


31. 
and a, = -(!- =) (39) 








612 F. D. Kahn 


For large values of 7 the solution therefore asym~totically tends to 


= | ay(K)et om dK 


k 

L piv SA/2 | (1 = Six ph(L--3K*/2)7 91K 

aoe my | 7 IE e™ dk 
iv 3N/2,7/2 Ker) 4 DiKA 

~ r Palace e | é e'““ dk 

me age HIE gtitg Ht (40) 
\/ (1277) ; 

The solution is valid only at points where .Y ~< 7, that is at all points 


where there is double streaming. 
If the condensation s,(X,7) is produced by an input function 6(€ — 2u) 
on the line 7 = 0, the corresponding point P in figure 3 has the coordinates 





€ = 2u, » = 0, or X —u,7=u. The expression for the condensation 
is Now 
. l IND 7—-u (X+u)* 
5,(X, 7) ~~ €X E (X+u)+ ——- - —— ]. (41) 
Jlleat) TL 2 2 37 
In general, let the input function be e“”) at the point X = —7, at time 


t= 7,, for 7, > 0; this instant is on the line X+7-=7=0. ‘The corre- 
sponding condensation will be 


~ iqry, 1N3 ._. T—T; (X + | 
—_~___________ ex — x +7,)+ - dr ye 
J mag V [122(7 —ex| 2 ( 0) 2 3(7 —7,) wit 
(42) 
The integrand in (42) has the form R(7,)exp[P(7,) +¢O(7,)], where P, 


O and R are real, and P(r,) is a decreasing function of 7,. ‘The integral 
may therefore be approximated by 


sf{A,7) = 








R(O)exp[ P(0)+1Q(0)] | exp[P’(0) +70'(0)]z, dz, 


.V3X) (; =) 
, exp43 ——— Cap 4 tla — aa 


V(l2rr) v3. 1 + 
“3 SF" 3-43, * SP 


This is large only when 7 >1, and in that case the greatest values 
of s; occur where X < 7, giving, at such points, 


and so 





(X,7) = 








(.v3X ae x" 
exp;? —— exp (5 7 5a) [ 
. ° T = —_—_—_ : - : - — ‘a 
54(X, ) 5 \ (1277) i(}v3+q)-3 


When the input function is given by 


fo) = | g(get dg (44) 


— & 


we tind 





i 
] 
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(.V3X 7 X?2) 











; re | 
| 8 45) 
[= i(g+$n3) 

As we shall show, the plasma oscillation in stream 1 before the collision 
generates an input function which may be described in the form of a Fourier 
integral, as in (44). The function g(q) is such that 

(9)8*(9')> = [eo(q)e0(9')]"* 8(q — 9’) (46) 
where 8(q) is once again the Dirac delta function, and the diamond-shaped 
brackets ¢...)stand for “‘ expectation value of ...’’. There is no phase relation 
between any g(q) and any g(q’), unless gq = q’. It follows that 


exp(t — 2X*/37) /’ <s(q)g*(q') > dqdq’ 








(as )= Tie I. | P=ia+ WIG Hg + DD] 
_ exp(t—2X?/37) ¢” —_ w(q) dq ” 
en V2 ee | 1+qv3+q"° (4) 


It only remains to express this result in terms of the internal motion of 
stream 1 before the collision. ‘This is best described by means of a position 
coordinate é at rest in the stream, and the time coordinate 7. ‘he linearized, 
dimensionless equations of the oscillation are 


dvjdr = E (motion), (48 a) 
ds/er + Ov/d& = 0 (continuity), (48 b) 
a6 /0& = s (Poisson). (48 c) 
Hence 6 lor? +é = 0, 
and 6 = G(éje” (49) 


in the complex notation. ‘The form of the solution shows the well-known 
fact that the oscillations in the single stream have a frequency independent 
of the wave number, and so have zero group velocity. ‘They cannot therefore 
propagate energy, and the part of stream 1 that has not yet collided remains 
unaffected by the part that has done so. 
From (48a) and (49) it follows that 
v= iG(é)e all (50) 
and from (48c) and (49) that 
s= G'(E)e aie 
The average electrostatic energy per unit volume is 
I 1 m?2U? 


— EE* = 7 —y— 66* = {Nm U°GG*(E), 


and the average kinetic energy per unit volume of the electron gas relative 


— 
vi 
—_ 

— 





to the stream as a whole is 
1Nm,uu* = }Nm, U?vv* = | Nm u?GG*(€). 
The energy per unit volume of the electron plasma is therefore 
! Nm, U?GG*(E). 
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The input function for the motion in the region of counterstreaming 
is now G'(é)e~’, evaluated on the line X(= 7-—&) = —7, or € = 27, and 
is equal to G’(2r)e-” = f(r). 

But the field in the incident stream may be represented by 


G(éje-* = e-** | T(1) ed: (52) 


if it is a random field there are no phase relations between any [(/) and 
I'(/’), unless / = I’. We have once again 


CPIDENL)) = WDII)P2A(I=1), (53) 
and (GGxéy= | | (ore(’yydiar 
={ Wd. (54) 


The average internal energy per unit volume present in the dimensionless 
wave number range (/,/+d/) is } Nm, U?W(1)d/; the total internal energy 
density of the plasma is 





Tp = 1Nm, U? ; W(I) dl. (55) 
Further Gi(tye* = ie~* (P(e dl 
and:60 fir) = ie~* | (he dl 
t | (tg+$)P (q+ 3) ce dq. (56) 


We can now identify 9(q) with 7(1¢g+1)({q+ 1). Since there are no 
Y & 1ar4 2 T » 

phase relations among the I’s, there are none among the g’s either, as 

required. It follows from (47) that 








oT exp(r—2X?2/37) 1; ; (¢+1)* 
(81 81(4,7)) a? a ee eee 
< (P(4q-+ SE (4g! + 4) dada’ 
exp(7 — 2.X?/3r) Mq+ip , ver 
192z7 | l+qv3+¢ Wie i) 


exp(7 —2X?/3r){ ” ‘ =a 
— Ptr ENE” yal +0(1/) 





; 57 

48rr (7) 
where / is the effective dimensionless wave number range of the internal 
oscillations. Now }Nm, U? = Txin is the initial kinetic energy of motion 


of the electrons in stream 1 relative to the (x, t) system of coordinates (which 
is symmetrical with respect to the two streams). It follows from (55) 
and (57) that 

exp(7—2X?/3r) Tp 





(s,5;(X,7)) = 48a7 Tran’ (58) 
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It was one of the initial assumptions in setting up the linearized equations 
of motion that the condensation s, is much less than unity. Equation (58) 
shows that this is most probably no longer true when 
pin iil hal Gem (59) 
48ar7 = Ty : 7 
However, when s, is of order unity, the energy per unit mass of the plasma 
oscillations in the double stream will be comparable with the initial relative 
kinetic energy per unit mass of the two streams; since the former must 
grow at the expense of the latter it follows that the relative motion will by 
that time have been stopped toa large extent. ‘The minimum value z,,;, which 
satisfies (59) therefore leads to an estimate of the length of time for which 
the two streams can interpenetrate, and of the depth to which they can do so. 
The smallest value of 7 for which (59) holds occurs at X = 0, that is 
at the position where the two streams made first contact, and here 


ExP(Terit) _ Ficin 
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or Torit = = (a+loga-—1), (60) 
where a = log{48z(7Txin/Tp1)}. (61) 


To take some examples, when 7'p)/7Tkin = 10-1, 10-* or 10-°, 7,,;, = 9°6, 
14-6 or 19-6, respectively. ‘The times corresponding to these values are 


given by — my \- 
crit. ~~ QO ~ ~erit 4a Ne2 ’ 


since 2U is the initial relative velocity of the streams the corresponding 
distances of penetration of stream 1 into stream 2 are given by 


: m, U? 
2L torit => 2r, rit J (Fxa) ° (62) 


Any plasma oscillation initially present in stream 2 will lead to a further 
contribution to the right-hand side of (58). If the oscillations in the two 
streams have roughly equal energies, the effect wiil be to multiply that side 
by 2, and this leads to a decrease by log,2 = 0-69 in the estimate for 7,,;+. 

‘The more accurate estimate for the distance of interpenetration of the 
streams, given in (62), only differs by the numerical factor 1/(2/7)r,,i1; 
from the earlier, cruder estimate given in (15). ‘The numerical examples 
that have been quoted show that in any particular case 7,,;, depends, but 
only rather insensitively, on the ratio of the energies of the plasma oscillations 
present in the streams before collision to the energy of their relative motion. 
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REVIEWS 


Surveys in Mechanics: G. I. Taylor 70th Anniversary Volume, 
edited by G. K. BATcHELoR and R. M. Davies. Cambridge University 
Press, 1956. 475 pp. 50s. or $9.50. 

It would be difficult to imagine a happier birthday present than a 
dedication volume of papers. The notion is not new, but few such volumes 
can have been so personally suited as this one is to Sir Geoffrey Taylor. 
Far from being an assortment of research progress reports by fellow scientists 
wishing to pay their respects, this is a limited collection of surveys of 
‘fields of knowledge which he has made his own”’, 

Following a biographical note by Sir Richard Southwell, the contents 
are as follows. 

‘The mechanics of quasi-static plastic deformation in metals, by R. Hill 


Dislocations in crystalline solids, by N. F. Mott 
Stress waves in solids, by R. M. Davies 
Rotating fluids, by H. B. Squire 
‘The mechanics of drops and bubbles, by W. R. Lane and H. L. Green 
Wave generation by wind, by F. Ursell 
Viscosity effects in sound waves of finite amplitude, by M. J. Lighthill 
Turbulent diffusion, by G. K. Batchelor and A. A. ‘Townsend 
Atmospheric turbulence, by ‘IT’. H. Ellison 
‘The mechanics of sailing ships and yachts, by K. S. M. Davidson 


Readers of this Journal will not need to be told that this set of topics 
falls short of those to which G. I. ‘Taylor has made central contributions. 
Isotropic turbulence dynamics (in fact, the very notion of isotropic tur- 
bulence), turbulent shear flow, propulsion at very small Reynolds numbers, 
flow past porous boundaries—these too could have been included. Like 
other fluid dynamicists, I have more than once had the experience of finding 
a pet idea stolen by G. I. Taylor decades before | even thought of it! 

The singular talents of G. I. Taylor as a researcher need no more 
testimony than his own published—and unpublished—works. In an age 
of increasing dichotomy between theoretical and experimental research, 
a teacher can point to no better illustration of the benefits to be gained 
from combining the two. 

‘The papers in this volume are of good quality and some are outstanding, 
though the group is variable in spirit and in purpose. My own limitations 
preclude pertinent comment on the solid mechanics work; perhaps a 
pardonable omission in this review. No single survey here is completely 
reportorial or completely critical, although the Lighthill, Ursell and 
Batchelor—Townsend contributions come closest to the latter. 

Before taking up the individual fluid mechanics papers, I can’t resist 
observing that Sir Geofttrey’s direct descent from the mathematician Boole 
(mentioned in the pleasant biographical sketch by Sir Richard Southwell) 
in a sense runs counter to the old adage that scientific talent is handed down 


from father to son-in-law. 
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Rotating fluids 

This is a clear, well-organized summary of work on disturbances in 
fluids which are basically in rigid body rotation. ‘The chief headings are 
“1. Small disturbances’, “2. A point source on the axis of rotation” 
and ‘‘ 3. Steady symmetric motions.” 

Part 2, apparently an improved version of a problem of Barua (Quart. 
J. Mech. Appl. Math. 8, 1955), still has an ad hoc flavour. No criticism 
of the original paper is offered. Part 3, which includes the non-linear 
aspects, may treat them a bit more briefly than some specialists might 
prefer, but it does include some stimulating discussion. 

Perhaps in the interest of brevity, the title seems more general than 
nage: 6 for example, there is much additional work on rotating fluids 


in the meteorological literature. 


The mechanics of drops and bubbles 

The topics treated here are: a particular kind of pendant drop, methods 
of producing drops, the behaviour of falling drops and their break-up, 
deposition of drops on obstacles, production of bubbles, and the dynamical 
behaviour of bubbles under various conditions. 

This is a good job of assembling and organizing some remarkably 
dispersed groups of papers. As many readers will realize, the research 
literature in some of these areas is not only spread throughout a great 
assortment of professional journals, but is also presented from divergent 
points of view. 

Since much of the work has been done with practical motivations, the 
publications often emphasize semi-empirical ‘ correlations’ among para- 
meters rather than a deductive approach based on fundamental fluid 
dynamics. In a sense these authors have continued the tradition by 
restricting themselves to reporting past work without pointing out the 
places where adverse criticism might be appropriate. ‘lypical of this 
approach is the presentation of small Reynolds number droplet dynamics 
in impingement problems. Here the equations of motion are simply set 
down with the quasi-steady Stokes drag for a sphere already inserted. ‘There 
is no hint of the considerable machinery of restriction and approximation 
which must be invoked. 


Wave generation by wind 

This is a gratifyingly critical article on a branch of hydrodynamics 
which seems on the verge of blossoming. The status of most of the older 
theoretical works is summarized in the following sentence. ‘‘ Some 
hypotheses have thus remained without proof or disproof for many years, 
and through lapse of time are now taken more seriously than their authors 
would have wished.’”’ Quite apart from its merit as a piece of diplomacy, 
this sentence refers to a phenomenon noticeable elsewhere. Many engineers 
of the last 20 or 30 years have had far more literal faith in the ‘mixing 
length’ than did ‘Taylor or Prandtl. 

The general topics covered are as follows. Simple theory of stability 
of an interface, semi-empirical theories of wave generation, semi-empirical 
systems of wave forecasting, sea roughness, and surface slope due to wind. 


F.M, 27 
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There is an occasional lack of explicitness in describing experimental 
arrangements, especially noticeable where the traditionally aerodynamic 
concept of ‘boundary layer’ is involved. 

The bibliography seems reasonable ; perhapsthe best known paper omitted 
is the extensive treatise of Hellstrém (Ing. Vegen. Akad. Proc. 158, 1941). 

Perhaps the most satisfying aspect of this paper is the emphasis on 
work that should be undertaken. It is not surprising to learn that well- 
defined laboratory experiments are given the highest priority. Research 
appears to be at the stage of evolution at which simplified laboratory studies 
aimed at an understanding of the basic mechanisms of phenomena are often 
dismissed as ‘unrealistic’ by professional oceanographers. ‘This is pre- 
sumably characteristic of any activity passing from a descripto-empirical 
stage to a scientific stage. In schools of oceanography in the United States, 
for example, it seems generally expected that even theoreticians should get 
their feet wet. How long is it since theoretical aerodynamicists were expected 
to hold pilot licences ? 


Atmospheric turbulence 

Closely related to the problem of wind generated water waves is that 
of atmospheric turbulence. ‘The general topics covered in this account 
are: mean profiles of wind, temperature and water vapour near the ground, 
experimental work on mean profiles, and turbulent fluctuations. 

This is a generally good presentation, and is oriented more from 
meteorology than from fluid dynamics. ‘lhe comparison of independent 
sets of data gives a critical aspect to the account. On the other hand, the 
invocation of traditional semi-empirical formulae is largely uncritical. 
lor example, in discussion of the diffusion equation often used for heat 
and water vapour transport studies, Dr Ellison choses to omit mention 
of the apparent invalidity of the simple gradient transport concept here. 
A necessary condition for its applicability is that the characteristic 
convective scale (mean free path in kinetic theory, Lagrangian length 
scale here) be sufficiently small relative to the mean concentration field, 
and it is well known that this condition is not satisfied by conventional 
turbulent flows. Yet gradient transport calculations can often ‘ correlate ’ 
empirical data quite well, so we might say that this approach is not useful 
in principle, but merely in fact. 

Although meteorologists seem more appreciative of laboratory experi- 
ments than do oceanographers, limitations of space apparently precluded 
coverage of the few but pertinent laboratory results on turbulent boundary 
layers over rough walls. 


The mechanics of sailing ships and yachts 

The first impression of this paper is that it has been written partly in 
some language other than English. Bewildered by the use of salty argot 
at crucial points in the exposition, I fortunately had a little-thumbed copy 
of Shanty Men and Shanty Boys by W. M. Doerflinger (New York: 
Macmillan, 1951). Having decoded the message, I found the story both 
fascinating and astonishing—fascinating because of what has been accom- 
plished and recorded in the art of designing and operating sailing ships, 
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astonishing because so little use has been made up to the present time of 
scientific progress in fluid mechanics. 

The latter fact suggests that the sailing ship is on the way out. Or 
perhaps the yacht racing set is less wealthy than is commonly believed. 
On the other hand, it is possible that the sea is a force for conservatism ; 
there is no great evidence that design of motor driven ships has taken 
much advantage of modern fluid dynamics either. 

In reference to the technical content, it might have been more appro- 
priate to title this paper “the performance of sailing ships and yachts”’. 
The problem of sailing against the wind is very well expounded. On the 
other hand there is a lack of dynamical explanation for the speed figures 
tabulated. In general one gets the impression of an intriguing field of 
engineering in which the most exciting problems remain to be solved. 


Turbulent diffusion 

This is a general and authoritative survey of the field of turbulent 
dispersion, with particular emphasis on the authors’ own contributions. 
Since this domain is somewhat closer to my line of work than the rest of 
the volume, I take the liberty of more explicit comments*. 

This exposition has such good qualities of lucidity and logical organiza- 
tion that newcomers to the subject may get the impression of a more 
satisfactory current status for the turbulent diffusion problem than actually 
exists. A possible reason is some tendency toward de-emphasis of the 
distinctionsamong empirical facts,intuitive conjectures,and theoretical results. 

Typical illustrations are as follows. (a) The authors point out (p. 356) 
that the transition probability O, “‘if it is like other probability density 
functions relating to turbulent motion..., will not be too far from the 
normal or Gaussian form.’”’ Nothing in this text vicinity warns the reader 
that the Gaussian form of these other p.d.f.’s is a partly empirical fact and, 
furthermore, that very few kinds have ever been measured. (5) On pp. 358, 
359, the unwary neophyte might get the impression that both particle dis- 
placement probability P and transition probability O have been shown to 
satisfy classical diffusion equations. Of course, this is not a deductive result 
—as the authors point out on p. 360. (c) The reader is not warned that 
the method used to present the data on diffusion in decaying isotropic 
turbulence, although successful for existing data, is restrictive: it requires 
the Eulerian and Lagrangian lengths to be simply proportional during decay. 
Although this may be a good approximation during the ‘initial period’ of 
decay, the authors certainly would not wish to convey this as a general 
property of turbulence. (d) Reference on p. 392 to later “deductions ” 
that the average product of the three principal strain rates is negative may 
imply to the literate reader that this is a purely theoretical result. In fact, 
it is partly empirical as of the publication date of the volume. 

Another section for cautious reading is that on fluid line and surface 
stretching. ‘The presentation here is as heuristic as the pioneering paper 
of Dr Batchelor. Having encountered a number of fluid dynamicists 
with the impression that the original paper proved that fluid lines stretch, 


*S. Potter, The Art of Reviewsmanship, New York Times, 8 May, 1955. 
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I should like to observe that the authors make no such claim. The fluid 
line stretching is a conjecture (and inferentially an empirical result from 
vorticity behaviour). ‘lhe exponential rate follows from quite intuitive 
assumptions. 

The discussion of lateral diffusion in shear flow is a mixture of careful 
formalism and ad hoc simplifications. ‘The derivation of equations through 
(3.19) is an improvement over earlier published work. ‘Then the quasi- 
Gaussian postulate is invoked in a situation where the non-Gaussian 
character of the fluctuation field may well be crucial. With a bit of hand 
waving, it is likely that the consequent approximation, uv? = 3u'v.v*, 
can be qualitatively justified : 7? is positive, and 7* has the same sign 
as v. Whether it can be given any quantitative justification seems unclear. 


9 


Somewhat more startling is the bland assumption that u’v? = u’v.v3/v 
It is not evident a priori that even the sign of this approximation is correct. 

The brief discussion of the behaviour of a diffusive scalar field embedded 
in an isotropic turbulence ignores much of the existing work on the subject 
(Obukhov, /svest. Akad. Nauk, S.S.S.R., Ser. Geofiz. 13, 1949, and 
Corrsin, J. Aero. Sct. 18, 1951; F. Appl. Phys. 27, 1951). 

‘There are at least two specific improvements over previous writings 
on these subjects by Drs Batchelor and ‘Townsend. (A) Equation (3.1), 
expressing the rate of growth of a fluid surface element, is in more elegant 
form than in the original paper of Batchelor. (B) The turbulence behind 
a regular grid is no longer referred to as being isotropic. ‘The non-isotropy 
of these flows has been known to turbulence researchers in the United 
States since perhaps 1942 (Corrsin, A.E. Thesis, Cal. ‘Tech., 1942), but the 
great importance and volume of the published grid turbulence data coming 
from the Cavendish Laboratory has given a rather general impression that 
virtual isotropy prevails within 16 or 20 mesh length of the grid. The 
latest Cavendish data are apparently in better agreement with those that 
have been taken at other laboratories. 

From a pedagogical viewpoint the only serious omission is perhaps a 
statement of the nature of the problem of relating Lagrangian dispersion 
statistics to Eulerian statistical information. ‘This, after all, is a possible 
definition of ‘the problem’ of turbulent diffusion. Given the statistical 
properties of the Eulerian (space-time) velocity field u,(.v;, t), the Lagrangian 
velocity field is v(a,t) = u,(X(a, t), t). (1) 


The fluid point displacement is 


“t 
X,,(a,t) = | v,(a, 0) dt, (2) 
and a; = X,(a,0). 


Each realization of the Lagrangian time variation following a tagged 
Huid point (a; = constant) is thus the variation along a trajectory through 
a member of the ensemble of Eulerian fields. ‘The central difficulty of the 
problem is indicated by noting that each irregular ‘sampling trajectory’ 
is analytically related to the field it samples. ‘Therefore we must deal 
with a stochastic integral equation. Putting (1) in (2) gives 


ej 
X,(a,t) = | u,(X(a, ty), ty) dey. 


“0 


_ 
w 
tte 
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The Lagrangian two-point correlation tensor is 
L;, = v(a, t)v,(a+e, t+7) 
= u(X(a, t), t)hu,(X(a+a,t+7),t+7). 
For stationary fields, this depends on (@,7) only. ‘Taylor expressed the mean 
square particle displacement in terms of the Lagrangian time correlations 
L,(9, 7). 

The theory of correlation functions obtained by sampling with random 
intervals over an ensemble of functions appears to be new to probability 
theory, even for the degenerate case in which the individual intervals are 
unrelated to the individual members of the ensemble being sampled. Yet 
we face here the still more complicated case in which each member of the 
function ensemble determines its own sampling interval through the 





integral operation (3). 


Viscosity effects in sound waves of finite amplitude 

This is a condensed 100 page review which treats virtually all phases 
of the title subject. Especially noteworthy are the careful clarification and 
correction of the various traditional approximations, and the use of Burgers’ 
equation (where it is really appropriate), designed originally as a ‘ model’ 
of turbulence. Equally appealing is the author’s consistent emphasis on 
the Reynolds numbers of the various kinds of waves. ‘The principal subject 
headings are as follows. ‘lhe physical mechanism of viscosity and other 
ditfusion effects in gases, attenuation of sound due to relaxation effects, 
classical theory of shock wave formation, structure of stronger shock waves, 
more general theory of plane shock wave formation, interiors of unsteady 
shock waves, formation and decay of non-plane shock waves, sound waves 
of moderate Reynolds number, relaxation effects in sound waves of finite 
amplitudes. 

The remarkable coverage of this review will undoubtedly make it a 
standard reference for a long time to come. ‘There are perhaps two other 
facets that could be included to round out the subject. Chapter two, 
on the physical mechanism, stops short of statistical mechanics. At the 
phenomenological level of the fluid dynamicist, the underlying physical 
mechanism resides in the molecular dynamics. Professor Lighthill has 
chosen to omit discussion of the challenging and poorly understood problem 
of how the probability density of molecular velocities changes from one 
form to another in a few mean free paths through a shock wave. It might 
have been appropriate to make special mention of the measurements of 
F. S. Sherman, giving profiles through a shock wave (Nat. Adv. Comm. 
Aero., Wash., Tech. Note no. 3298, 1955). 

In conclusion, this volume is an impressive tribute to one of the great 
engineering scientists of this century. ‘The only ‘Taylor-associated writings 
I expect to enjoy more are the research papers he will publish himself in 
the future! 


S. CoRRSIN 
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